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TANKWARS  for  the  Parametric  Consideration  of  System  Concepts 

by 

Harry  L.  Reed,  Jr. 


INTRODUCTION 


The  core  of  the  effort  was  to  develop  and  exercise  a  methodology 
for  the  analysis  of  future  armament  concepts  for  armored  vehicles  to  allow 
the  assessment  of  the  application  of  new  technology  to  those  weapon 
systems.  TANKWARS  was  picked  as  a  good  tool  for  starting  this  effort. 

Thus  TANKWARS  as  adapted  and  used  to  carry  out  numerous  parametric  studies 
of  tank  armament  concepts.  The  data  from  these  studies  were  furnished  as 
developed  to  Dr.  Howe. 

The  data  developed  are  of  little  use  to  anyone  who  is  not  aware  of 
the  classified  issues  involved,  and  in  many  instances  the  data  represent 
stages  in  the  thinking  about  what  the  methodology  should  do  and  what  issues 
should  be  given  attention  for  further  analysis;  thus  Dr.  Howe  decided  that 
no  useful  purpose  would  be  served  in  including  those  data  in  this  report. 

Of  course,  the  data  are  being  used  by  Dr.  Howe  to  develop  his 
recommendations  for  the  Army's  Technical  Base  Program  in  armaments. 

OBJECTIVE 

This  report  then  will  discuss  aspects  of  TANKWARS  relevant  to  the 
above  work,  changes  that  were  made  in  the  TANKWARS  program,  ideas  for  the 
future  application  of  TANKWARS  (or  possibly  its  successor)  in  further 
parametric  analyses,  and  some  thoughts  on  the  development  of  more  advanced 
tools  for  modeling.  This  report  should  be  specially  useful  as  an  aajunct 
to  documentation  on  TANKWARS  or  on  GROUNDWARS. 

1.  TANKWARS  -  A  Computer  Simulation 

TANKWARS  is  a  Monte  Carlo  computer  simulation  of  engagements 
between  two  homogeneous  mechanized  forces.  The  model  simulates  individual 
weapon  systems, and  the  engagements  include  search,  detection,  selection, 
firing,  impact,  functional  destruction,  disengagement,  and  reengagement. 
Nominally  it  can  handle  up  to  20  armored  vehicles  on  each  side.  The 

computer  program  was  written  by  Mr.  Fred  Bunn  of  the  Ballistic  Research 

Laboratory  (BRL) ,  Aberdeen  Proving  Ground,  Maryland  (reference  1). 

GROUNDWARS  is  an  outgrowth  of  the  TANKWARS  model  and  was  developed 
by  the  Army  Material  Systems  Analysis  Activity  (AMSAA) .  Some  of  the 
features  of  GROUNDWARS  are  addressed  in  Section  3.  At  the  time  we  first 
discussed  GROUNDWARS  with  AMSAA,  version  3.97  was  running,  version  4.0  was 
being  worked  on,  and  a  new  concept  using  the  "C"  language  was  in  the  early 

stages  of  development.  We  have  had  some  very  useful  discussions  with  AMSAA 

people  (specially  Barry  Burns)  which  we  have  found  most  useful  and  which  we 
hope  were  of  some  use  to  AMSAA. 
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We  felt  that  TANKWARS  provides  a  reasonable  m  on  n  model  that 
includes  a  level  of  detail  about  the  systems  under  consideration,  that  it 
could  allow  sensible  engineering  conclusions,  and  that  it  represents  a 
tactically  meaningful  environment  for  giving  military  significance  to  any 
such  conclusions. 

An  important  consideration  was  running  time  since  the  parametric 
analyses  require  the  running  a  large  numbers  of  cases.  TANKWARS  appeared 
to  offer  a  reasonable  balance  between  system  detail  and  tactical  context 
with  practical  run  times.  The  significance  of  running  time  became  even 
more  obvious  when  the  computer  originally  picked  for  use  (a  Gould  9000) 
proved  to  be  too  slow,  and  the  work  was  moved  to  a  Cray  X/MP-48.  The 
latest  version  4.0  of  GROUNDWARS  appears  interesting  in  this  regard  as  a 
considerably  speeded  up  and  improved  successor  to  TANKWARS.  Initially, 
GROUNDWARS  3.97  appeared  to  be  more  complicated  than  necessary,  and  it  was 
too  large  to  fit  on  the  Gould  computer  (or  so  the  compiler  kept  saying). 

Of  course,  the  use  of  the  larger  Cray  and  the  recent  availability  of  the 
faster  version  of  GROUNDWARS  are  reason  for  reopening  that  decision. 

Other  advantages  of  TANKWARS  were  felt  to  be: 

*  It  is  well  enough  written  to  allow  change  -  an  important 
consideration  since  new  weapon  concepts  often  require  changes  in 
employment  for  full  realization  of  their  potential.  This  is 
specially  the  case  in  modern  fire  control  concepts  which  offer 
concurrent  advantages  in  target  acquisition  and  in  command  and 
control . 

*  TANKWARS  is  in  use  in  the  analysis  community  and  is  the  basis 
for  AMSAA's  GROUNDWARS.  Transition  to  GROUNDWARS  4.0  should  be 
relatively  easy  if  that  proves  to  be  desirable.  In  fact,  many  of 
the  subroutines  are  the  same. 

*  Mr.  Bunn,  the  developer  of  TANKWARS,  was  easily  accessible  for 
numerous  discussions  on  nuances  of  the  program  and  changes. 

One  feature  of  TANKWARS  that  was  not  utilized  in  the  effort 
described  in  this  report  was  that  of  modeling  sequential  engagements  and 
looking  at  resupply  on  one  side.  AMSAA  also  did  not  include  this  feature 
when  they  created  GROUNDWARS. 

As  was  mentioned  before,  the  general  acceptance  of  TANKWARS  is 
testimony  to  its  usefulness.  But  there  are  always  needs  for  more.  For  the 
purposes  at  hand  the  shortcomings  of  TANKWARS  include: 

Only  one  type  of  armored  vehicle  (including  one  type  of  weapon  and 
one  type  of  ammunition)  is  allowed  on  each  side.  Thus  a  mix  of  weapons 
that  might  have  application  at  different  ranges  in  the  engagement  cannot  be 
played.  Neither  can  different  systems  in  a  cooperative  arrangement  such  as 
tanks  advancing  with  other  vehicles  in  protective  overwatch  with  guided 
missile  systems. 
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The  available  scenarios  only  include  an  advance  by  one  side  or  a 
meeting  engagement  in  which  both  sides  are  stationary  throughout  the 
engagement.  We  see  the  need  for  the  ability  to  conduct  engagements  between 
two  moving  (and  occasionally  stopping)  forces. 

Likewise  there  is  a  need  for  simulating  more  complicated  overwatch 
tactics  such  as  leapfrogging.  As  will  be  noted  further  on,  a  simple 
overwatch  tactic  has  been  added  in  which  the  overwatching  force  never  moves 
and  has  the  same  weapon  as  the  advancing  force. 

The  decision  process  for  engaging,  disengaging,  and  reengaging  also 
needs  to  be  reworked.  For  example,  the  ability  might  be  added  to 
disengage  from  a  low  priority  target  if  a  higher  priority  target  were  to 
become  available.  At  present  the  disengagement  criteria  are  uniformly 
applied  based  on  time,  or  based  on  rounds  fired,  or  based  on  damage  to  the 
target,  but  are  not  related  to  the  relative  importance  of  the  target.  The 
importance  of  the  target  comes  into  play  only  in  the  initial  selection. 
Faking  disengagement  by  firing  only  a  small/fixed  number  of  rounds  and 
using  the  subsequent  target  selection  process  to  represent  the  search  for  a 
more  important  target  is  not  an  adequate  surrogate.  The  decision  to 
reengage  the  old  target  would  not  be  the  same  as  the  simple  decision  to 
continue  the  firing  sequence  since  the  time  to  fire  the  "next"  round  would 
be  the  time  to  fire  a  "first"  round  rather  than  the  shorter  time  to  fire  a 
subsequent  round.  When  runs  were  made  in  which  target  information  was 
shared  among  the  tanks  in  a  force,  the  result  was  often  negative.  This 
indicates  that  the  additional  information  was  not  being  used  properly  by 
the  decision  process  associated  with  engagement  and  disengagement. 

Even  when  done  on  the  X/MP-48,  the  calculations  for  vast  parametric 
arrays  are  somewhat  slow.  Most  of  the  computing  time  is  spent  in  the 
search  and  related  target  acquisition  process  which  is  repeated  every 
second.  The  version  4.0  of  GROUNDWARS  treats  this  differently  and  achieves 
a  significant  speed  increase. 


2.  Changes  to  TANKWARS 


In  talking  about  the  changes  made  to  TANKWARS  during  this  effort, 
we  often  refer  to  the  FORTRAN  source  program  listed  in  Appendix  A.  That 
listing  is  only  given  to  allow  one  to  get  an  idea  of  what  the  program  does 
and  more  importantly  to  show  the  changes  that  were  made.  Anyone  who  is 
interested  in  acquiring  the  code  for  use  should  contact  Mr.  Bunn  at  BRL  for 
documentation. 

The  program  as  shown  includes  a  variety  of  changes  by  the  author 
that  are  often  rather  crude  patches.  No  uniform  attempt  has  been  made  to 
clean  up  the  modified  code,  clear  out  any  parts  that  have  been  rendered 
useless  by  the  changes,  rework  the  code  to  make  the  additions  more 
efficient,  and  blend  the  changes  into  the  fabric  of  the  code  more  smoothly. 
Since  there  appears  to  be  good  reasons  to  move  to  a  more  advanced 
program,  there  would  be  no  good  purpose  served  by  such  an  effort.  Also,  a 
very  large  portion  of  the  run  time  is  spent  in  the  process  of  searching  for 
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targets  and  updating  the  positions  of  the  vehicles;  so  any  small 
inefficiencies  introduced  into  most  of  the  code  are  not  a  matter  of  great 
significance. 

The  UNIX  environment  was  used  for  the  minicomputer  operation  and 
UNICOS  for  the  operation  on  the  Cray  computer. 

As  cases  were  run,  we  found  hiccups,  bugs,  and  needs  for  changes  in 
the  program.  We  also  developed  a  sense  of  what  we  shouT-d  be  looking  at  and 
looking  for. 

The  parametrics  considered  include  accuracy,  lethality,  acquisition 
(regular  and  pinpoint),  rates  of  fire,  time  to  lay  the  weapon,  decoys, 
overwatch,  decision  algorithms,  and  fire  control  options.  Questions  have 
proliferated  more  rapidly  than  answers,  and  that  trend  will  probably 
continue. 

TANKWARS  uses  a  number  of  input  files.  The  game  file  and  the 
miscellaneous  file  are  discussed  in  considerable  detail  in  reference  1: 

*  the  game  file  sets  up  the  numbers  of  players,  the  ranges  of 
engagement,  the  types  of  scenarios  (e.g.  who  is  defender  and  who  is 
aatacker) ,  and  other  data  including  the  names  of  files  that  define 
the  characteristics  of  the  weapon  systems  being  played.  For  the 
version  of  TANKWARS  discussed  in  this  report  there  are  two  files 
for  each  of  the  two  armored  vehicles,  a  misc  file  and  a  vul  file. 

*  the  misc  file  for  a  side  gives  a  variety  of  data  about  the 
vehicle  on  that  side  such  as  its  dimensions,  speed,  ability 
to  detect  targets,  rate  of  fire,  etc. 

*  the  vul  file  for  a  side  gives  the  probabilities  of  hit  and  of 
various  types  of  kill  given  a  shot  for  the  weapon  used  by 

that  side  against  the  vehicle  on  the  other  side.  As  noted  below, 
the  original  version  of  TANKWARS  used  two  files  here,  one  for 
accuracy  and  one  for  conditional  probability  of  kill  given  a  hit. 


2.1  Changes  for  Running 

For  the  convenience  of  running  numerous  parametric  cases,  a  shell 
file  approach  was  created  for  overall  management  of  the  calculation.  A 
few  game  files  were  created  each  of  which  used  the  same  names  for  the  other 
input  files  -  namely  blue. misc,  blue. vul,  red. misc,  and  red. vul. 

Then  a  UNIX  shell  file  (called  runtw  for  example)  was  created  such 
as 


cp  $2  blue. misc 
cp  $3  blue. vul 
cp  $4  red. mi  sc 
cp  $5  red. vul 
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echo  Defender:  $2";"  $3  >  $6 
echo  Attacker:  $4";"  $5  »  $6 
. ./source/twxl  <  $1  »  $6 

This  shell  program  was  called  by  a  command  line  such  as: 

runtw  filel  file2  f i 1 e3  file4  file5  file6 

This  command  line  produces  what  we  will  call  a  case  in  this  report. 
A  case  is  a  series  of  runs  each  with  the  number  of  Monte  Carlo  replications 
specified  in  filel  (the  game  file)  and  each  with  one  of  the  opening  ranges 
specified  in  filel.  The  mi  sc  file,  file2,  describes  the  blue  force 
vehicles.  The  blue  vul  file,  file  3,  contains  data  on  the  ability  of  the 
blue  vehicles  to  kill  the  red  vehicles.  The  misc  file,  file4,  is  for  red. 
The  vul  file,  f i 1 e5  is  for  red  kills  of  blue.  The  file,  f i 1 e6 ,  is  that 
into  which  the  output  data  are  to  be  stored. 

Finally  a  number  of  cases  can  be  run  by  setting  up  a  shell  file 
containing  a  number  of  these  command  lines. 

As  a  matter  of  convenience  we  used  the  blue  force  as  the  defending 
force  and  the  red  force  as  the  attacking  force  rather  than  using  these 
terms  as  indicators  of  political  affiliation. 

TANKWARS  prints  a  considerable  amount  of  detailed  data  about  each 
run  made  at  each  opening  range.  Except  for  occasional  diagnostics,  these 
data  were  more  than  were  needed,  and  vast  outputs  were  being  produced  for 
the  large  number  of  runs.  Therefore  the  indicated  change  in  the  subroutine 
NXWAVE  was  introduced  and  used  to  produce  one  line  of  data  for  each  opening 
range  in  a  case.  As  that  change  is  shown,  that  output  line  contains: 

R  Ndef  Ndefdec  Natt  Noverw  Exch  Defrds  Attrds 

where 

R  =  the  opening  range 

Ndef  =  the  average  number  of  defenders  (not  decoys)  killed 

Ndefdec  =  the  average  number  of  defender  decoys  killed 

Natt  =  the  average  number  of  advancing  attackers  killed 

Noverw  =  the  average  number  of  overwatching  attackers  killed 

Exch  =  the  exchange  ratio  *  Ndef  /  (Natt  Noverw) 

Defrds  =  the  average  number  of  rounds  fired  by  a  defender  vehicle 

Attrds  *  the  average  number  of  rounds  fired  by  a  defender  vehicle 

Other  formats  could  be  easily  tailored  as  needed. 

2.2  Modifications  for  Decoys  and  Overwatch 

We  were  concerned  with  decoys  that  were  stationary  and  flashing 
(simulating  firing).  We  also  included  the  ability  to  have  tanks  on  the 
side  of  the  attacking  force  that  remained  at  the  opening  range  in  hull 
defilade;  these  were  called  overwatching  tanks. 
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The  changes  needed  for  decoys  and  overwatching  tanks  are 
intertwined  in  the  code  and  are  contained  in  the  subroutines  CANGO,  DEATHS, 
DEPL02,  FINISH,  INIT2,  INPUT,  and  SEARCH.  Decoys  were  already  in  TANKWARS, 
but  the  changes  were  needed  to  remove  them  from  win  statistics  (changed  to 
say  that  if  all  the  real  tanks  are  killed  on  a  side,  that  side  loses  even 
if  the  decoys  are  alive)  and  to  count  only  real  tanks  in  the  exchange 
ratios. 


The  changes  in  INIT2  include  ordering  the  real  and  decoy  tanks  so 
that  if  the  attacking  force  (called  red  in  the  program)  has  any  decoys,  the 
decoys  will  be  in  the  overwatching  role  as  long  as  there  are  no  more  decoys 
than  there  are  overwatching  tanks.  Since  we  didn't  have  a  good  idea  on 
how  to  treat  moving  decoys,  we  only  considered  stationary  ones.  The 
stationary  decoys  are  handled  much  the  way  that  tanks  are  handled.  They 
acquire  and  "shoot"  at  tanks  but  no  rounds  fly  to  the  target,  and  they  are 
acquired  and  killed  as  are  tanks. 

The  change  in  SEARCH  allows  the  defender's  acquisition  ability  to 
account  for  both  hull  defilade  stationary  overwatchers  and  fully  exposed 
moving  attackers.  We  further  decided  that  an  attack  would  be  stopped  (and 
thus  the  attacking  side  would  lose  that  case  if  all  the  advancing  tanks 
were  killed  regardless  of  whether  the  overwatching  tanks  were  killed  or 
not.  A  more  realistic  model  would  include  the  possibility  for  the 
advancing  tanks  and  the  overwatching  tanks  to  trade  places  (leapfrogging). 

2.3  Change  to  Ensure  Loading  of  "First"  Rounds 

This  became  a  problem  when  some  of  the  time  constants  associated 
with  laying  the  gun  became  less  than  the  time  needed  to  load  a  round.  As 
it  was  written,  TANKWARS  accounted  for  the  time  needed  to  reload  the  gun 
assubsequent  rounds  were  fired  at  a  target.  However,  no  such  loading  time 
was  "required"  from  the  time  the  last  round  was  fired  at  one  target  until 

the  next  round  was  fired  at  another  target.  This  was  not  a  problem  for  the 

very  first  round  since  it  was  assumed  that  the  tank  would  go  into  battle 
with  a  round  in  the  tube.  Changes  were  made  in  the  subroutines  ENGAGE, 
HALT,  and  INIT2.  The  change  in  INIT2  was  to  put  in  a  bogus  value  for 
tfire2  which  was  tmin  seconds  before  the  battle  would  begin.  The  variable 
tfire2  is  the  time  that  the  tank  last  fired  and  tmin  is  the  time  to  reload. 

Thus  the  tank  could  fire  if  needed  as  early  as  t  =  0. 

2.4  Stop  to  Fire  Changes 

Most  of  the  runs  that  have  been  made  under  this  project  were  made 
with  an  attacker  who  could  shoot  on  the  move,  and  the  attacker  rarely  fired 
from  a  stationary  position  (only  when  he  was  mobility  but  not  fire  power 
killed).  Some  excursions  were  made  into  the  question  of  the  advantages  (if 
any)  of  stopping  to  shoot.  Some  changes  were  required  in  the  subroutine 
HALT  which  accommodates  firing  when  stopped.  The  change  adds  the  need  to 
load  for  the  "first"  round  as  mentioned  above,  and  the  change  also  deletes 
the  feature  that  removes  3  seconds  from  the  laying  sequence  (which 
presumably  was  included  since  the  laying  process  might  have  been  started 


before  the  tank  had  come  to  a  halt).  However  it  is  possible  for  a  tank 
which  is  stopped  and  has  just  finished  firing  at  a  target  to  sta.  i.  to 
move  and  immediately  find  another  target  and  halt  in  less  than  3  seconds. 
This  is  exacerbated  by  the  often  short  times  we  have  been  using  for  the 
time  to  lay  the  first  rounds,  and  it  is  further  complicated  by  the  fact 
that  we  have  been  using  times  to  lay  the  gun  that  are  different  for 
stationary  firing  and  for  moving  firing  for  our  baseline  fire  control 
system.  It  is  unclear  how  to  pro  rate  these  two  laying  times  during  the 
stopping  process.  While  not  thoroughly  satisfying,  we  decided  to  let  the 
vehicle  come  to  a  halt  and  then  use  the  shorter  lay  time  associated  with  a 
stationary  firer. 

Note  that  the  piece  of  programming 
tlast  =  t  -  3. 

remains  as  a  relic  in  the  modified  subroutine. 

2.5  Changes  for  Cardioid  Distribution 

The  original  TANKWARS  handled  the  consideration  of  the  azimuthal 
orientation  of  hits  on  the  vehicles  by  creating  a  distribution  of  attack 
geometries.  The  advancing  force  proceeded  on  a  course  that  for  each  Monte 
Carlo  replication  had  an  orientation  that  was  randomly  chosen  (usually 
cardioid)  about  the  axis  that  joined  the  center  of  the  target  array  and  the 
center  of  the  attacking  force.  Thus  the  tanks  often  marched  in  a  direction 
that  resulted  in  the  forces  never  coming  very  close.  This  led  to  a  number 
of  indecisive  repl ications^within  a  sequence  of  Monte  Carlo  runs. 

AMSAA  does  not  have  this  feature  in  GROUNDWARS.  The  attackers 
advance  toward  the  target,  and  they  apply  a  random  angle  at  the  time  of 
each  impact. 

TANKWARS  and  GROUNDWARS  use  two  tables  -  a  table  of  dispersions  for 
different  cases  and  ranges  and  a  table  of  conditional  kill  probabi 1 ities 
for  different  kill  criteria,  dispersions,  ranges,  exposures,  and  angles  of 
attack.  The  probability  of  hit  is  calculated  using  normal  distributions 
for  hits  on  the  turret  and  chassis  and  then  the  conditional  kill 
probabilities  are  multiplied  by  the  hit  probability  to  obtain  the 
probabilities  of  kill  for  that  instance.  These  kill  probabilities  are 
compared  against  a  random  number  and  one  (or  none)  is  selected. 

We  have  adopted  a  scheme  that  has  the  attackers  advancing  directly 
toward  the  defenders  as  does  GROUNDWARS;  but  we  also  do  a  lot  of 
preprocessing  of  the  data  to  account  for  the  angular  distribution  of  hits. 
The  rationale  follows: 

The  conditional  kill  tables  are  large,  having  some  6000  entries. 

The  hit  probability  calculations  involve  Gaussian  function 

calculations. 
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We  had  hopes  that  we  eventually  could  add  more  than  one  type  of 
ammo  for  each  side;  so  that  the  memory  requirements  would  get 
large. 

Almost  all  the  shooting  was  either  by  stationary  vehicles  shooting 
at  moving  vehicles  or  vice  versa.  The  issue  of  the  difference 
between  the  probability  of  hit  for  first  and  for  subsequent  rounds 
in  stationary  fire  against  stationary  targets  was  moot  (this 
involves  random  bias  and  random  dispersion). 

The  approach  was  to  create  a  table  of  hit  and  kill  probabilities 
(given  a  shotjas  functions  of  range  for: 

*  the  three  cases  of: 

stationary  shooting  at  stationary  targets 
stationary  shooting  at  moving  targets 
moving  shooting  at  stationary  targets 

*  two  target  conditions  of: 

hull  defilade 
fully  exposed 

*  five  levels  of  kill 

hit 

mobility  kill 
fire  power  kill 
mobility  and  firepower  kill 
catastrophic  kill 

This  along  with  nine  range  values  (0  to  4000  meters  by  500  meters) 
gives  a  table  with  270  entries  and  puts  a  lot  of  mathematical  calculation 
outside  the  Monte  Carlo  replications.  The  only  drawback  seems  to  be 
the  need  to  do  something  about  the  correlation  of  impacts  for  the 
stationary  /  stationary  case  should  that  become  important. 

The  changes  to  accommodate  this  kill  model  are  in  the  subroutines 
DAMAGE,  INPUT,  KILL,  MAYHIT,  and  RDPKH.  Note  that  the  subroutines  ACCERR, 
ACCMS,  ACCSM,  ACCSS,  IZHIT,  and  RDER0R  are  not  needed. 

Most  of  the  changes  are  straightforward  and  relate  to  shortening 
the  calculation  by  avoiding  specific  efforts  to  calculate  hit  probability. 

One  piece  of  the  code  is  worth  discussing.  While  it  is  very 
similar  to  the  original  code,  it  must  be  understood  to  understand  the  BASIC 
program  given  below  for  the  calculation  of  the  vulnerability  table. 

It  is  in  the  subroutine  KILL: 

temp  =  ranu(C.O)  ranu  returns  a  random  number 
IF  (temp  .gt.  p ( 1 ) )  THEN 
c  no  hit  and  no  kill 

hit  =  .false. 
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injury  =  ALIVE 

ELSEIF  (temp  .gt.  p(2) )  THEN 
c  a  hit  and  a  "k"  kill 

hit  =  .true, 
injury  =  KKILL 

ELSEIF  (temp  .gt.  p(3) )  THEN 
c  a  hit  and  an  "m&f"  kill  but  no  "k" 
hit  =  .true, 
injury  =  MFKILL 

ELSEIF  (temp  .gt.  p(4))  THEN 
c  a  hit  and  an  "f"  kill  but  no  "m"  kill 
hit  =  .true, 
injury  =  FKILL 

ELSEIF  (temp  .gt.  p(5))  THEN 
c  a  hit  and  an  "m"  kill  but  no  "f"  kill 
hit  =  .true, 
injury  =  MKILL 

ELSE 

c  a  hit  but  no  kill 

hit  =  .true, 
injury  =  ALIVE 

ENDIF. 

In  this  section  of  code  the  random  number  (temp)  is  located  within 
a  collection  of  segments  in  the  unit  interval  which  represent  various 
mutually  independent  outcomes  of  the  hit  or  miss: 

p(l)  is  the  probability  of  a  hit;  so  if  1.0  >  temp  >  p(l)  the  round  missed, 


p(2)  is  the  probability  of  any  kill  less  that  a  K  (catastrophic)  kill;  so 
if  p ( I)  >  temp  >  p (2)  the  round  hit  and  achieved  a  K  kill, 

p(3)  is  the  probability  of  any  kill  that  is  not  both  a  mobility  and  a 
firepower  kill;  so  if  p (2)  >  temp  >  p (3 )  the  round  hit  and  achieved  both  a 
mobility  and  a  firepower  kill, 

p(4)  is  the  probability  of  a  mobility  kill  but  not  a  firepower  kill;  so  if 
p(3)  >  temp  >  p(4)  the  round  hit  and  achieved  a  firepower  kill  but  not  a 
mobility  kill, 

p(5)  is  the  probability  that  the  hit  produced  no  kill;  so  if  p(4)  >  temp  > 
p(5)  the  round  hit  but  did  not  kill. 

With  some  thought  the  reader  should  be  able  to  convince  himself  that: 

1.0  >=  p(l)  >=  p(2)  >=  p(3)  >=  p(4)  >*  p(5)  >=  0.0. 

A  BASIC  program  has  been  written  to  produce  the  vulnerability  table 
from  accuracy  data  and  from  the  BRL  vulnerability  data  format.  This 
program  considers  the  turret  and  the  hull  to  be  two  shoeboxes.  The 
particular  code  shown  in  Appendix  B  treats  bcth  boxes  with  the  same 
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cardioid  distribution  for  all  cases  (fully  exposed,  hull  defilade,  moving, 
stationary).  This  can  be  easily  expanded  to  consider  different 
distributions  for  different  cases,  and  could  be  extended  to  different 
distributions  for  the  hull  and  the  turret  for  the  same  case.  The  latter 
would  require  some  rework  of  the  original  vulnerability  data  which  only 
treat  the  hull  and  the  turret  together  for  one  aim  point  and  the  turret 
separately  for  another  aim  point.  We  would  need  the  data  for  the  hull  and 
the  turret  separately  for  the  one  aim  point  associated  with  the  fully 
exposed  target.  Such  data  are  basically  available,  but  probably  not  in  the 
right  form. 

The  BASIC  program  was  written  for  and  run  on  a  PC.  It  could  of 
course  be  easily  written  in  FORTRAN  or  left  in  BASIC  and  run  on  a  more 
capable  computer  if  desired.  The  program  takes  about  10  minutes  to  run 
under  interpreted  BASIC  on  a  Tandy  1000SX.  This  is  no  problem  since  it 
need  run  very  infrequently. 

2.6  Detection  and  Pinpoint  Changes 

Changes  were  made  in  the  two  types  of  detection  addressed  by 
TANKWARS  -  pinpoint  detection  and  the  detection  of  targets  without  using 
their  firing  signatures. 

Here,  pinpoint  detection  is  defined  as  the  detection  of  the  firing 
of  an  opponents  gun  and  the  subsequent  location  of  the  firing  vehicle  well 
enough  for  the  person  detecting  to  be  able  to  aim  his  gun  effectively  at 
the  firer.  Two  probabilistic  concepts  are  involved  -  the  probably  that  the 
observer  can  achieve  a  pinpoint  detection  as  defined  above  and  the  median 
time  it  takes  him  to  accomplishment  that  pinpoint  detection. 

Pinpoint  is  handled  in  the  subroutine  PINPNT.  One  change  was  to 
allow  the  median  time  for  the  pinpoint  process  to  be  input  from  the 
appropriate  misc  input  file.  TANKWARS  has  this  value  as  a  fixed  number  and 
a  recompilation  of  the  program  was  required  each  time  the  value  was 
changed.  There  is  a  related  change  in  the  subroutine  RDMJSC  that  gets  this 
median  time  from  the  misc  file. 

We  found  (as  have  others)  that  pinpoint  was  an  important  capabi  lity 
for  the  attacker  who  is  looking  for  otherwise  often  hard  to  find  targets 
that  are  in  hull  defilade.  One  thought  for  improved  systems  was  to  add  the 
capability  for  a  tank  to  tell  other  tanks  on  his  side  that  he  saw  a  firing 
and  give  them  the  location  of  that  targets.  To  test  this  idea,  we  added  a 
change  to  the  subroutine  PINPNT  that  allowed  the  other  tanks  to  locate  the 
target  at  the  same  time  that  the  observer  saw  it.  So  far  we  have  been 
working  with  this  optimistic  version  of  the  idea.  This  could  be  easily 
changed  in  the  PINPNT  subroutine  by  adding  additional  time  to  the  time  used 
to  schedule  the  XFER  event;  the  others  would  then  locate  the  target  later 
than  the  original  observer. 

XFER  is  the  subroutine  that  has  been  added  to  accomplish  this 
transfer. 
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The  misc  file  contains  an  additional  logical  value  that  is  read  in 
the  subroutine  RDMISC  and  used  to  set  the  value  of  the  variable  xxfer  which 
is  used  as  a  flag  to  call  or  not  call  transfer  in  PINPNT. 

The  other  type  of  detection  is  what  one  gets  with  looking  for  the 
targets  with  binoculars,  IR  devices,  etc.  without  the  aid  of  a  firing 
signature. 

Since  the  transfer  of  pinpoint  seemed  to  be  a  help  in  some 
instances  and  since  VHSIC  seem  to  offer  considerable  possibilities  for  fire 
control  improvement,  we  considered  the  possibility  that  any  target  that  is 
detected  could  be  made  available  to  everyone  on  the  side  of  the  observer. 

This  was  accomplished  simply  enough  by  adding  a  change  to  the 
subroutine  DETECT  that  is  similar  to  that  added  to  PINPNT.  In  the  version 
of  the  program  in  Appendix  A,  the  variable  xxfer  either  initiates  the 
transfer  of  both  types  of  detection  or  inhibits  the  transfer  of  any 
detections. 

As  alluded  to  above  in  Section  1,  the  availability  of  all  these 
targets  did  not  produce  the  Nirvana  anticipated.  Transfer  of  pinpoint 
detections  in  situations  with  low  pinpoint  ability  seems  desirable,  but  the 
gross  exchange  sometimes  seems  to  produce  problems  and  sometimes  seems  to 
help.  Also  as  mentioned  above,  the  prioritization  of  targets  and  the 
ability  to  keep  one  target  from  saturating  the  attention  of  the  force  need 
much  further  consideration. 

Issues  such  as,  the  prioritization  of  targets,  the  disengagement 
algorithms,  the  ability  to  transfer  targets  to  vehicles  that  can't  see  the 
targets  at  that  time,  and  the  ability  of  a  moving  vehicle  to  remember 
stationary  targets  when  the  moving  vehicle  loses  line  of  sight  contact  all 
need  further  consideration.  The  few  limited  results  we've  gotten  in  this 
regard  have  not  been  satisfactory  in  the  sense  of  providing  understanding. 

A  good  bit  more  work  is  still  needed. 

2.7  Changed  Handling  of  Shared  Targets 

TANKWARS  has  an  option  that  keeps  a  tank  from  firing  at  a  target  if 
another  tank  is  already  firing  at  that  target.  To  accomplish  this  option, 
one  sets  a  variable  called  share()  to  true  to  prohibit  firing  at  the  same 
target.  We  felt  that  it  was  extreme  to  prohibit  such  firing,  but  that  is 
was  still  desirable  to  be  able  to  limit  it.  A  change  was  added  to  the 
subroutine  PRIORN  so  that  if  share()  =  true,  the  priorities  are  adjusted  so 
that  a  target  being  attacked  has  a  lower  priority  than  anyone  not  being 
attacked.  The  targets  being  attacked  are  ordered  in  the  same  way  among 
themselves  as  are  the  targets  not  being  attacked.  Again  the  results  of 
using  this  option  have  been  mixed,  and  more  thought  about  the 
prioritization  scheme  is  indicated. 

2.8  Busy  Bug 


There  is  a  bug  in  TANKWARS  in  the  subroutine  VANTER  that 
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handles  the  vanishing  of  vehicles  in  terrain.  If  the  tank  has  made  the 
decision  to  engage  a  target  and  then  schedules  a  firing,  it  has  set  a  flag 
busy()  =  true.  When  the  scheduled  firing  routine  is  reached,  the  busy  flag 
is  set  to  false,  and  a  firing-on-target  flag  is  set.  If  the  firing  tank 
loses  line  of  sight  by  vanishing  in  terrain  before  the  scheduled  firing, 
VANTER  cancels  the  firing  but  does  not  reset  the  busy  flag  to  false.  Since 
the  firing  that  would  have  reset  the  flag  has  been  canceled,  the  busy  flag 
is  never  reset  during  the  engagement,  and  the  tank  can  no  longer  select  any 
targets  after  is  reappears.  We  added  a  fix  that  sets  busy  to  false  in 
VANTER.  This  resulted  in  a  sizable  improvement  in  the  general  performance 
of  the  attacking  force.  The  defender  has  little  occasion  to  use  VANTER. 

2.9  Know  Bug 

There  is  another  bug  in  version  2  of  TANKWARS  which  seems 
to  be  caught  in  version  2a.  The  variable  know()  is  used  in  the  subroutine 
PRIORT  to  distinguish  between  a  target  that  has  been  hit  and  one  that  has 
been  missed.  However,  the  value  know()  =  1  was  not  set  in  the  subroutine 
MAYHIT  if  the  target  was  missed  (knowQ  =  2  was  set  if  the  target  was  hit). 
We  added  a  correction  for  this  in  MAYHIT. 

3.  Possible  Advantages  of  GROUNDWARS 

GROUNDWARS  is  more  actively  maintained  and  is  among  other  things  a 
newer  version  of  TANKWARS.  It  also  has  some  additional  features  including 
the  availability  of  a  very  large  output  option  for  detailed  data  on  the 
combat  that  might  be  useful  for  some  diagnostics,  the  ability  to  handle 
subgroups  of  attacking  vehicles,  the  inclusion  of  the  effects  of  artillery 
fire,  an  improved  smoke  model,  and  most  of  all  an  improved  acquisition 
model  which  is  claimed  to  produce  a  factor  of  six  or  seven  improvement  in 
runtime.  This  acquisition  model  replaces  the  step-by-step  model  in 
TANKWARS  with  a  model  that  predicts  when  the  detections  will  occur  if  at 
all.  Note  that  this  is  not  a  different  physical  model  but  rather  a 
different  mathematical  formulation  of  the  solution. 

On  the  other  hand,  as  with  TANKWARS,  GROUNDWARS  also  needs  an 
ability  to  handle  leapfrogging,  an  ability  to  handle  moving  meeting 
engagements,  the  ability  to  handle  multiple  weapon  types,  and  some  rework 
of  the  assignment  of  priorities  for  engagement  and  disengagement.  Also  Dr. 
Howe  wishes  to  look  at  the  synergistic  effects  among  multiple  hits  on  a 
target  which  will  require  a  new  vulnerability  model  for  either  TANKWARS  or 
GROUNDWARS. 


4.  What  Next 


The  next  major  step  that  is  planned  is  to  develop  a  more  complete 
model  which  can  handle  combined  arms  in  essentially  the  same  detail 
as  TANKWARS  or  GROUNDWARS  and  address  the  additional  issues  mentioned 
above.  This  could  be  a  natural  (but  very  large)  step  forward  from 
developing  the  ability  to  have  more  than  one  armored  weapon  system  on  each 
side,  the  ability  to  do  leapfrogging,  and  the  ability  for  handling  more 
dynamic  meeting  engagements. 
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While  there  is  that  longer  term  goal,  there  is  also  a  shorter  term 
set  of  questions  to  which  Dr.  Howe  needs  answers  which  include: 

*  The  ability  to  model  advanced  C3  concepts  including  the  ability 
to  remember  targets,  the  ability  to  transfer  targets,  and  a  rework 
of  the  assignment  of  priorities  for  engagement  and  disengagement. 

*  A  new  vulnerability  model  that  can  account  for  possible 
synergistic  effects  among  impacts  on  the  target. 

*  The  ability  to  handle  a  meeting  engagement  between  two  forces 
that  are  moving. 

*  The  ability  to  have  more  than  one  weapon  systems  on  each  side. 

In  particular,  the  ability  to  have  both  missiles  and  guns  on  each 
side. 

There  are  two  interesting  approaches  to  the  longer  term  goal: 

(1) .  to  start  from  scratch  and  build  a  model  using  the  "C" 
programming  language  and  particularly  using  "structures"  to  represent  the 
fundamental  elements  (weapon  systems,  vehicles,  etc.) 

(2) .  to  build  on  existing  code  (here  it  would  seem  desirable  to  use 
GROUNOWARS,  specially  in  view  of  the  new  version  4.0,  it's  state  of 
maintenance,  and  it's  acceptance  by  AMSAA.).  Note  here  that  this  would 
most  certainly  require  major  surgery. 

CONCLUSIONS 


It's  interesting  that  both  AMSAA  people  and  the  author  arrived  at 
the  conclusion  that  structures  were  particularly  appropriate  in  this 
application.  AMSAA  is  presently  developing  yet  a  newer  combat  model  using 
"C".  In  terms  of  the  final  (if  one  could  say  that  there  ever  is  a  final 
version  of  such  a  code)  code,  the  use  of  structures  would  be  the  more 
powerful  approach  and  would  be  more  amenable  to  the  modifications  that  are 
always  required  to  consider  new  weapon  and  tactical  concepts.  The 
difficulty,  of  course,  with  such  an  approach  is  that  there  will  be  little 
short  term  benefit  of  such  an  effort.  And  the  questions  to  which  the  above 
changes  could  produce  answers  would  have  to  wait. 

The  second  approach  has  the  advantage  that  the  things  to  be  added 
to  the  code  could  produce  useful  results  as  they  are  added.  However,  where 
major  surgery  is  necessary  this  will  not  be  so  likely.  On  the  other  hand, 
the  use  of  FORTRAN  and  the  basic  architecture  of  GROUNOWARS  which  relies  on 
numerous  globally  declared  arrays  could  make  generalization  of  the  code 
very  difficult  and  could  make  the  resulting  code  very  opaque  and  very 
awkward  to  adapt  to  changing  concepts. 

It  is  the  author's  opinion  that  the  first  approach  is  the  proper  one.  At 
the  same  time  it  may  be  possible  to  provide  some  short  term  ability  without 
too  much  difficulty  by  adding  a  new  prioritization  scheme  and  a  more 
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complete  vulnerability  model  to  existing  code.  The  exercise  of 
particularly  the  prioritization  scheme  might  be  valuable  in  the  development 
of  a  new  code  by  providing  a  mechanism  for  discussing  tactics  with  the  User 
with  a  view  toward  a  better  representation  of  tactics  in  the  new  code. 

Certainly  even  if  the  first  approach  is  taken,  TANKWARS  and 
GROUNDWARS  will  valuable  sources  of  ideas. 
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APPENDIX  A 


TANKWARS  PROGRAM 


Intentionally  left  blank. 


Appendix  A 
TANKIARS  Prograa 


c  VXiX 

c  coaeon. h  file 

c  - 

c  coaaon  block  declarations  for  Extended  Coibat  Sieulation  (Tanklars  II) 

iaplieit  integer(i-n) ,  real  (a-h,o-z) 

paraaeter  (NN=20) 
character«4  color 
character^  kview 
integer  ALL,  NULL,  FLS  TGT 
integer  FD,  HD,  FE 
integer  TURRET,  HULL 
integer  BLU,  RED 
integer  UEETNG,  RATTAK,  8ATTAK 
integer  ALIVE,  UK ILL,  FKILL,  UFKILL,  IKILL,  KKILL 
integer  SLOING,  STATNY,  ACCELG,  UAXYL 
integer  scene,  tactic,  prevrd,  aray 
logical  cansee,  fot,  aot,  kneels 
logical  busy,  eapty,  foes,  los,  seen,  serchg,  repeat 
c  Change  by  HLReed 

logical  istest,  share,  xxfer,  aeaory 
real  INFINT 

c  Change  introduced  by  H.L.Reed  on  3  Uar  89  to  allow  overwatch 

c  tanks  to  be  added  to  the  attacking  force.  See  also  changes  in 

c  subroutines  cango,  deplo2,  in i t2,  and  input. 

logical  inwatch 
c 

coaaon  /aspekt/  angles(15),  pangle(IS),  iangd 
coaaon  /chare  /  color(2),  kview(2) 
coaaon  /consts/  PI,  TIOPI,  DEG,  VN0RTH(3) 
coaaon  /eonst2/  ALL,  NULL,  FLS  TGT, 

1  FO,  HD,  FE,  TURRET,  HULL,  BLU,  RED,  UEETNG,  RATTAK, 

2  BAnAK,  ALIVE,  UKILL,  FKILL,  UFKILL,  IKILL,  KKILL  , 

3  SLOING, STATNY, ACCa.G,UAXVL, INFINT 

coaaon  /contrl/  nreps,  keyd(20),  keya(20),  scene,  taax,  aeth  sa 
coaaon  /epath  /  naaxt(2) , accel (2) , decs  I (2) , ishtfs(2) , 

1  speed (2) ,  angle(2),  accaax(2),  wy I th (2) ,  aapl(2) 
coaaon  /cranda/  iranda,  jranda 
coaaon  /eshot/  kshot(2,20) 
coaaon  /ctrace/  trace 
logical  trace 

coaaon  /endgaa/  sysdia(2,8),  nang,  ndisp 
e  coaaon  /errors/  ssrgs(2,10),  sargs(2,10),  velas(2,20), 

c  1  sserrs(2, 18, 10) ,  saerrs(2, 16, 10) ,  addons (2, 2, 20) ,  nadds(2) 
coaaon  /error2/  rex,  rey,  reliab(2) 
c  Change  by  HL  Reed  to  allow  transfer  added  xxfer.  Also  added  the 
c  variables  which  follow  xxfer  but  are  not  being  used  at  this  tiae. 

coaaon  /f eye  I  e/  nrds(2),  nrpt(2),  nipods(2),  nrpb(2),  tactic(2), 

1  tof (2,8) ,trelod(2) ,  tfirst(2,8),  taedin(2),  tfixed(2,8), 

2  rof(2),kind  rd(2) ,tbuap(2) ,nbuap(2) ,thide(2) , tain (2) , 

3  npr i or (2) ,nchans(2) ,share(2) ,xxfer(2) , aeaory (2) , 

4  tpopl (2) , tpop2 (2) , tpop3 (2) , nstatn , naove, rstatn , 

5  raove,  ttotal,  nstatf,  naovef,  rstatf,  raovef 
coaaon  /n  sys/  ntanks(4,6),  nblu,  nred 

c  Change  by  HLReed  added  pntiae  for  aedian  pinpoint  tiae 

coaaon  /sensor/  psense(2,8),  pinf in(2,3, 10) ,  tbar(2,3, 10) , 

1  ndets(2) ,  tlook(2),  pinp(2),  repeat,  recknz(2),  pfalse(2,2), 

2  pntiae (2) 

coaaon  /states/  aray(NN),  aspect(HN,2) ,  cansee(NN), 

1  busy(NN) ,  eapty (NN) ,  fot(NN,NN) ,  foes(NN.NN),  ichg(NN), 

2  kneeal(NN),  kneels(NN.NN) ,  know(NN.NN) , 

2  life(NN),  los(NN.NN),  aot(NN.NN),  aotion(NN),  aslf ly(2,NN,S) , 

3  nhot(NN) ,  nbrst(NN),  ndet(NN),  nrd(NN),  nrib(NN),  mpod(NN), 

4  nrot(NN) ,  nrtgt(NN),  nchan(NN),  prevrd(NN),  rgvis(3,NN), 

5  seen(NN.NN),  serchg(NN),  tf i re(NN.NN) ,  tfire2(NN),  vbx(NN), 

6  vby(NN),  t0 (NN) ,  x0(NN),  y0(NN),  vx0(NN),  vy0(NN), 

7  xp(NN) ,  tlast(NN) 

coaaon  /state2 /  idecoy(NN),  iflash(NN),  ndecoy(2),  nflash(2) 
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c  Change  for  output  HLReed 

coaeon  /stats/  statb(8),  aystat(4) 

coaaon  /tstore/  a(1000),  iholy 

coaaon  /v ars8/  irginc,  rgincr,  rginc2 

coaaon  /csaoke/  tsaoke(20) ,psaoke(20,3) ,  invisb 

coaaon  /saokel/  touti 1(21, 5) , toutvl(21, S) , tout i (21,5) , 

1  toutv  (21 . 5) ,  t  i  n  i  (21 , 5) ,  t  i  nv  (21 , 5} ,  ptb  I  (21) ,  rtb  I  (S) 

coaaon  /vl8a /  krep 
coaaon  /vl7a/  istest 

coaaon  /v23/  nwave,  nwaves,  nsurv,  naval,  nused(3000),  nreps3, 

1  statc(8},  noaaao,  loaaao,  noaao2,  loaao2 
coaaon  /v2 4/  nstats(5,2) 
coaaon  /whore/  ain  rg,  aax  rg,  inc  rg 
coaaon  /where2/  nrg,  rg0,  rg,  s(3),  vt(3),  vf(3) 
coaaon  /fitnes/  qu i t (2 , 2) ,  alloc(2,8),  fit(NN,6) 
c  Change  introduced  by  H.L.Reed  on  8  Uar  89  for  overwatch 

coaaon  /ovrwtch/  nwatch (3) , inwatch (NN) 

c  V7.2 

c  MAIN  ROUTINE 

c  9  Main:  read  input  and  siaulate  scenarios, 

include  'coaaon. h’ 

1  foraat  (’  The  Sustained  Coabat  Model:  Tank  Vars  II.’,/, 

1  1  Written  by  Fred  Bunn  (ph  (301)  278-6648,  autovon  298-6648)',/, 

2  ’  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MO',/, 

3  '  Version  7.2  Created  10/24/88.',/) 

c 

call  input 
00  20  scenes 1,3 
call  forces 
20  CONTINUE 

ENO 


c  V7.1 
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SUBROUTINE  A80RT  (t.f irer, tgt) 

c  6  Abort:  abort  asl  fro»  fir or  to  tgt  (to  ail  tgts  if  tgt=0) 
include  ’coeaon.h1 
logical  defndr 

integer  arayf,  arayt,  firer,  tgt 

1  fonat(f8.2,lx,a4,  i  3 ,  ’  asl  for  ' ,  a4,  i3 , 1  aborted.’) 

2  foraat  (’A80RT:  f i rer, tgt, i ,aslf ly=’ ,4i3) 

3  foraat  (’ABORT:  asl  approaching  tgt’,*i3) 
c 

if  (trace)  print  «,’>abort’ 
arayf  *  aray(firer) 
arayt  =  3-arayf 
00  20  i=l , 5 

c  Check  all  S  aissile  pointers  for  this  firer 

asl  =  asl  f I y (arayf .firer, i) 
if  (keys (19) .gt.0)  print  2,  firer,  tgt,  i,  asl 
IF  (asl. gt.0)  THEN 

c  Uissile  found  (pointer  is  non-zero) 

asl  tgt  =  a  (as  1-1) 
if  (keya(19) .gt.0)  print  3,  asl  tgt 
IF  (tgt.eq.0  .or.  tgt. eq. as  I  tgt)  THEN 
c  Abort  this  aissile 

kshot(arayf ,3)  *  kshot(arayf ,3)»i 
call  cancel  (asl,  elUPCT,  NULL) 
asl  fly(arayf, firer, i)  *  0 

If  (asl  tgt.ne.FLS  TGT)  aot(f irer,asl  tgt)  =  .false, 
c  Release  area  for  storage  of  aissile  data 

a  (asl)  =  -a (asl) 

if  (keyd(l) .ge.2)  print  1, t, color (arayf ) ,f irer, 

1  co I  or  (arayt) , as  I  tgt 

c  Pop-down  to  reload  if  defender,  pod  eapty,  i  fully  alive 

defndr  *  (scene. eq.BATTAK  .and.  arayf .eq. RED)  .or. 

1  (scene. eq.RATTAK  .and.  arayf .eq.BLU) 

if  (defndr. and. eapty(f irer) .and. (I ife(firer) . le. ALIVE)) 

2  call  skedul(t,f irer, ePOPON, NULL) 

ENDIF 

ENDIF 

20  CONTINUE 

if  (trace)  print  «,’<abort’ 

END 
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SU8SWUTINE  ACCB.F  (t,  firer) 
c  9  Acc elf:  s isolate  tank  starting  to  accelerate, 
include  ’ coaeon. h' 
integer  firer 

1  foreat  (f8.2,lx,a4, i3, '  speed  up\9x,’(was  slowing)’) 

2  foreat  (f8.2,l*,a4,  i3,  ’  speed  up’^x.’^was  halted)’) 

3  foreat  (f8.2,lx,a4,  i3, '  speed  up\9x,'(was  speeding  up)') 

4  foreat  (f8.2,lx,a4, i3, ’  speed  on\9x,’(is  cruising)') 

c 

if  (keyd(4)  .gt.0)  print  *,’>accel’ 

if (I ife(f irer) .ne.FXILL.and. invisb.eq. l.and.knceal (f irer) .ne.FD) 
1  call  skedul  (t,f  irer, ’vanish’.NUl) 
narey  a  arey (firer) 

IF  (eot i on (f i rer) . eq . SLQWNG)  THEN 
c  Previous  eot ion  was  slowing 

if  (keyd(l) .ge.2)  print  I,  t,  color(naray),  firer 
cal  I  path(f irer,t,eotion(f irer),8.01xly1vx1vy) 
dt  =  (speed(nareyj-vy)/accel (narey) 
call  skedul  (t»dt,  f  i  rer,  'eaxveP  .NULL) 
eot  ion  (firer)  =  ACCELG 
ELSE  IF  (aotion(f irer) .eq.STATNY)  THEN 
c  Previous  eotion  was  stationary 

if  (keyd(l)  ge.2)  print  2,  t,  color(naray) ,  firer 
cal  I  path(f i rer, t, sot  ion (f i rer) , 0. 0, x,y , vx, vy) 
c  schedule  tiee  full  velocity  reached  (aax  vet) 

dt  a  speed (narey) /acce I (arey (firer)) 
cal  I  skedul (t»dt,f irer, 'eaxvel ’.NULL) 
sotion(f  i  rer)  a  ACCELC 
ELSE  IF  (eot i on (f i rer) . eq . ACCELC)  THEN 
c  Previous  eotion  was  accelerating 

if  (keyd(l) .ge.2)  print  3,  t,  color(narny),  firer 
ELSE  IF  (eotion(f irer) .eq.UAXVL)  THEN 
c  Previous  eotion  was  cruising  at  aax  velocity 

If  (keyd(l) .ge.2)  print  4,  t,  eolor(naray) ,  firer 
ENDIF 

if  (keyd(4) .  gt.0)  print  «(’<aceer 
END 


c  V7.1 


FUNCTION  ANCLF  (a,  b) 

c  9  Anglef:  find  angle  between  two  vectors. 

diaension  a(3),  b{3) 
c 

vabsa  =  sqrt(  dot(  a, a  )  } 
vabsb  =  sqrt(  dot(  b,b  )  ) 
dotab  a  dot C  a,b  ) 

CHANGED  1  Apr  86.  Next  line  replaced  by  3. 
c  da  =  acos(dotab/(vabsa*vabsb)) 

da  a  dotab/ (vabsawvabsb) 
da  =  aainl(l. ,aaaxl(-\. ,da)) 
da  =  acos(da) 

r3  *  a(l)«b(2)  -  a(2)eb(l) 
anglef  a  -sign(da,r3) 

END 

c  V7.1 


A  -  6 


FUNCTION  ANGSUU  (a,  b) 

c  3  Angsua:  add  2  angles  and  adjust  answer  to  lie  between  *-PI. 
c  *  a»b 

10  IF  (c. It. -180.)  THEN 
c  =  c-380. 

GOTO  10 

B.SE  IF  (c .gt . 180.)  THEN 
e  =  e-360. 

GOTO  10 
ENDIF 

c  angle  is  adjusted, 

angsua  =  c 
END 

c  V7.3 


SUBROUTINE  APPEAR (t, tgt, f irur) 

Appear:  if  tgt  appears  treat,  otherwise  reschedule  appearance 

include  ’coieon.h' 

integer  tgt.firer,  arayf,  anyt 

coeeon  /terane/  d (40) ,  xold(20),  yold(20),  d i at (20) ,  iseg(20) 
rss(x.y)  =  sqrt(x«x»y*y) 

forsat(f8.2, lx,a4, i 3 , ’  appears  1 ,9x , ’ (x=* ,f8. 1, 1  y=’ ,f8.1, ’) ') 

foraat(f8.2,lx,a4, i3, ’  LOS  to  ',14,13,’  starts.’) 

if  (trace)  print  «,’>appear' 
arayt  =  aray(tgt) 
arsyf  =  3-arsyt 
IF  (invisb.eq.l)  THEN 

if (speed (arayt) . I e.0.) print  *, 'APPEAR:  arsyt,speed=’ ,arsyt, 
speed (arayt) 

if  (speed (arayt) . I e.0.)  stop 

cal  I  path  (tgt, t, sot  ion (tgt) ,0.2,x,y,vx,vy) 

Terrain  causes  interaittent  LOS. 
travel  =  rss(x-xold(tgt) ,  y-yo I d (tgt) ) 

IF  (travel .gt.dist(tgt))  THEN 
Tgt  is  no  longer  aasked  by  terrain 

if  (keyd(l) ,gt. 1)  print  l,t,color(ar»yt) ,tgt,x,y 
xold(tgt)  =  x 
yold(tgt)  =  y 
iseg(tgt)  =  iseg(tgt)»l 

if  (iseg(tgt) .gt.40)  iseg(tgt)=iseg(tgt)-40 
dist(tgt)  -  d(iseg(tgt)) 
call  aprter(t, tgt, f Irer, FE) 

Schedule  next  disappearance 
dt  a  dist (tgt) /speed (arayt)  ♦  0.01 
call  skedul  (t»dt,  tgt,  ’vanish'  ,NU-L) 

B.SE 

Still  aasked  by  terrain,  so  reschedule  aask  end 
IF  (life (tgt) .eq. ALIVE)  THEN 
dt  =  (dist(tgt)  -  travel)  /  speed(arayt)  ♦  0.01 
call  skedul  (t«dt, tgt, ’appear’ , NULL) 

END  IF 
ENDIF 
S.SE 

Tgt  is  no  longer  aasked  by  saoke 

if  (keyd(l) .gt.l)  print  2,t,color(3-arayt) ,f irer, 
color (arayt), tgt 
call  aprsak(t, tgt.firer) 

Schedule  next  disappearance 
r  =  rgf(t,f irer, tgt) 
p  =  ranu(0) 
pout  a  ranu(0) 

IF  (kview(RED) .eq.kview(BLU))  THEN 
IF  (arayf.eq.8LU)  THEN 
IF  (kview(arayf) .eq. ’I’)  THEN 
dtin=tdintp(ptbl , rtb I ,tini ,p,r,21,S) 
dtout=tdintp(ptbl , rtbl , touti , pout,  r,  21, 5) 

EiE 

dtin=tdintp(ptbl ,  rtbl  ,tinv,p,r,21,5'. 
dtout=tdintp(ptbl , rtbl ,toutv,pout,r,21,5) 

ENDIF 

ca  1 1  skedu I (t«dt i n , tgt , 1  van i sh ’ , f i rer) 
cal  I  skedu I (t*dtin,f irer, ’vanish’, tgt) 
cal  I  skedul (t*dtin»dtout, tgt, 1  appear ’,f irer) 
cal  I  skedul  (t«dtin«dtout,fi rer, ’appear’ , tgt) 

ENDIF 

BJE 

IF  (kvise(iriyf).eq.’V’)  THEN 
dtinv*tdintp(ptbl , rtbl ,tinv,p,r,21,5) 
dtoutv*td i ntp (ptb I , rtb I , toutv , pout , r , 21 , 5) 
dtout i =td i ntp (ptb I , rtb I , tout i , pout , r , 21 , S) 
dt i n i =dt i nv« (dtoutv-dtout i ) «0 . 5 
ca 1 1  skedu I (t*dt i nv , tgt , 1  van i sh ’ , f i rer) 
cal  I  skedul (t*dtinv«dtoutv,tgt, ’appear’, f irer) 
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call  skedul  (t»dtini ,f irer, ’vanish’.tgt) 
call  skedu I (t»dtout i »dt i n i ,f irer, ’appear’ ,tgt) 
ENDIF 
ENDIF 
£N  OIF 

if  (trace)  print  «,’<appear’ 

END 

V7.1 
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SUBROUTINE  APRSMK(t,tgt,f irer) 
e  9  Aprsak:  Tgt  appears  out  of  saoke,  reset, 
include  ’coaaon.h' 
integer  tgt, f irer 

coaaon  /terane/  d(40),  xo I d (20) ,  yold(20),  dist(20),  iseg(20) 
c 

if  (trace)  print  *,,>apprsak‘ 
naray  =  any  (tgt) 

c  Restore  all  I ines-of-sight  involving  tgt 

los(f irer, tgt)  =  aray(f irer) .ne.naray 
c  Turn  search  on  if  it  is  off 

IF  (.not. repeat)  THEN 
repeat  =  .true. 

call  skedu I  (t» .01,0, 'search' , NULL) 

ENDIF 

if  (trace)  print  «,,<aprsak‘ 

END 

c  V7.1 
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e  0 


1 

c 


c 


20 

c 


e  V7.1 


SUBROUTINE  APRTER(t,tgt,f i rer, jexpos) 

Apprter:  Tgt  has  appeared  fros  behind  terrain,  reset, 
include  'coaaon.h' 
integer  tgt,firer 

coaeon  /terane/  d(40),  xold(20),  yold{20),  dist(20),  iseg(20) 
foraat(f8.2,lx,a4, i3, ’  aprters  ’ ,9x, ' (x=‘ ,f0.1, ‘  y*' ,f8. 1 , ’) ’) 

if  (trace)  print  e/Japrter1 
narsy  =  any  (tgt) 
knceal(tgt)  =  jexpos 

Restore  all  I ines-of -sight  involving  tgt 
00  20  i-l,nblu*nred 
IF  (knceal (i) .ne.FD)  THEN 
los(tgt,i)  =  arsy(i)  .ne.naray 
los(i,tgt)  =  arey(i) .ne.narsy 
ENDIF 
CONTINUE 

Turn  search  on  if  it  is  off 
IF  (.not. repeat)  THEN 
repeat  =  .true. 

call  skedu I (t* .01,0, 'search' .NULL) 

ENDIF 

if  (trace)  print  «1,<aprter' 

END 
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SLOCK  OATA  8LKDAT 
include  'coaaon.h' 

data  color,  pi,  twopi,  deg 

1  /‘Slue’,  ’Red  *,  3.141S92854,  8.283185308,  57.29577951/ 
data  VNORTH  /».,  1.,  0./ 
data  ALL,  NULL,  FLS  TGT  /0,  0,  -1/ 
data  FD,  HO,  FE  /l,  2,  3/ 
data  TURRET,  HULL  /l,  2/ 
data  8LU,  RED  /l,  2/ 
data  MEETNC,  RAnAX,  BATTAK  /l,  2,  3/ 

data  ALIVE,  UKILL,  FKILL,  IFKILL,  KILL,  KKILL  /l, 2, 3,4,5, 6/ 
data  SLOINC,  STATNY,  ACCH.G,  UAXYL,  INFINT 

1  /  1,  2,  3,  4,  l.e3S/ 

data  keyd,  keya  /4 0*0/ 

ENO 

c  V7.1 
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SUBROUTINE  BOUNDS  (naray,  box,  angl I,  rl,  r’) 
c  S  Bounds:  find  the  horizontal  bounds  of  hull  or  turret, 
c  Def i n i t i ons : 

c  angl  I  •  angle  off  the  nose  of  the  box  (rad), 

c  box  -  1  seans  turret  box,  2  scans  hull  box. 

c  naray  -  1  scans  blue  firers,  2  seans  red  firers. 

c  c,  s2,  s3  -  tesporary  variables, 

c  rl,  r2  -  left  and  right  boundaries  of  boxes  (a), 

include  'coaaon.h1 
integer  box 
c 

if  (trace)  print  «,,>bounds‘ 
c  initialize 

teap  *  (angl  Utwopi)/twopi 
theta=  (teep-aint(teap))*twopi 
c  theta  *  aaod  (angl  Utwopi ,  tsrop i ) 

c  *  sysdis(naray,4«(box-l)'*2)  »  cos(theta) 
s2=  sysdia(naray,4«(box-l)*3)  «  sin(theta) 
s3=  sysdis(narsy,4»(box-l)*4)  *  sin  (theta) 
c 

IF  (theta . I e . 0 . 25*t«op i )  THEN 
c  case  0  <  theta  <*  90 

rl  »  -s2  -  e 
r2  *  s3  ♦  c 

ELSEIF  (theta. Ie.0.5«t*opi)  THEN 
c  case  90  <  theta  <3  180 

rl  3  -s2  ♦  c 
r2  *  s3  -  c 

Q.SEIF  (theta . I e . 0 . 75« tsop i )  THEN 
c  case  180  <  theta  <3  270 

rl  3  s3  •  c 
r2  3  -s2  -  e 
ELSE 

c  case  270  <  theta  <3  360 

rl  a  s3  -  c 
r 2  3  -s2  ♦  c 
END  IF 

if  (trace)  print  t.^bounds1 
END 


c  V7. 1 
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SUBROUTINE  CANCQ.  (I.  ict,  it) 
c  0  Cancel:  cancel  ’act'  events  for  *1*  entity, 
c  (all  events  if  acts”) 

c  Definitions  of  local  variables: 

c  ■  -  pointer  to  previous  event 

c  n  -  pointer  to  current  event  being  considered 

include  ’clock.h' 
logical  is  >hat,  is  oho,  is  whoa 
characters  act 

1  foreat(9x,  ’cancel  * , i 3 , *  * , afi , i 3, '  at  tiee’,fB.2) 

c 

13  0 

n  =  nxevnt 

10  IF  (n.ne.0)  THEN 

c  Continue  unti I  n=0 

is  who  3  I  . eq.who(n) 

is  «hat  3  act.eq.ehat(n)  .or.  act.eq.’all  1 
is  «hoe  3  it.eq.whoe(n)  .or.  it.eq.0 
IF  (is  *ho  .and.  is  what  .and.  is  »hoa)  THEN 
c  Then  reeove  event 

if  (prflag  )print  1,  I,  act,  it,  vhen(n) 
if  (a.eq.0)  nxevnt  3  next(n) 
if  (a.ne.0)  next(e)  3  next(n) 
next(n)  =  nxidle 
nxidle  3  n 

if  (a.eq.0)  n  =  nxevnt 
if  (a.ne.0)  n  =  next (a) 

ELSE 

c  Don't  reaove  event.  Shift  to  next  event, 

a  3  n 

n  »  next(n) 

ENDIF 
GOTO  10 
ENDIF 
END 


c  V7.1 
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LOGICAL  FUNCTION  CAN  GO  (firer,  t) 
c  6  Can  go:  True  iff  is  stationary  and  can  aove. 
include  'coaaon.h1 
integer  firer 

logical  is  atkr,  a  alive,  faster 
c 

naray  =  aray (firer) 

is  atkr  =  (naray.eq.8LU  .and.  scene. eq.BATTAK)  .or. 

1  (narey.eq.RED  .and.  scene. eq.RATTAK) 

a  alive  -  I i f e (f i rer) .eq. ALIVE  .or. 

1  I ife(f irer) .eq.FKILL 

faster  =  (aotion(f irer) .eq.STATNY  .or. 

1  aot i on (f i re  r) . eq . SLOVNG) 

c  Change  introduced  by  H.L.Reed  on  8  Man  89  to  a  I  toe  overwatch  tanks 

c  to  be  added  to  the  attacking  force.  See  also  changes  in  subroutines 

c  deplo2,  init2,  and  input  and  in  coaaon.h. 

can  go  =  is  atkr  .and.  a  alive  .and.  faster  .and. 

1  (.not.  inwatch (f i rer)) 

END 

c  V7.1 
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SUBROUTINE  CREATE  (n,  ient) 

c  8  Create:  create  a  temporary  entity,  (a  bullet  or  esl) 

c  purpose  -  this  routine  'creates'  a  temporary  entity,  chat  it 

c  actually  does  is  find  space  to  store  the  attributes  of  that 

c  entity  and  reports  back  the  index  of  the  first  storage  word 

c  as  the  entity  nueber. 

c  definitions  - 

c  a  is  the  vector  used  to  store  attributes  in. 

c  i  this  is  the  index  of  the  storage  space  ve  are  currently 

c  looking  at. 

c  ient  is  the  index  of  the  storage  space  there  the  attributes 

c  till  be  stored,  it  is  also  the  nuaber  that  till  be  used 

e  to  identify  the  teeporary  entity  created, 

c  isdone  the  routine  is  done  if  isdone  *  1  and  it  is  not  done 

c  if  isdone  =  0.  the  only  reason  for  the  routine  to  be  done 

c  is  if  it  finds  space  to  store  the  attributes  in. 

c  istart  this  is  the  starting  point  for  the  search,  if  «e  get 

c  back  to  istart  eithout  a  find,  >e  have  a  storage  overload 

c  and  we  error  off. 

c  j  this  is  the  index  of  the  next  storage  space,  we  want 

c  to  look  at  it  with  the  possibility  of  catenating  it  to  the 

c  storage  space  beginning  at  i.e  n  the  nuaber  of 

c  attributes  to  be  stored. 

c  nreq  this  is  the  nuaber  of  spaces  required,  it  equals  the 

c  nuaber  of  attributes  plus  one  word,  this  one  word  is  used 

c  for  searching  purposes,  if  it  is  negative  (-abs(a)),  that 

c  indicates  that  the  -ext  a  words  are  being  used  to  store 

c  the  a  attributes  of  an  entity,  if  it  is  positive,  then  the 

c  next  a  word.*  -r  available  for  use. 

c  note  -  the  aaouit.  . .  storage  space  is  1000  words,  if  you  want  to 
c  increase  ‘.hi*  ju’ll  have  to  change  all  occurences  of  1000. 

c  also  no‘.^  in  initx  these  aust  be  set  -  a=0,  a(l)=1000., 

c  i=l. 

c 

l<-;,cal  trace 

coaaon  /contrl/  nreps,  key d (20) ,  key* (20),  scene,  taax,  aeth  sa 
coaaon  /ctrace/  trace 
*  coaaon  /tstore/  a (1000),  i 

1  foraat  (10x, ’CREATE:  Not  enuf  space  to  store’, 

1  iS,  ’  attributes.’) 

2  foraat  (10x, ’CREATE:  i,  j,  a(i),  a(j)  =’, 

1  /10x,2i5,2f 10.3) 

c 

if  (trace)  print  •,’>createl 
c  Initialize 

isdone  =  0 
istart  -  i 
nreq  =  n*l 

c  Find  eapty  space  in  the  a-array 

10  IF  (isdone.ne.l)  THEN 

c  Try  next  eapty  space 

c  Catenate  eapty  spaces  if  possible 

20  CONTINUE 

c  Find  next  space  (and  error  off  if  we’re  back  at  start) 

j  =  i»iabs(int(a(i))) 
if  (j  . gt.  1000)  j*l 

IF  ((j.eq.l)  .or.  (a(i).lt.0)  .or.  (a(j).lt.B))  THEN 
c  test  this  space  for  size. 

IF  (a(i) . It. float (nreq))  THEN 
c  aove  to  next  space 

i  *  j 

if  (i .eq. istart)  print  1,  n 
IF  (i .eq. istart)  STOP 
ELSE 

c  reserve  space 

isdone  *  1 
iteap  *  i«nreq 

if (a(i) .ne.f loat(nreq))a(iteap)  *  »(i)-f loat(nreq) 
a(i)  *  -nreq 
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ient  *  i 

'  *  j 

ENDIF 

B.SE 

c  do  catenation 

*(')  =  *(i)*a(j) 

■  IF  (a ( i ) . gt . 0 . 0  .and.  a(j).gt.0.0)  GOTO  20 
print  2,  i,  j,  a(i) ,  a(j) 

STOP 
ENDIF 
GOTO  10 
ENDIF 

if  (trace)  print  «,,<create’ 

END 

c  V7.1 
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SUBROUTINE  CRESET 

c  0  Creset  -  Reset  variables  used  by  create, 

coeaon  / tstore/  a(1000),  iholy 
paraaeter  (NN=20) 

00  20  i =2 ,1000 
a(i)=0.0 

20  CONTINUE 

a(l)=-NN 

a(NN.l)  *  1000-NN 

iholy=NN.l 

END 


e  V7.1 


A  -  18 


C  l 


c 


1 

2 

c 


SUBROUTINE  DAMAGE  (t,  I,  it,  injury) 

Damage:  schedule  effects. 

Changed  May  18,  1989  for  simplified  hit  and  kill  sodel,  HL  Reed 
include  'coeaon.h' 
character«2  kt(6) 

data  kt  / ’no‘ , 'U- 1 , !F- ’ , 'UF' , ’I- 1 , 'X- ’/ 
foraat(f8.2,lx,a4, i3,lx, 'Hits  1  ,a4,i3,'  (no  damage) . ') 
foraat(f8.2,lx,a4, i3,lx,a2, '-kills  ,,a4,i3) 

if  (trace)  print  •,’>daaage’ 
n=aray(I) 
a  =  3-n 

IF(keyd(l) .ge.2)  THEN 

if  (injury. eq.l)  print  l,t,color(n) , I,eolor(a) , it 
if  (injury. gt.l)  print  2,t,color(n) ,I,kt(injury) ,color(a) , it 


injold  =  I ife(it) 

IF  (injury. eq.KKILL  .and.  injury .ne. injold)  THEN 
c  Treat  first  catastrophic  kill. 

I ife(it)  =  KKILL 

call  daaagf (t, it, a) 

call  daaaga  (t, it) 

call  cancel  (it, 1 iki 1 1  ’.NULL) 

call  newtgt  (t,I,it) 

call  deaths (t) 

ELSEIF  (injury. ne. injold  .and.  injold. It. UFKILL)  THEN 
c  Treat  new  damage  (less  than  catastrophic). 

IF  (injury. eq.UKILL)  THEN 
if  (injold. eq.FKILL)  life(it)  =  UFKILL 
if  (injold. eq. ALIVE)  life(it)  =  MKILL 

if  (injold. eq. ALIVE  .or.  injold. eq.FKILL)  call  daaaga  (t , it) 
ELSE  IF  (injury. eq.FKILL)  THEN 
if  (in jold.eq. ALIVE)  1  ife(it)*FKllJL 
if  (injold. eq.UKILL)  I ife(it)=UFKILL 
call  daaagf (t, it, a) 

ELSE  IF  (injury. eq.VFKILL)  THEN 
if  (injold. It. UFKILL)  life(it)  =  UFKILL 
if  (injold. ne.UKILL)  call  daaaga  (t,  it) 
if  (injold. ne.FKILL)  call  daaagf  (t,it,a) 

ENOIF 

if  (I ife(it) .eq.kFKILL. and. injold. It. UFKILL) 

1  call  skedul  (t*tbuap(n) , it, ' iki 1 1  'NULL) 

ENOIF 

if  (trace)  print  «,'<daeage' 

END 

c  V7.1 
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SUBROUTINE  DAMAGE  (t,it,«) 

c  Dasagf  -  Discard  activities  due  to  firepoeer  kill, 

include  ’cotson.h' 

nrtgt  (it)  *  0 
nchan(it)  =  0 

c  Clear  any  guidance  channels  in  use  by  target, 

if  (kindrd(i/ .eq.4)  call  abort(t, i t , 0) 

00  40  j=l,nblu*nred 
fot(it,j)  =  .false. 

40  CONTINUE 

call  cancel (it, 'fire  ’.NULL) 
call  cancel (it, 'reload' ,NILL) 
call  cancel  (it, 'select' , NULL) 

IF  (I ife(it) .eq.FKILL  .and.  speed(a) .gt.0.0)  THEN 
call  cancel  (it, 'slowup' ,NU_L) 
call  cancel  (,t,'haU  '.NULL) 
call  cancel  (it,' acce l  ’.NULL) 
call  skedul  (t,  it,  'accel  ’.NULL) 
dt  =  thide(a) 

if  (a.eq.BLU  .and.  scene. eq.RATTAX)  dt  =  S.0 
if  (a.eq.RED  .and.  scene. eq. BA TTAK)  dt  =  5.0 
call  skedul  (t*dt, it, 'hide  ’,NU-L) 

ENDIF 

END 

c  V7  .I 
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SUBROUTINE  DAkAOl  (t,  it) 

c  9  Daaaga  -  Siaulate  aobility  kill  on  tha  tgt. 
include  ’coeaon.h’ 
logical  sos 

c  sos  -  stopped  or  s lowing 

if  (trace)  print  •,’>dasaga’ 

call  cancel  (it,  'aaxvel',  NULL) 

call  cancel  (it,  ’accel  ',  NULL) 

call  cancel  (it,  ’hide  ',  NULL) 

sos  =  vabs(vt) . ie.0.8  .or.  aotion(it) .eq.SLOWNO 

if  (.not. sos)  call  skedul  (t,  it,  ’slowup1,  NILL) 

if  (trace)  print  *,><daaaga’ 

END 

c  V7.1 
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SUBROUTINE  DEATHS  (t) 

c  0  Deaths:  Find  death  toll  on  each  side.  A  tank  is  considered 

c  dead  if  it  is  I-killed,  K-killed,  or  F-kiiled  k  hidden, 

include  ’coaaon.h* 
logical  deadl,  dead2 
integer  dead (2) 

1  foraat  ( i  3 ,  ’  Blu  dead, ',13,'  Red  dead.'} 

c 

if  (trace)  print  »,’>deaths' 
dead (BLU)  =  3 
dead (RED)  =  0 

c  Change  aade  by  H.L.  Reed  on  Uarch  31,  1989  to  keep  decoys  out 

c  of  the  sin  decision  and  to  keep  thea  out  of  the  exchange  ratio 

c  (see  also  finish. f) . 

DO  20  i =1 , (nblu-ndecoy(BLU)) 
deadl  -  I ife(i) .ge.IKILL 
dead2  =  knceal (i) .eq.FD  .and.  I ife(i) .ge.FKILL 
if  (deadl  .or.  dead2)  dead (BLU)=dead (BLU) ->1 
20  CONTINUE 

DO  30  i=nblu»l, (nblu^nred-ndecoyJREDJ-nsatchfseene)) 
deadl  =  I  if e(i) .ge.IKILL 
dead2  =  knceal (i) .eq.FD  .and.  I ife(i) .ge.FKILL 
if  (deadl  .or.  dead2)  dead(RED)=dead(RED).l 
30  CONTINUE 

if  (keyd(l) .go. 2)  print  l.dead 

if  ( (nb I u-ndecoy (BLU) ) . eq . dead (BLU)  .or. 

1  (nred-ndecoy (RED) -nwatch (scene)) .eq. dead (RED)) 

1  call  skedul  (t-5.  .NULL,  'finish1  .NULL) 
if  (trace)  print  (deaths’ 

END 

c  V7.1 
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SUBROUTINE  DEPL02(i south, j south, nsouth, i north, jnorth.nnorth) 
c  0  Deploy:  position  i  orient  all  tanks  at  beginning  of  engagement, 
include  1 coaeon. h’ 

2  foreatCtank’ ,  i2, 1  xi’.ffi.B,’  y=* ,f6.0, '  heading=’,f5.0, 

1  ’  speeds’ ,fS.l, 'a/s') 

if  (trace)  print  »,,)deplo2’ 
spacng  =  100.0 

c  Position  southern  tanks  on  the  x-axis 

vsouth  =0.0 

if  (scene. eq.BATTAK)  vsouth  =  speed (BLU) 

if  (scene. eq.RATTAK)  vsouth  =  speed (RED) 

x0(isouth)  =  -0.5«(nsouth-l)»spacng 

xp(i south)  =  x0(i south) 

y0(isouth)  =  0.0 

vy0(i south)  =  vsouth 

aspect (i south, TURRET)  =0.0 

aspect (i south, HULL)  =0.0 

n  =  i south 

c  Related  to  8  Uar  89  change 

inwatch (n)  =  .false. 

if  (keyd(l) .ge.2)  print  2,  n,  x0(n),  y0(n),  0.0,  vy0(n) 

DO  20  n=isouth»l,j south 
x0(n)  =  x0(n-l)  •  spacng 
xp(n)  =  x0(n) 
y0(n)  =  0.0 
vy0(n)  =  vsouth 

c  Change  introduced  by  H.  L.  Reed  on  8  Uar  89  to  allow  overwatch 

c  tanks  to  be  added  to  the  attacking  force.  See  also  changes  in 

c  the  subroutines  cango,  in it2,  and  input  and  in  coaaon.h. 

inwatch(n)  =  .false. 

IF (n. gt. (j south  -  naatch (scene)))  THEN 
inwatch(n)  =  .true. 
vy0(n)  *  0.0 
aotion(n)  *  STATNY 
knceal(n)  =  HD 
END  IF 

c  End  of  8  Uar  89  change 

aspect (n, TURRET)  =0.0 
aspect (n.HULL)  =0.0 

if  (keyd(l) .ge.2)  print  2,  n,  x0(n),  y0(n),  0.0,  vy0(n) 

20  CONTINUE 

e  Position  northern  tanks 

c  find  center  of  northern  line  of  tanks 

c  Change  by  HLReed  thetal  is  no  longer  used  to  provide  the  cardioid 

c  distribution.  That  is  done  in  the  new  vulnerability  aodel . 

thetal  =  0. 

xcen  =  rg0*s in (thetal) 
ycen  =  rg0*cos(thetal) 

c  place  northern  tanks  on  northern  line 

c  Change  by  HLReed  theta2  is  no  longer  used  to  provide  the  cardioid 

c  distribution.  That  is  done  in  the  new  vulnerability  aodel. 

theta2  =  0. 

headng  =  thetal«theta2»PI 
cosa  =  cos (headng) 
sina  =  sin (headng) 
dx  =  -cosawspacng 
dy  *  sinawspacng 

x0(i north)  *  xcen-0.S«(nnorth-l)*dx 
y0(i north)  *  yeen-0.5*(nnorth-l)»dy 
vy0(i north)  *  0.0 
da  =  headngwdeg 
n  *  i north 

if  (keyd(l) .ge.2)  print  2,  n,  x0(n),  y0(n),  da,  vy0(n) 

CHANGED  31  Uar  85  Following  2  lines  added. 

aspect (i north, TURRET)  =  headng 
aspect (i north, HULL)  =  headng 
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00  30  n=inorth»l, jnorth 
c  Related  to  6  Uar  89  change 

inwatch (n)  =  .false. 
x0(n)  *  x0(n-l)  ♦  dx 
y® (n)  =  yfl(n-l)  ♦  dy 
yy0(n)  =0.0 

aspect (n, TURRET)  =  headng 
aspect  (n , HULL)  =  headng 
dw  =  headngtdeg 

if  (keyd(l).ge.2)  print  2,  n,  x0(n),  y0(n),  dw,  vy0(n) 
30  CONTINUE 

I ast=nnorth*nsouth 
if  (trace)  print  *,,<depio2’ 

END 

c  V7.1 
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SUBROUTINE  DEPLOY 

c  0  Deploy:  position  i  orient  all  tanks  at  beginning  of  engageaent. 
include  ’coaaon.h' 

if  (trace)  print  #,,>deploy’ 

DO  20  n3l,nblu«nred 
t0(n)  =  0.0 
20  CONTINUE 

IF  (scene. eq.BATTAK)  THEN 
c  position  blue  tanks  on  the  x-axis 

cal  I  deplo2(l,nblu1nblu,nblu»l,nblu-nred,nred) 

ELSE 

c  position  red  tanks  on  the  x-axis 

cal  I  deplo2(nblu*l,nblu«nred,nred,l,nblu,nblu) 

ENDIF 

if  (trace)  print  «,‘<deploy’ 

END 

c  V7.1 
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SUBROUTINE  DETECT  (t,  firer,  tgt) 

c  3  Detect:  find  if  tgt  detected  ind  schedule  subsequent  events, 
include  'coaaon.h' 
integer  firer,  tgt,  srayf,  arayt 
1  forast  (f8.2,lx,a4, i3, 1  detects  ',lx,a4,i3) 

c 

if  (trace)  print  «,’>detect’ 
arayf  *  aray (firer) 
arayt  *  3-arayf 

IF  (los(f irer.tgt)  .and.  .not.seen(f irer,tgt)  .and. 

1  ndet(f irer) . It. ndets (arayf))  THEN 

if (keyd(l).ge.2)print  1 , t,co lor (arayf ),f irer, color (arayt), tgt 
.ndet(firer)  a  ndet(f irer)»l 
seen  (f i rer, tgt)  a  .true. 

c  Set  thuaan  to  zero  as  does  GROUNDfARS  -  HLReed 

t  huaan  -  0.0*exp(rol ln(0.5)) 
call  selecs(t,f irer, thuaan) 
c  Change  by  HLReed  to  alio*  transfer  of  targets 

IF(xxfer(areyf))  THEN 
i  =  1 

if(firer  .GT.  nblu)  i  *  nblu»l 
call  skedul (t*2. , i, 'xfer  ’ , tgt) 

ENDIF 
END  IF 

if  (trace)  print  v/Cdetect' 

END 


c  V7  3 
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SUBROUTINE  DET  RC  (naray) 

c  Det  rg :  Find  the  aax  ranges  at  which  each  firer  in  'naray'  detects, 

include  'coaaon.h' 

integer  naray,  first,  last,  tank,  cond,  krg 
real  pi,  p2,  r,  rl,  p 

1  foraat  ('  Range  to  ahich  tank  can  see',/, 

1  *  Tank  HO  FE-S  FE-U  ranu') 

2  foraat  ( i 5 , 3f 8 . 1 , f 8 . 4) 
c 

if  (trace)  print  »,’>detrgl 
if  (keyd(l)  .ge.2)  print  1 

c  Find  first  and  last  firers  on  this  side  (naray). 

IF  (naray.eq.8LU)  THEN 
first  =  1 
last  a  nblu 
ELSE 

first  =  nblu*l 
last  =  nblu»nred 
ENDIF 

c  Loop  thru  all  tanks  on  the  side 

00  80  tank  a  first, last 
p  -  ranu (0.0) 

DO  70  cond=l,3 
pi  =  1.0 

c  Search  for  P- infinity  values  bounding  x 

DO  60  krg=l,8 

p2  =  pinf infnaray, cond, krg) 

IF  (p2  .It.  p)  GOTO  65 

pi  ,  p2 

60  CONTINUE 

p2  =  0.0 

65  CONTINUE 

c  Interpolate  on  p- infinity  to  find  range, 

rl  a  irgincv(krg-l) 
r  *  rl  *  i rg i nc« (pl-p) / (pl-p2) 
rgvis(cond, tank)  a  r 
70  CONTINUE 

if  (keyd(l) .ge.2)  print  2, tank, (rgvis(cond.tank) ,cond=l,3)  ,p 
80  CONTINUE 

if  (trace)  print  «,'<detrg’ 

END 


c  V7.1 
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FUNCTION  DEVIC2  (ittn,  range) 
c  Devie2:  find  resolvable  cycles  for  device  2. 

real  eofs(8),  s(7) 
save  cofs,  fovia,  s,  teepd 
data  teepd,  foviw  /2.S,  .54/ 
data  cofs  /. 41089327,  -.0892577,  .026008138, 

1  -.004097143,  .0003209958,  -.00000982049/ 

data  s/4.798,  .938,  -.236,  .011,  .013,  -.001,  0./ 
c 

c  Calculate  contrast 

y  =  -.45  ♦  1.19*alogl0(attn) 
extcof  =  10.**(y) 
ails  =  cofs  (6) 

DO  30  i =5, 1,-1 
aels  =  aalsvrange  -  cofs(i) 

30  CONTINUE 

attn2  s  extcof-aals 

cntrst  *  abs(teapd*exp(-attn2»range)) 

IF  (cntrst. gt. 0.0112)  THEN 
c  Target/background  is  sufficient  to  detect 

clog  »  a  log (cntrst) 

c  Ire  =  sua  froa  i=l  to  7  {s  sub  i  clog  sup  { i -1}}# 

rc  =  s(7) 

DO  40  i =8 ,1,-1 
re  =  rcvclog  ♦  s(i) 

40  CONTINUE 

devic2  =  aainl(rc/3.44, ,9/foviw) 

ENDIF 

END 

c  V7.1 
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SUBROUTINE  DIS  ENG  (t,  firer,  tgt, drop, take) 
c  7  Diseng:  atteapt  to  disengage  1  firer  froa  1  target, 

e  Oiseng  is  called  by  iapact  if  firer  condition  warrants, 

c  When  I  include  guns,  other  routines  nay  call  it. 

include  ’coaaon.h’ 
integer  arayf,  arayt,  tgt,  firer 
logical  in  brst,  hav  aao,  on  tgt,  drop,  take,  eango 
3  foraat  (f8.2,lx,a4, i3, ’  dis-engs  ’ ,a4, i3,20x, *|tgts=* , i2) 

c 

if  (trace)  print  »,'>diseng’ 
c  Set  useful  local  variables 

ay  tgt  =  nrtgt (firer) 
arayf  *  aray (firer) 
arayt  =  3-arayf 

hav  aao  =  nrd(f irer) . It.nrds(arayf) 

inbrst  ■  nrpb(arayf) ,gt. 1  .and.  (0.ne.aod(nrib(f irer) , 

1  nrpb(arayf))) 

if  (tgt. eq. FIS  TGT)  on  tgt  =  .true, 
if  (tgt.ne.FLS  TGT)  on  tgt  =  fot(f irer, tgt)  .or. 

1  (kindrd (arayf) .eq.4  .and.  eot(f irer, tgt)) 

IF  (on  tgt)  THEN 
c  Firer  on  this  target 

kind  *  kindrd(arayf) 

IF  (kind. I e. 2  .or.  kind.eq.S)  THEN 
IF  (nrpb(arayf) . le.l)  THEN 

c  Single  shot  gun  systea  or  STAFF  fire  i  forget  systea. 

IF  (tgt.ne.FLS  TGT)  THEN 

if  (fot(f irer, tgt))  call  cancel  (firer, ’fire  ’,tgt) 
fot(f irer, tgt)  »  .false. 

ENDIF 

hav  aao  *  nrd(f irer) . It.nrds(arayf) 

IF  (hav  aao)  THEN 
thuaan  *  0.*exp(rol In(B.S)) 
call  selecs(t,f irer.thuaan) 

&SEIF  (can  go(firer,t).and.ishtfs(arayf).gt.B)  THEN 
c  Firer  aoves  on. 

if (keyd(l).ge.2)print  3,  t.color(arayf), firer, 

1  color(arayt),tgt,nchan(f irer) 

call  skedul (t,f irer, 'accel  ' .NULL) 

ENDIF 

n rot (firer)  *  0 
nrtgt  (firer)  *  0 

h^e 

c  Burst  fire  gun  systea. 

print  *,'DISENG:  Not  iapleaented  for  burst  fire  guns.’ 
STOP 
ENDIF 

BLSEIF (kind. eq.4)  THEN 
c  Guided  aissile  systea. 

if  (drop)  nchan(firer)  =  nehan(f i rer)-l 
IF  (tgt.ne.FLS  TGT)  THEN 
IF  (fot(f irer, tgt))  THEN 
call  cancel (firer, 'fire  ’,tgt) 
fot(f irer, tgt)  =  .falsa. 

ENDIF 

ENDIF 

if (keyd(l) .ga.2)print  3,  t, col or (arayf) ,f irar, 

1  color(arayt) , tgt, nehan(f irer) 

c  Firer  atteapts  to  select  a  new  target 

IF  (take)  THEN 

call  frdasl  (t,f irer, tgt, arayf) 

c  The  firer  begins  to  select  a  new  target  right  now  and 

c  finishes  the  selection  in  a  few  seconds. 

ENDIF 

ENDIF 

ENDIF 
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c  V7. 1 


IF  (.not. repeat)  THEN 
repeat  *  .true. 

call  skedul  (t*. 01,0, ’search’ , NULL) 
ENDIF 

if  (trace)  print  «,,<diseng’ 

END 
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FUNCTION  OOT  (a,  b) 

c  9  Dot:  find  dot  product  a  dot  b. 
dimension  a(3),  b(3) 

dot  *  a(l).b(l).a(2)*b(2).a(3).b(3) 

END 


c  V7.3 


SUBROUTINE  ENGAGE  (tl,  t2,  firer,  tgt} 
c  ?  Engage:  begin  engagement  of  a  nee  tgt  by  this  firer. 
include  'coeeon.h' 
integer  areyf,  areyt,  firer,  tgt 
1  foreatC  ENGAGE:  areyf,  ishtfs.fi rer,eotion,STATNY’,8i3) 

c 

if  (trace)  print  «,’>engage’ 
areyf  *  a  ray  (firer) 
areyt  =  3-areyf 

IF  (I ife(f i rer) . lt.FKILL.AND.nrd(f irer) . It.nrds(areyf)  )THEN 
if  (keye(18) ,gt.l)print  1, 

1  areyf, ishtfs (areyf) ,f i rer,eotion(f i rer) .STATNY 
nbrst (firer)  =  1 

IF  (ishtfs(areyf) .gt.0  .AND.  eotion(fi rer) .ne. STATNY 
1  .AND.  speed  (areyf ) .gt. 0.0)  THEN 

c  halt  to  fire 

call  cancel  (f  irer,  ’aaxvel' ,NUJ.) 
call  cancel  (f  irer, 'aceel  ',NU_L) 
cal  I  skedul  (tl, firer,  'sloeup'  ,NlXl) 

ELSE 

c  Schedule  a  fire  event  otherwise 

c  find  range  to  target 

IF  (tgt.eq.-l)  THEN 
rg  =  rg0 

CHANGED  1  Apr  86.  Next  line  changed, 
c  nrg  =  int  ((250 . -rg)» . 002) 

nrg  =  int(l.S-rg/i rginc) 

CHANGED  11  Jun  86  Preceding  line  changed  to  next  line, 
nrg  a  i nt (0 . 5- rg/ i rginc) 

B.SE 

da  =  rgf(tl, tgt, firer) 

ENDIF 

nrg  a  ein0(8,nrg) 

dt  a  tf  irst(arey(firer),nrg)  »  exp(rol  ln(0.5)) 
prev  rd (firer)  a  1 
nrib(firer)  a  0 
nrot(firer)  a  0 

CHANGED  16  Jul  86  Next  line  added, 
c  if (kindrd  (areyf) .eq.4)  dt=0.1 

change  23  Nov  89  by  HLReed  to  aake  sure  a  round  has  been  loaded 
t3  =  aaaxl(tf ire2(f irer)»tein(arayf),t2*dt) 
call  skedul  (t3, firer, 'fire  1 , tgt) 

ENDIF 

ENDIF 

IF  (trace)  print  *,’<engagel 
END 
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SUBROUTINE  EVENTS 

c  9  Events:  call  each  event  in  sequence, 
include  'coeeon.h1 
charactered  iehat 

1  foreat  ('  EVENTS:  No  such  event  type.  Events' ’,a6, 1,1 

1  ’  Who**, i2. *  Iho*=’ , i2, '  Tiae=’,f7.2) 

c 

if  (trace)  print  «,’>events’ 
c  Initial ize  for  battle 

call  reset  (keyd(S) . gt.0) 
call  e reset 
call  init 
ta  1st  =  0.0 

c  Perfore  all  events  in  the  battle 

10  CONTINUE 

call  nextev  (iwho,  iwhat,  iwhoa,  t) 

IF  ( i what. eq. ’search’)  THEN 
call  search  (t) 

ELSEIF  (iwhat. eq. ’vanish’)  THEN 
call  vanish  (t, i who, i whoa) 

ELSEIF  ( i what. oq. ’appear’)  THEN 
call  appear  (t, i who, i whoa) 

ELSEIF  (iwhat. eq. ’detect')  THEN 
call  detect  (t, iwho, iwhoa) 

ELSEIF  ( i what. eq. ’select')  THEN 
call  select  (t.iwho) 

c  Change  by  HLReed  to  allow  transfer  of  targets 

ELSEIF  (iwhat. eq. 'xfer  ')  THEN 
call  xfer(t, iwho, iwhoa) 

ELSEIF  (iwhat. eq. 'fire  ')  THEN 
call  fire  (t, iwho, iwhoa) 

ELSEIF  (iwhat. eq. ’ iapact')  THEN 
call  iapact  (t, iwho) 

0-SEIF  (iwhat. eq. 'slowup')  THEN 
call  slowup  (t, iwho) 

B-SeF  (iwhat. eq.'halt  ’)  THEN 
call  halt  (t, iwho) 

0-SEIF  (iwhat. eq.'accel  ’)  THEN 
call  accelf  (t,iwho) 

0-SEIF  (iwhat. eq. 'aaxvel ’)  THEN 
call  aaxvel  (t.iwho) 

0-SEIF  (iwhat. eq. '  iki 1 1  ’)  THEN 
call  latekl  (t, iwho, iwhoa) 

B-SEIF  (iwhat. eq. ’hide  ')  THEN 
call  hide  (t,iwho) 

0.SEIF  (iwhat. eq. 'reload')  THEN 
call  reload  (t.iwho) 

0-SEIF  (iwhat. eq, ’popdn  ’)  THEN 
cal  I  pop  dn  (t, iwho) 

0-SEIF  (iwhat. eq. 'finish')  THEN 
call  finish  (ta  1st) 

GOTO  99 
0L.SE 

print  1, iwhat,  iwho,  iwhoa,  t 
STOP 
ENDIF 
ta  lst=t 
GOTO  10 

99  if  (trace)  print  «, ’(events’ 

END 


c  V7.1 
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FUNCTION  EYE  (aluenc,  attn,  range,  visrg) 

c  Eye:  find  resolvable  cycles  for  the  huaan  eye.  (Device  1) 

real  a (4, 7) 
save  a,  aeon 

c  sunset  o'cast  heavy  o'cast  overcast  day  clear  day 

data  a / 

1  1.2378091942,  1.7176918034,  1.990992801S,  2.0892716525, 

2  0.4694720809,  .4739034812,  .4484981232,  .2813866389, 

3  .0493317078,  -.2102695514,  -.4084256747,-1.0084578626, 

4  -.0601756751,  -.4161055149,  -.6856409935,-1.4323484287, 

5  -.0558327470,  -.2696921300,  -.4318233767,  -.8450225947, 

6  -.0174190671,  -.0756229822,  -.1197712507,  -.2235482536, 

7  -.0018530403,  -.0077222394,  -.0121729428,  -.0218136690/ 

data  aeon  / . 4/ 

c 

c  Find  sky-to-ground  ratio 

sog  =  (visrg*1.0)/3. 
sog  »  aainl(3. , a«axl(l. ,sog)) 
entrst  =  aeon/(1.0-sog«(exp(attn«range)-1.0)) 

IF  (entrst. ge. 0.02)  THEN 

c  Target/Background  contrast  is  sufficient  to  detect 

i  =  ■in0(4,l*int(alogl0(a!u«ne))) 
ack  a  10.**i 
j  ■  ain0(4, i*l) 
clog  =  a  log (entrst) 
rlo  =  a(i  ,7) 
rhi  =  a  (j ,  7) 

DO  20  k=6, 1,-1 

rlo  =  rlo  «  clog  ♦  a(i,k) 

rhi  =  rhi  «  clog  ♦  a(j,k) 

20  CONTINUE 

c  Interpolate  1  coepute  cycles  across  target 

eye  =  rlo-(rhi-r(o)«(alu«nc-ack/10.)/(ack*.9) 

ENDIF 

END 

c  V7.2 
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SUBROUTINE  FINISH  (t) 

c  3  Finish:  update  statistics  at  and  of  a  single  engagement, 
include  'coaaon.h' 

integer  balive,  dal  ire,  ralive,  brds,  rrds 
di sens  ion  state  (8) 

1  foraat(i6,2(Si3) ,4i3,lx,2f5.1, a  9) 

2  foreat(’  Rep  Status  of  Coebatants  ', 

I  '  Rds  Used  Used/Tank 1 ,/ 

1  '  | - Blue— -|  | . Red— -|  *, 

1  'by  System  Blue  Red  seed’,/ 

1  6x,2(lx,'AL  UO  FO  kF  K’),2x,*l  2  3  4’) 
c 

if  (trace)  print  «,'>finish' 
c  Count  surviving  blues  and  rounds  fired 

balive  =  0 
brds  -  0 
dal ive  -  0 

c  Change  aade  by  H.L.  Reed  on  March  31,  1989  to  keep  decoys  out  of 
c  win  ratios  and  exchange  ratio.  See  also  deaths. f 
DO  10  i=l, (nblu-ndecoy (BLU)) 
k  *  life(i) 
if  (k.ge.S)  k=k-l 
nstats(k.BLU)  =  nstats(k,8LU)*l 
if  (I ife(i) . It.FKILL)  balive  =  balive*! 
brds  *  brds*nrd(i) 

10  CONTINUE 

DO  11  i  =  nblu  -  ndecoy(BLU)  *1,  nblu 
if  (I ife(i). It.FKILL)  dal ive  =  dalive.l 

11  CONTINUE 
call  finsh2 

c  Count  surviving  reds  and  red  rounds  fired, 

ralive  *  0 
rrds  =  0 

DO  20  i*l,  (nred-ndecoy (RED)) 
j  *  i*nblu 
k  =  I ife(i»nblu) 
if  (k.ge.5)  k=k-l 
nstats(k,RED)  a  nstats(k,RED)*l 
if  (I ife(j) . It.FKILL)  ralive  =  ralive*l 
rrds  *  rrds*nrd(j) 

20  CONTINUE 

■ystat(l)  a  aystat(l)  *  nblu-balive  -  ndecoy(BLU) 

aystat(2)  a  aystat(2)  *  ndecoy(BLU)  •  dative 

naov  =  0 
ndow  a  0 

DO  23,  i  a  nblu»l,  nblu-nred  -  ndeeoy(RED) 
if (I ife(i) .ge.FKILL.and. inwatch(i))  ndow  a  ndow  ♦  1 
if (I ife(i) .ge.FKlLL.and. .not. inwatch(i))  naov  a  naov  *1 
23  CONTINUE 

aystat(3)  a  aystat(3)  *  naov 

aystat(4)  a  aystat(4)  ♦  ndow 

c 

00  30  i *1 , 4 
30  stata(i)  =  0.0 

j  a  nred-nwatch(scene)-ndecoy(RED)-naov 
if  (balive.gt.0  .and.  j.eq.0)  stata(l)al. 

if  (balive.eq.0  .and.  j.gt.0)  stata(2)=l. 

if  (balive.gt.0  .and.  j.gt.0)  stata(3)al. 

if  (balive.eq.0  .and.  j.eq.0)  stata(4)ai. 

stata(S)  a  (nblu-ndecoy (BLU) )-ba I ive 
stata(S)  a  (nred-ndecoy (RED))-ra I ive 
stata(7)  a  float (brds) /float (nblu) 
stata(8)  a  f loat(rrds)/f loat(nred) 
excha  a  0.0 

if  (stata(S) .gt.0.0)  excha  a  stata(8)/stata(5) 
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00  40  1=1,8 

statb(i)  =  statb(i)«stata(i) 

40  CONTINUE 

c  Update  aeao  consuaption  for  F-alive  Blue  tanks  if  eoabat  occurred. 

IF  (brds»rrds  .gt.0)  THEN 
DO  50  i=l,nblu 
j=nrd(i)»l 
50  CONTINUE 

ENDIF 
c 

if  (keyd(l) .ge.2  .or. 

1  (krep.eq.l  .and.  keyd('l) .eq.l))  print  2 

if  (keyd(l) .gt.0)  print  1,  krep,  nstats,  (nrd(i) , i =1 , 4) , 

1  stata (7) ,  stata(S),  iranda 

if  (trace)  print  e.’Cfinish’ 

END 

c  V7.1 
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c  9 
c 

c 

c 

90 

c  V7.1 


SUBROUTINE  FINSH2 

Finsh2:  update  statistics  at  «nd  of  a  single  engageeent. 
include  'coeeon.h' 

if  (trace)  print  *,,>finsh2' 

DO  90  i -1 , nb I u 

select  blues  for  further  coabat 
IF  (I ife(i) .eq. ALIVE  .and.  nwa*es.gt.l)  THEN 
nsurv  -  nsury*l 
nused(nsurv)  =  nrd(i) 

ENDIF 

IF  (nrd(i-)  .gt.nrds(BLU)-S)  THEN 
count  systeas  with  no  i  low  aaao 

if  (nrd(i) . It.nrds(BLU))  loaaao  =  loaiio»l 
if  (nrd(i) .ge.nrds(BLU))  noaaao  =  noano*l 
ENDIF 
CONTINUE 

if  (trace)  print  *1’<finsh2' 

END 
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SUBROUTINE  FIRE  (t,f irer.tgt) 

c  7  Fire:  Simulate  firing  of  a  round  k  schedule  effects, 
include  ’coeaon.h1 

integer  bullet,  areyf,  arayt,  firer,  tgt 

1  forsat(f8.2,  lx,  a4,  i3,  ’  fires  at  1 ,  a4,  i3) 

2  foraat(f8.2,  lx,  a4,  i3,  ’  ran  out  of  aeso.’) 

c 

if  (trace)  print  •,'>fire’ 
busy(firer)=. false. 

IF  (life(firer).ge.FKILL)  THEN 
print  «,’FIR£:  f i rer' , f i rer, 1  is  F-killed  or  worse.’ 

STOP 

ELSEIF  (tgt.eq.0)  THEN 

print  «,’FIRE:  f irer’, firer, '  has  no  target.’ 

STOP 

ELSE 

c  Find  nrs  for  tgt,  arsy  of  firer,  arty  of  tgt 

areyf  =  arey (firer) 
areyt  =  3-arsyf 

if  (keyd(l).ge.?)  print  l,t,color(areyf), firer, 

1  eo I  or (arayt) ,tgt 

c  Update  last  firing  tiae  for  firer  i  for  firer  at  this  tgt 

if  (tgt.gt.0)  tfi re (firer, tgt)  =  t 
tfire2 (firer)  =  t 

c  Update  positions,  velocities  and  turret  orientation 

IF  (tgt.eq.-l)  THEN 
rg  =  rg0 

nrg  *  aax0(l,  int(0.5-rg/irginc)) 
s(l)  =  0.0 
s(2)  *  rg0 
s(3)  »  0.0 

if  ((arayf .eq.BLU  .and.  scene. eq.BATTAK)  .or. 

1  (arayf .eq. RED  .and.  scene . ne . BATTAK) )  s(2)  *  -rg0 

B_SE 

da  »  rgf(t, tgt, firer) 

ENDIF 

aspect (f irer .TURRET)  =  anglef (VNORTH.s) 
c  Schedule  any  pinpoint  detections 

call  pinpnt  (t, firer) 

IF  (if lash(f irer) .eq.0)  THEN 

c  Branch  for  real  firer  (do  nothing  if  firer  is  flashing  decoy) 

c  Create  round  with  various  attributes 

call  create  (10,bullet) 
a(bul let-1)  =  tgt 
a(bul let-2)  =  firer 
tfly  =  tof (arayf , nrg) 
t2  =  t-tfly 

a (bu I  let-3)  =  s(l)-tf lywvt(l) 
a  (bullet-4)  *  s(2)-tf ly«vt(2) 
a  (bullet-7)  =  psense (arayf, nrg) 
a  (bullet-9)  *  vabs(vf) 
if  (tgt.eq.-l)  a(bullet-10)  =0.0 
if  (tgt.gt.0)  a(bul let-10)  ■  vahs(vt) 
a (bu I let-10)=  vabs(vt) 
kshot(arayf ,1)  *  kshot(arayf ,1)  -  1 

c  Schedule  iapact  for  bullet 


call  skedul  (t-tfly, bul let, 'iapact', tgt) 

IF  (kind  rd(arsyf) .eq.4)  THEN 
if  (tgt.gt.0)  aot(f irer.tgt)  * 

.true. 

00  20  i>l,S 

IF  (asl  f ly(arayf, firer, i)  .eq. 

0)  GOTO  25 

20 

CONTINUE 

25 

print  «,'FIRE:  Too  aany  aissiles' 
STOP 

asl  f ly(arayf, firer, i)  *  bullet 

ENDIF 
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e  Update  stowed  rounds  and  expenditure 

nrd(firer)  ■  nrd(firer)*l 
nrib(firer)  *  nr i b(f i rer) -*1 
if (nrib(f irer) .gt.nrpb(areyf))  nrib(firer)*l 
nrot(firer)  »  nrot(f irer)«l 
e  llove,  fire,  or  switch  targets  as  required 

IF  (kind  rd(arsyf) . Ie.2  .or.  kind  rd(arsyf) .eq.S)  THEN 
IF  (nrpb(arsyf)  .le.l)  call  frd  ssg(t,f irer.tgt.areyf) 
ELSEIF  (kind  rd(areyf)  ,eq.  4)  THEN 
c  Sieultaneous  aissiles  branch 

IF  (nehan(f irer). It.nchans(arsyf))  THEN 
call  frd  esl (t,f irer,tgt,areyf) 
c  ELSE 

c  All  guidance  channels  busy,  fait  until  iipact. 

ENDIF 

B.SE 

print  «,’FIRE:  kind  rd’.kind  rd(firer),’  unknown.* 
ENDIF 
ENDIF 
ENDIF 

if  (keyd(l) .ge.2  .and.  nrd(f irer).ge.nrds(arsyf))  print  2, 

1  t, color (arwyf) ,f i rer 
if  (trace)  print  *,*<fire‘ 

END 
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SUBROUTINE  FORCES 

c  7  Forces:  loop  through  desired  blue/red  force  ratios, 
include  'coaeon.h' 
integer  ranged 

1  formate  SCENff: ' ,  i3) 

2  foraat('  Meeting  engagement.  fBlues  ■ ’ ,  i 3 , '  (Reds  =  ’ , i 3) 

3  formate  N«d  attack.  (Blues  = 1 ,  i 3 ,  *  (Reds  «*,i3) 

4  formate  Slue  attack.  (Blues  ** ,  i3,  *  (Reds  =  ’ ,  i 3) 

c 

if  (trace)  print  »,'>forces‘ 
min  blu  *  ntanks (scene, 1) 
min  red  =  ntanks (scene, 4) 

IF  (min  blu  ,gt.  0  .and.  min  red  .gt.  0)  THEN 
max  blu  *  ntanks (scene ,2) 
max  red  *  ntanks (scene, S) 
inc  blu  =  max0(l,ntanks(scene,3)) 
inc  red  =  max0(l, ntanks (scene, 8)) 

DO  50  nblu  3  sin  blu, max  blu, inc  blu 
DO  40  nred  3  min  red, max  red, inc  red 

if  (scene  .eq.  UEETNG)  print  2,  nblu,  nred 

if  (scene  .eq.  RATTAN)  print  3,  nblu,  nred 

if  (scene  .eq.  8ATTAK)  print  4,  nblu,  nred 

DO  30  ranged  *  min  rg,  max  rg,  inc  rg 

rgd  3  ranged 
ca 1 1  saves  (scene) 

30  CONTINUE 

40  CONTINUE 

50  CONTINUE 

ENDIF 

>f  (trace)  print  (forces’ 

END 
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SUBROUTINE  FRD  USL  (t,  firer,  tgt,  arayf) 

c  0  Frd  asl:  Fired  a  aissile.  nos  schedule  effects, 
include  'coaaon.h' 
logical  done,  tactc3,  prflag 
integer  arayf,  firer,  tgt 

1  foraat('FRD  USL:  t, firer, tgt, arayf*', f7. 2, 3i3) 

2  foraat(f8.2,  1*.  a4,  i3,  1  begins  to  reload.') 

3  foraat(’FRD  USL:  f irer=',i3, ’  tactic*' , i2, '  |rds  fired** 

1  *  frds  to  f ire*', i2) 

c 

if  (trace)  print  »,’>frd  asl* 
prflag*. false. 

if  (prflag)  print  1 , t , f i rer, tgt, arayf 
IF  (nrd(f irer) . It.nrds(arayf))  THEN 

c  Systea  has  aore  rounds  on  board. 

if  (prflag)  print  «,  'FRD  USL:  nrd(f irer)*’, nrd(f irer) 
IF  (aod(nrd(f irer) .nipods(arayf)) .gt.0  .or. 

1  nrd(f  irer)  .eq.0)  THEN 

c  Systea  has  aore  rounds  in  pod. 

if  (prflag)  print*,'  FRD  USL:  No  reload,  nrd,  n i pods 
1  nrd(f  i  rer)  ,nipods (arayf) 

tactc3  *  tactic(arayf) .eq.3 
done  =  nrot(f i rer) .eq.nrpt(areyf) 
if  (prflag)  print  3,  firer,  tact ic (arayf) , 

1  nrot (firer),  nrpt(areyf) 

IF  (taetc3  .and.  done)  THEN 

c  Seitch  targets  after  firing  a  fixed  nr  of  rds  at  it 

if  (tgt.ne.FLS  TGT)  fot(f i rer, tgt)  =  .false, 
call  selecs  (t, firer, 0.0) 

ELSE 

c  Schedule  next  round  fired 

tiaea  *  tain(arayf) 

if  (prflag)  print  «,  'FRD  SSG:  shoot  again.' 

tiaeb  *  tf ixed(arayf ,nrg) 

tiaec  *  tied  in (arayf)  *  exp(rol I n (0 . S) ) 

dt  *  aaaxl(tiaea,tiaeb«tiaec) 

call  skedul  (t-dt,f irer, 'fire  ’ , tgt) 

'  ENDIF 

ELSE 

c  Treat  eapty  aissile  pod 

if  (prflag)  print  FRD  USL:  Reloading’ 
eapty (firer)  *  .true, 
call  cancel (firer, ’fire  ',tgt) 
call  cancel  (fi rer, 'select' .NULL) 
nrot(firer)  *  0 

c  shud  htf  that  is  slowing  to  engage  speed  up  nov? 

dill  skedul  (t*trelod(arayf)  ,fi  rer ,  'reload' ,  NLRJL) 
if  (keyd(l) .ge.2)  print  2,t,color(arayf),f irer 
ENDIF 
ENDIF 

if  (tgt.gt.0)  fot(f i rer, tgt)  *  .false. 

C  ABOVE  LINE  GOOD  FOR  UUSL  THAT  IS  NOT  LOADED  »/  TARGETS 

if  (trace)  print  *,’<frd  asl’ 

END 
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SUBROUTINE  FRO  SSG  (t,  firer,  tgt,  areyf) 
c  6  Frd  ssg:  Schedule  effects  after  firing  single  shot  gun. 
include  ’coaeon.h’ 
logical  can  go,  done,  tactc3 
integer  areyf,  firer,  tgt 

1  foreat('FR0  SSG:  t,f irer,tgt,areyf=’,f7.2,3i3) 

2  foreat(f8.2,  lx,  a4,  i 3 ,  'is  out  of  aiao.  Vi  1 1  atteapt’, 

1  ’  to  hide  if  aobi ie. ' j 

c 

if  (trace)  print  «,'>frd  ssg' 

IF  (nrd(f  irer) .  It.nrds(areyf))  THEN 
c  Have  aaao  branch 

tactc3  =  tactic(arayf) .eq.3 
done  =  nrot(f irer) .eq.nrpt (areyf) 

IF  ((tactc3  .and.  done))  THEN 

c  Switch  targets  after  firing  a  fixed  nr  of  rds  at  it 

busy (firer)  =  .false. 

call  dis  eng  (t,  firer,  tgt, .true. , .true.) 
c  If  no  other  tgt  and  can  aove,  skedul  acceleration 

if  (can  go (firer, t)  .and.  ishtfs(arayf) .eq. 1) 

1  call  skedul  (t, firer, 'acce I  ’.NULL) 

nrot  (firer)  =  0 
ELSEIF  (tgt.gt.0)  THEN 
c  Schedule  next  round  fired 

tiaea  =  tain(arayf) 
tiaeb  =  tf ixed (areyf ,nrg) 
tiaec  =  taedin(arayf)  •  exp(rol I n (0 . 5) ) 
dt  *  aaaxl(tiaea,tiaeb-tiaec) 
call  skedul  (t*dt,f i rer, *f i re  ’ , tgt) 

ENDIF 

ELSE 

c  Out-of-aeeo  branch 

eapty (firer)  a  .true. 

IF  (cango(f irer.t))  THEN 
call  skedul  (t,f irer, 'accel  ‘.NULL) 
call  skedul  (t-thide (areyf) ,f i rer, 'h ide  ’.NULL) 
ENDIF 
ENDIF 

if  (trace)  print  «,'<frd  ssg' 

END 
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SUBROUTINE  HALT  (t,  firer) 
c  7  Halt:  siaulate  tank  halting, 
include  'coaaon.h' 
logical  cango,  threat 
integer  arayf,  firer,  tgt 

1  foraat  (f8.2,lx,a4, i3, '  halts',12x, ' (x=' ,f8. 1, '  yx’.fB.l',)’) 

2  foraat(’  HALT:  firer,  Oturrer,  Ohull  =  ',  i3,  2f*8 . 1) 

3  foreat('  HALT:  firer,  tgt,  arayf,  nrg  =',  4i3) 

4  foraat('  HALT:  rxl,rxe,tf irst.dt  =',5fl0.3) 

5  foraatC  HALT:  t,  tlastx,  dt  =* ,5f  10.3) 
c 

if  (trace)  print  * , 1 >ha I t 1 

if (invisb.eq.l)cal I  cancel  (firer, 'vanish' , NULL) 
naray  =  aray  (firer) 

call  path  (f irer, t,aotion(f irer), 0.0, x,y,vx,vy) 

if  (keyd(l) .ge.2)  print  1,  t,  color(naray) ,  firer,  x,  y 

■ot ion  (firer)  =  STATNY 

tlastx  =  t-3. 

arayf  x  naray 

c  see  if  fire  is  a  halt-tc-f ire-systea  and  can  still  shoot 

IF  (ishtfs (arayf) .eq.l  .and. 

1  I i f e (f irer) . It.FKILL  .AND.  nrd(f i rer) . It. nrds (arayf))  THEN 
c  This  is  a  halt-to-fire  systea.  schedule  firing  if  tgt 

c  still  available. 

threat  *  .false. 

IF  (nrtgt(firer).eq.FLS  TGT)  THEN 
threat  =  kneeal (f irer) .ne.FD 
B.SEIF  (n rtgt (firer).gt.B)  THEN 
threat  *  fot(f irer, nrtgt(f irer)) 

ENDIF 

IF  (.not. threat)  THEN 

c  firer’s  tgt  has  vanished,  firer  aay  aove 

if (cango(f irer,t))cal I  skedul  (t,  firer,  'accel  ’,  NULL) 

b_se 

if  (keyd(l) .ge.2)  print  e.'HALT:  tlastx,  aspect  needs  vk!' 
rxl=rol In (0. 5) 
rx2=exp(rxl) 
tgt=n rtgt (firer) 
duaay  *  rgf (t,f i rer, tgt) 
dt  x  tf irst (arayf ,nrg)«rx2 
e  change  Dec  89  by  HLReed 

dt  =  aaaxl(dt,tf ire2(f irer)-tain(arayf)  -t) 
prev  rd (firer)  *  1 
nrib(firer)  »  8 
nrot(firer)  =  0 

call  skedul  (t-dt,  firer,  'fire  1 ,  tgt) 

ENDIF 

ENDIF 

if  (trace)  print  •,’<halt' 

END 
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SUBROUTINE  HIDE  (t,  tgt) 
c  S  Hide:  Simulate  tank  hiding, 
include  ’coaaon.h’ 
integer  firer,  tgt 

1  foraat  (f8.2,x,a4, i3,  ’  goes  into  full  defilade.') 

e 

if  (trace)  print  », 1 >h i de  ' 

if  (keyd(l) .gt.l)  print  1,  t,  eolor(aray(tgt)) ,  tgt 
kneeal(tgt)  *  FD 

c  Cancel  all  activities  involving  this  tgt 

c  except  discard  rounds-in-f I ight  in  the  iapact  routine 

firer  =  1 

if  (tgt. le.nblu)  firer=nblu*l 
last  *  nblu 

if  (tgt. le.nblu)  last=nblu*nred 
00  20  i*firer,last 
los(i,tgt)  =  .false. 

los(tgt.i)  -  .false. 

20  CONTINUE 

call  newtgt  (t,  firer,  tgt) 
call  cancel  (tgt, 'all  '.NULL) 
call  skedul (t, tgt, 'slowup ', NULL) 
call  deaths (t) 

if  (trace)  print  *,'<hide  ' 

END 
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SUBROUTINE  IMPACT  (t,  bullet) 
e  0  lepact:  find  ahat  bullet  does  ft  ahat  firer  does, 
include  ’coaeon. h' 
logical  loaded,  hit 
integer  bullet,  expose 
atr(i)  *  a(bul let«i) 

if  (trace)  print  *,  ’>iapact‘ 
c  Find  useful  variables, 

it  =  atr(l) 

I  =  atr(2) 
n  *  aray(I) 
k  x  kindrd(n) 
expose  x  knceal(it) 
rgx  x  0.0 

c  Find  ahat  bullet  does. 

IF  (it.eq.FLS  TGT)  THEN 
c  Round  does  nothing. 

kshot(n,4)  =  kshot(n,4)«l 
H.SEIF  (expose. eg. FD  .and.  k.le.2)  THEN 
c  Count  round  hitting  bera. 

kshot(n,5)  =  kshot(n,5)-l 

if  (keyd(l) .ge.2)  print  «,  ’Tgt  in  full  defilade.' 

ELSE 

c  See  if  round  hits. 

call  aayhit(t,I, it,n,k,atr(9),atr(10) ,expose,hit) 

END  IF 

a(bullet)  x  -a(bullet) 
c  Find  ahat  firer  does. 

IF  (k.eq.4)  THEN 
c  Missile 

c  Clear  guidance  channel . 

DO  20  jxl.S 

IF  (aslf ly(n, I, j).eq. bullet)  GOTO  30 
20  CONTINUE 

print  *,  'IMPACT:  Msl  not  assigned  a  channel.' 
print  «,  'Channels  assigned  to’ , (aslf ly (n,I, j) , jsl( 5) 
print  «,  'Msl  fx^bullet,'  Contact  Fred  Bunn' 

STOP 

30  CONTINUE 

aslf ly (n, I, j)  x  NULL 
loaded  x  nchan(I) .ge.nchans(n) 
call  diseng  (t, I, it, .true. , loaded) 
aot(I, it)*. false, 
fotfl, it)*. false. 

if  (knceal (I) .eq.HD  .and.  nehan(I) .eq.0  .and.  eapty(I)) 

1  call  skedul (t,I, 'popdn  '.NULL) 

c  KE,  HEAT,  or  STAFF  [rethink  this  for  STAFF) 

IF  (it.eq.FLS  TGT  .or.  hit.and. tactic(n) .eq.2  .or. 

1  rgx. gt. 4000.0)  THEN 

c  Saitch  targets  if  false  target  or  rd  hit  ft  I  saitch  on  a  hit. 

c  Von't  go  here  if  I  hit  the  bera;  fls  tgts  don’t  go  behind  the 

c  bera,  and  if  true  tgts  do,  the  rd  aon't  hit. 

ndet(I)  =  ndet(I)-l 
nrtgt(I)*0 

cal  I  diseng(t, I, it, .true. , .true.) 

ENDIF 

ENDIF 

if  (trace)  print  «,  '<iapact' 

DID 
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INTEGER  FUNCTION  INOEXX(a,  n,  x) 


e  Find  the  index  j,  ehere  a(j)  <=  x  <  a(j*l) 

c  Adapted  froe  Nueerical  Recipes,  p90.  The  array  ii  eust  be 

integer  n  , j I ,  ju,  je 
,  logical  incres,  above 
real  a(n),  x 
c 

incres  =  a  (n) .gt.a(l) 

j  1=0 

ju=n»l 

10  IF  (ju-jl.gt.l)  THEN 
ja=(ju-jl)/2 
above=x.gt.a(je) 

IF  ((incres. and. above)  .or.  .not. (incres. or. above))  THEN 

j  l*j» 

H.SE 

ju=j» 

ENDIF 
GOTO  10 
ENOIF 
indexx=j I 
END 


increasing. 
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SUBROUTINE  INIT 

c  4  Init:  Initialize  scenario  t  schedule  search  at  tiae  zero, 
include  ’coeaon.h’ 
integer  firer,  tgt 
logical  regard 
coeaon  /cdetrg/  tdet(2) 
coaeon  /cregrd /  regard(NN) 
c 

if  (trace)  print  * , 1 >In i t ’ 

call  skedul (taax,0, ’f inish'.NULL) 

call  deploy 

last  =  nred«nblu 

cal  I  in i t2  (1,  nblu) 

call  init2  (nblu*l,  last) 

da  a  rgf  (0.0,l,l»nblu) 

c  Set  state  variables  for  both  red  and  blue  systeas 

DO  30  firer=l,last 
busy (firer)  =  .false. 
eapty(firer)  a  .false. 
serehg(f irer)  a  .true. 
ndet(f  i  rer)  a  0 
naray  a  aray (firer) 
ncol  =  kneeal (f irer)-l 

if  (naray  . eq.BLU  .and.  scene . eq . BA TTAK)  ncol  =  3 
if (naray. eq. RED  .and.  scene. eq.RA TTAK)  ncol  =  3 
cansee(f irer)  a  rg. It. rgvis(ncol , firer) 
regard (firer)  a  .true. 

DO  20  tgt=l , I ast 

foes(f irer,tgt)=  aray(f irer) .ne.aray(tgt) 
know(f irer, tgt)  =  0 

I os (firer, tgt)  a  foes (firer, tgt)  .and.  invisb.ne.2 
aot (firer, tgt)  *  .false. 
fot(f irer, tgt)  »  .false. 
seen(f irer, tgt)  a  .false. 

20  CONTINUE 

30  CONTINUE 

c  Hardwired  values 

accaax(BLU)  a  2.5 
accaax(RED)  a  2.5 
wvlth(BLU)  =  50. 
wvlth(RED)  a  50. 

aapl(BUJ)  =  accaax(BLU)/(TI0PI»speed(BLU)/wvlth(BLU))»»2 
aapl(RED)  =  aceaax(RED)/(T*0PI«speed(RED)/wvlth(RED))»*2 
call  serchl 

if  (trace)  print  •, '<Init' 

END 
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SUBROUTINE  INIT2  (hirst,  last) 
c  0  I n i 1 2 :  initialize  jach  tank  on  one  side, 

include  * coaeon. h 
diaension  iexposv2,3) 
data  iexpos  /3,  3,  2,  3,  3,  2/ 

1  foraat('  INIT2:  nexal,  nrd(l-3)=' ,4i5) 

c 

if  (trace)  print  «,'>init2’ 
naray  =  BLU 

if  (if irst.gt.l)  naray=RED 

Iast2  =  nblu*nred 

jscene  =  iexpos(naray,  scene) 

DO  10  i= if irst,  last 
aray(i)  =  naray 
life(i)  =  ALIVE 
nrd(i)  =  0 
nrtgt(i)  =  0 
nchan(i)  *  0 
nrot(i)  =  0 
knceal(i)  *  jscene 

c  Change  introduced  by  HL  Reed  8  liar  89  to  alloa  oxerwatch  tanks  to 

c  be  added  to  the  attacking  force.  See  also  subroutines  d'<lo2, 

c  Input,  and  cango  and  coaaon.h. 

if  (inwatch (i ) )  knceal(i)  =  HO 
iehg(i)  =  0 
aotion(i)  =  UAXVL 

if  (knceal (i) .eq.HD  or.  scene. eq.MEETNG)  aotion(i)  =  STATNY 
c  End  of  8  Uar  89  changes. 

nhot(i)  s  0 
DO  8  j=l,S 

asl  f I y (naray, i,j)  =  0 
nstats(j ,  naray)  =  0 
8  CONTINUE 

DO  10  j*l,last2 
tfire(i,j)  *  0.0 

c  Change  by  HLReed  to  aake  sure  first  round  is  loaded  see  ENGAGE  also 

tfire2(i)  *  -  tain(naray) 
knee  Is  (i , j)  *  .false. 

10  CONTINUE 

IF  (ndecoy(naray) .gt.0)  THEN 

c  Change  by  H.L.  Reed  on  18  March  89  to  correct  decoys  for  Red  Force 

c  and  to  aake  flashing  decoys  usually  be  overvatching  tanks. 

Idecoy  =  last  -  ndecoy(naray)  ♦  1 
DO  20  i  =  Idecoy,  last 
iflash(i)  =  0 

if  (i.gt.last  -  nf I  ash  (naray))  i f I  ash  ( i )  =  1 
if (if lash(i) .eq.0)nrd(i)=999 
c  End  of  16  March  89  change. 

20  CONTINUE 

ENDIF 

call  detrg(naray) 

IF  (invisb.eq.l)  THEN 
if  (naray. eq. RED  and.  scene . eq . RATTAN) 

1  call  terain  (if irst, last) 

if  (naray. eq. BLU  and.  scene. eq.BATTAK) 

1  call  terain  (if irst, last) 

ELSE 

if  (naray .eq. BLU)  call  saoke 
ENDIF 

c  Correct  the  nr  of  rounds  used  by  blue  systeas 

IF  (naray. eq. BLU  .and.  naaves.gt.l)  THEN 
DO  40  i>l,nblu 
nrd(i)  =  nused(nevaUi) 

<0  CONTINUE 

c  print  1,  neval , (nrd(i), i =1 , 3) 

neval  >  nevaUnblu 
ENDIF 

if  (trace)  print  • , ’ < i n i t2 ’ 

END 
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SUBROUTINE  INPUT 
c  9  Input:  read  s  i  sc  inputs 
include  'coeeon.h' 
character*32  fnase 
integer  indx(S) 

1  forest (il,a32} 

4  forest (a32) 

c 

c  Change  introduced  by  H.L.Reed  on  8  Uar  89  to  allow  overwatch 

c  tanks  to  be  added  to  the  attacking  force.  See  also  changes  in 

c  subroutines  deplo2,  init2,  and  cango  and  in  coeeon.h. 

read  (5, •  )  (ntanks(l, j) , j  =  l,6) ,  nwatch(l) 
read(5,*)  (ntanks(2, j) , j=l,6) ,  nwatch(2) 
read (5 , «) (ntanks (3 , j ) , j =1 , 8) ,  nwatch (3) 
c  End  of  8  Uar  89  change 

read(S,»)(keyd(i),i=l,S) 
trace=keyd(4) .gt.0 
read (5,*) indx 
DO  20  i=l,5 

if  (indx(i) .gt.l  .and.  indx(i) . la. 20)  keya(indx(i))=l 
20  CONTINUE 

read(S,*)ein  rg,  aax  rg,  inc  rg,  irginc 
rgincr  =  irginc 

rginc2  =  0.5* irginc 

read(5,*)  nreps,  nwaves,  iangd,  eeth  se,  irande 

read (5,*)  teax 

read (5,4)  fnaee 

call  rdei sc (fnaee ,BLU) 

c  Read  pkh  data  for  Blue.  Change  by  HLReed  for  new  vulnerability  node  I 

read  1,  i pkh,  fnaee 
call  rdpkh (fnaee ,BLU) 
read(S,4)  fnaee 
call  rdei sc (fnaee, RED) 

c  Read  pkh  data  for  Red.  Change  by  HLReed  for  new  vulnerability  eodel 

read  1,  i pkh,  fnaee 
call  rdpkh (fnaee, RED) 
read(5,«)  invisb.n 
IF  (invisb.ne.l)  THEN 
print  *,'  Seoke  causes  intervisibility.' 
read  • 

read  *,  ((touti l(i , j) , j=l, 5) , i=l,21) 

read  • 

read  *,  ((toutvl(i , j) , j=l, 5) , i*l,21) 

read  « 

read  «,  ((touti (i , j) , j«l,5) , i*l,21) 

read  « 

read  »,  ((toutv(i , j) , j*l,S) , i*l,21) 
read  • 

read  «,  ( (t ini (i , j) , j=l,8) , i=l,21) 

read  • 

read  «,  ( (t i nv ( i , j) , j=l , 5) , i=l,21) 

ENDIF 

if  (trace)  print  *,'< input’ 

END 
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FUNCTION  IZHIT  (nbox,  ndia,  naray,  x,  y,tgt,  theta) 
c  6  Iz  hit:  find  if  the  target  is  hit. 
include  ’ coaeon. h' 
integer  tgt 

1  foraat  ('  IZHIT:  the  round  is  high,  y,  ylieit,  x  a’,3f7.3) 

2  foraat  (’  IZHIT:  the  round  is  low.  y,  yliait,  x  =’,3f7.3) 

3  foraat  (’  IZHIT:  the  round  is  wide,  y,  yliait  =’,  2f7.3,./ 

1  1  x,  xleft,  xright  =  3f7.3) 

4  foraat  ('  IZHIT:  the  round  hits.  y,  yliait  =’,  2f 7 . 3 , / 

1  1  x,  xleft,  xright  =  3f7.3) 

c 

if  (trace)  print  • , *> i zh i t * 
izhit  a  0 

yliait  a  sysdia (naray ,nd ia) 

IF  (yliait. Ie.abs(y))  THEN 
c  Too  high  or  too  low 

IF  (  keya (6) .gt.0)  THEN 

if  (y.gt.0.0)  print  1,  y,  yliait,  x 

if  (y.le.0.0)  print  2,  y,  yliait,  x 

ENDIF 
ELSE 

c  Height  ok 

c  Find  theta 

vtgt  a  vabs(vt) 
theta  a  aspect (tgt, nbox) tdeg 
c  Select  tgt  orientation  between  0  1  360  deg. 

if(vtgt.gt.0  .and.  nbox.eq.HULL)theta=anglef (VNQRTH,vt)*deg 
phi  a  anglef (VNORTH,s)»deg 
thetal  =  angsua(phi , -theta) 
theta  a  thetal/deg 

call  bounds  (naray,  nbox,  theta,  xleft,  xright) 
if  (x.gt. xleft  .and.  x. It. xright)  izhit  a  I 
IF  (keys (0) .gt.0)  THEN 

if (izhit. eq.0)  print  3,  y,y I iait.x, xleft, xright 
if  (izhit. eq.l)  print  4,  y,yl iait,x, xleft, xright 
ENDIF 
ENDIF 

if  (trace)  print  »,'<izhit’ 

END 
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SUBROUTINE  KILL  (firer,  tgt,  hit,  injury, r) 
c  9  Kill:  find  kill  type  for  a  hit  on  a  tgt. 
c  The  routine  is  called  by  aayhit. 

c  changed  May  18,1989  for  simplified  hit  and  kill  aodel,  HL  Reed 
include  ’eosaon.h1 
logical  hit 
integer  firer,  tgt 
coaaon  /cpkh2/  pki 1 1 (2,3,2,S,9) 
save  /cpkh2/ 

c  Change  for  interpolation  on  range  for  pki 1 1 

p(i)  =  (1.0  -  r)  «  pki 1 1  (narsy,ncase,nhdfe, i,jrg) 
l  ♦  r  «  pki 1 1  (naray,ncase,nhdfe, i , jrg*l) 

if  (trace)  print  «, * >k i  1 1 ' 
nhdfe  =  knceal  (tgt)-l 
narey  =  army  (firer) 

IF  (sot  ion  (tgt)  .ne.  STATNY)  THEN 
ncase  =  3 

ELSE  IF  (motion (firer)  .ne.  STATNY)  THEN 
ncase  =  2 

ELSE 

ncase  =  1 
END  IF 

c  Find  kill  level 

c  Change  12-9-89  by  HLReed  for  interpolation  on  range  for  pki 1 1 
c  Get  ratio  based  on  500  eeter  intervals 
r  =  r/500. 

c  Is  range  >  4000  eeters  if  so  then  use  3999.5 
if (r  .GE.  8.)  r  =  7.999 
c  Get  integer  part 
jrg  =  int(r) 
c  and  fractional  part 

r  »  r  -  float(jrg) 

c  Correct  for  the  fact  that  indices  start  at  1  rather  than  0 
jrg  *  jrg  ♦  1 
teap  s  ranu(0.0) 

IF  (teap  .gt.  p(l))  THEN 
c  no  hit  and  no  kill 

hit  ■  .false, 
injury  *  ALIVE 

S-SEIF  (teep  .gt.  p(2))  THEN 
c  a  hit  and  a  'k*  ki  1 1 

hit  *  .true, 
injury  =  KKILL 

ELSEIF  (teep  gt.  p (3) )  THEN 
c  a  hit  and  an  'aif'  kill  but  no  'k* 
hit  =  .true, 
injury  =  MFKILL 
ELSEIF  (teap  .gt.  p (4) )  THEN 
c  a  hit  and  an  *f*  kill  but  no  'a*  kill 
hit  =>  .true, 
injury  *  FKILL 

ELSEIF  (teap  .gt.  p(5))  THEN 
c  a  hit  and  an  'a'  kill  but  no  *f*  kill 
hit  e  .true, 
injury  >  UK ILL 
ELSE 

c  a  h i t  but  no  kill 

hit  »  .true, 
injury  ■  ALIVE 
ENDIF 

if  (injury . eg. ALIVE)  kshot(naray,10)  s  kshot(naray,10)*l 
if  (injury .eg. MKILL)  kshot(naray, 11)  *  kshot(naray, 11). 1 
if  (injury .eg. FKILL)  kshot(narey, 12)  »  kshot(naray,12)»l 
if  (injury. eg. IFKILL)  kshot(naray,13)  *  kshot(naray,13)*l 
if  (injury .eg. KKILL)  kshot(naray, 14)  *  kshot(naray, 14) « 1 
if  (trace)  print  *, *<ki II* 

END 
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FUNCTION  KILLS  (firer,  tgt,  x,  y) 
c  9  KillS:  find  kill  type  for  a  STAFF-1  ike  round, 
c  sector  -  angle  band 

c  band  -  range  band  (distance  above  tgt). 

include  ’coaeon. h’ 
integer  firer,  tgt 
integer  fan,  fans,  band,  sector 

eoaaon  /cpkh5/  anglia(4),  pkh5(2,7,4,4,12) ,  yl,  y2,  y3,  fans 
save  /cpkhS/ 
diaension  pksave(4) 

1  foraat  ('  KILLS:  narey,  ht,d,yl,y2,y3  =',  i 4 , Sf 8 . 1) 

2  foraat  ('  kill  leveU’.il,1  ran=',f5.3,‘  p(a,f  ,af  ,k)  =  ’ ,4fS.2) 

3  foraat  (’  KILLS:  band,  angles,  sector,  fan=’, i 4 , f 8 . 1 , 2 i 4) 

4  foraat  (’  KILLS:  x,y=’,2f8.3) 

c 

i f  (keyd (4) . gt . 0)  print  • , 1 >k i 1 1 S 1 
fan  =  nang 

if  (nang. gt. fans)  fan=14-nang 

pk  »  ranu(0.0) 

naray  =  aray (firer) 

ht  =  y-y2-y3 

d  =  sqrt(x«»2  ♦  ht*»2) 

IF  (d. Ie.yl-y2)  THEN 
kill  =  ALIVE 

ELSEIF  (d.ge.yl-y2-70.0)  THEN 
kill  x  ALIVE 

a.sE 

band  =  7-int((d-yl-y2)/10.0) 
anglex  =  atan2(abs(x) ,ht)»0Efl 
angles  =  90. -anglex 
sector  *  i ndexx (angles, 4, anglia) 
pksave(l)  *  pkh5(naray, band .sector, 1, fan) 

pksave(2)  =  pksaveflJ-pkhS^naray.band.sector^.fan) 
pksave(3)  =  pksave(2)*pkh5(naray, band, sector, 3, fan) 
pksave(4)  =  pksave(3)-pkh5(naray, band, sector, 4, fan) 
c  Find  «Hich  kill  type  occurs. 

if  (pk. It.pksave(4))  kill  *  KKILL 

if  (pk. It.pksave(3))  kill  =  VFKILL 

if  (pk. It.pksave(2))  kill  *  FKILL 

if  (pk. It.pksave(l))  kill  *  UK ILL 

if  (pk.ge.pksave(4))  kill  at  ALIVE 

if  (ki 1 1 .eq. ALIVE)  kshot  (naray, 10)  *  kshot (naray, 10) -1 
if  (ki 1 1 .eq.UKILL)  kshot(naray, 11)  *  kshot(narey, 11)-1 
if  (kill. eq .FKILL)  kshot (naray, 12)  =  kshot(naray,12)«l 
if  (ki 1 1 ,eq.kFKILL)kshot(naray, 13)  =  kshot(naray,13)-l 
if  (kill. aq . KKILL)  kshot(naray, 14)  =  kshot(narey,14)-l 
ENDIF 
kill 5=k i 1 1 

if  (keyd (4) .gt.0)  print  • , 1 <k i 1 15’ 

END 
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SUBROUTINE  UTE  KL  (t,  tgtjj) 

c  3  Late  kl:  Siaulata  recognition  of  aif  kill  after  period  of  inactivity, 
include  ’coaaon.h' 
integer  firert  tgt 

1  foraat(f8.2,lx,a4,  i 3 , 1  I-killed.’) 

c 

if  (trace)  print  •,’>latekll 

if  (keyd(l) .gt.l)  print  1,  t,  color(aray(tgt)),  tgt 
firer  =  1 

if  (tgt. le.nblu)  firer=nblu»l 
I ife(tgt)  =  IKILL 
call  cancel  (tgt,  ’ikill  ’,NLLL) 
call  nevtgt  (t, firer, tgt) 
call  deaths (t) 
if  (trace)  print  » ,  ’ < I atek I  * 

END 

c  V7.1 
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SUBROUTINE  MAX  VH_(t,  firer) 
c  6  Max  vel:  siaulate  tank  reaching  cruise  speed, 
include  ’coeson.h’ 
integer  firer,  tgt 

1  foreat  (f8.2,lx,a4, i3, *  at  full  speed.') 

c 

if  (trace)  print  «,’)eaxver 

if  (keyd(l).ge.2)  print  1,  t,  color(arsy(f irer)),  firer 
call  path (f irer, t,sotion(f irer), 8.0, x,y,»x,vy) 
sot  ion  (firer)  =  UAXYL 
tgt  *  nrtgt (firer) 

IF  (tgt.gt.0)  THEN 

if  (I ife(tgt) . It.IKILL)  call  engage(t,t,f irer, nrtgt(f irer)) 
ENOIF 

if  (trace)  print  (.'(eaxyel1 
END 

c  V7.4 
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SUBROUTINE  UAYHIT  (t, I. it, n,k,vl,v2, expose, hit) 
c  0  U ayhit:  Find  what  the  round  does. 

c  Changed  May  18,  1989  for  simplified  hit  and  kill  aodel,  KL  Reed 

include  1 coaeon. h‘ 
logical  hit 
integer  expose 

if  (trace)  print  »,  ’>aayhit’ 
kshot(n,6)  =  kshot(n,6)*l 
c  Find  whether  a  hit  occurs, 

hit  =  .false, 
rgx  =  rgf (t,I, it) 
r  *  rgx 

call  ki 1 1  (I, it, hit, injury, r) 

IF  (hit)  THEN 
c  Treat  hit. 

kshot(n,8)  =  kshot(n,8)*l 

if  (I ife(it) .eq.UFKILL)  nhot(it)=nhot(it)*l 

if  (nhot(it) .gt.nbuap(n))  call  skedu I  (t, it, '  ik i 1 1  '.NULL) 

know (I, it) =2 

prerrd(I)  =  2 

IF  (reliab(n)  .ge.  ranu(0))  THEN 
call  daaage(t,  I,  it,  injury) 

B.SE 

c  Round  is  a  dud. 

kshot(n,9)  =  kshot(n,9)»I 
ENDIF 
ELSE 

c  Treat  aiss. 

kshot(n,7)  *  kshot(n,7)-l 

IF  (psense(n,nrgf (rgx, rgincr}}  .gt. ranu(B.B))  THEN 
prevrd(I)  «  4 

c  Changed  by  HLReed.  seeas  to  be  needed 

know(I, i t)=l 

if  (keyd(l) .ge.2)  print  Uiss  is  sensed.’ 

e_sE 

prewrd(I)  *  3 

if  (keyd(l) .ge.2)  print  Uiss  is  not  sensed.' 

ENDIF 

ENDIF 

c  Careful.  If  either  aoving,  aake  sure  nx  rd  is  treated  as  1st 

c  round  if  SS  case  occurs. 

if  (vabs(vt) .gt.0  .or.  vabs(vf) . gt . 0)  pre»rd(I)=l 
if  (vl.gt.0  or.  v2.gt.0)  prevrd(I)=l 
if  (trace)  print  *,  '<eayhit' 

END 

c  Needs  data  checking, 
c  V7.1 
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SUBROUTINE  UK  TBL  (onlye.onlyf ,fand«,«disp,k,nro») 
e  3  UV  tb I :  aaka  head-on  pkh  table  for  echo, 
include  'coeaon.h1 

coeaon  / cpkh/  table  (4,12),  echo(2,7,7),  jrg(2,7),  jdisp(2,7) 
c 

if  (trace)  print  *,’>iktbll 
jrg(k.nroe)  *  (nro«-l)*irginc 
00  11  j =1 , 4 

echo(k,nro»,  j)  =  table(j.l) 

11  CONTINUE 

echo (k , n roe, S)  =  onlya 

echo(k,nrow,6)  =  onlyf 

echo(k,nro«,7)  =  fande 

jdisp(k,nro»)  =  edisp 
90  if  (trace)  print  «,’<ektbl‘ 

END 

c  V7.2 
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SUBROUTINE  NETIGT  (t,  firer,  tgt) 
e  3  New  tgt:  redirect  all  ’attackers'  of  tgt  to  a  new  target, 

c  New  tgt  called  for  non-fa  I se  tgts  only  and  only  if  tgt  condition 

c  warrants  it.  It  should  only  be  called  if  tgt  is  V-killed, 

c  vanishes,  or  hides. 

c  Maybe  it  should  be  called  if  the  tgt  is  I-kilted  by  a  gun  systea. 

include  ' coaeon. h1 

integer  first,  firer,  tgt,  arayf,  arayt 
logical  loaded,  hav  aao,  cango 

1  foraat(f8.2,lx,a4, i 3 , 1  dis-engs  1 ,a4, 13,20*, *ftgts=’ , i 2) 

2  foraat(f8.2,  lx,  a4,  i 3 ,  1  begins  to  reload.’) 
c 

if  (trace)  print  •,‘>newtgt‘ 
c  Find  first  and  last  ’attacker’ 

first  =  1 

if  (firer. gt.nblu)  first  =  nblu»l 
last  »  nblu 

if  (firer. gt.nblu)  last  =  nblu-nred 
arayf  *  a ray (first) 
arayt  *  3- arayf 
kind  a  kindrd(arayf) 
nrpb2  a  nrpb(arayf) 

DO  20  j=first,  last 

IF  ((aot(j,tgt)  .or.  fot(j.tgt))  .and.  I ife(j) . It.FKILL)  THEN 
IF  (kind. I e. 2  .or.  kind. eg. S)  THEN 
c  Single  snot  gun  systea  or  other  fire  i  forget  systea. 

IF  (nrpb(arayf) . le.l)  THEN 
c  Single  shot  gun  systea. 

call  cancel (j, 'fire  ’,tgt) 
if  (nrtgt(j) .eq.tgt)  busy(j)  =  .false, 
if  (nrtgt(j). eq.tgt)  nrtgt(j)  =  0 
hav  aao  a  nrd(j) .  It.nrds(arayf) 

IF  (hav  aao)  THEN 
thuaan  a  0.«exp(rol In (0.5)) 
call  se I ecs(t,j, thuaan) 

B.SEIF  (can  go (j ,t) .and. ishtfs (arayf ) .gt.0)  THEN 
c  Move  out 

call  skedul (t, j , 'accel  ’.NULL) 

END  IF 

nrot(j)  a  0 
fot(j.tgt)  =  .false. 

if  (keyd(l) .ge.2)  print  l,  t,  col  or (arayf) ,  j, 

1  color(arayt) ,  tgt,  nchan(j) 

ELSE 

c  Burst  fire  gun. 

print  «,'NEWTGT:  Not  iapleaented  for  burst  fire.’ 

STOP 

ENDIF 

ELSEIF  (kind. eq . 4)  THEN 
c  Guided  aissile  branch. 

if  (fot(j.tgt))  call  cancel (j , 'fi re  ’,tgt) 
if  (not (j, tgt))  call  abort (t , j , tgt) 
loaded  a  nchan(j)  .eq.  nchans(arayf) 

IF  ((.not.eapty(j)  .and.  aot(j,tgt)  .and.  loaded)  .or. 

1  (.not.eapty(j)  .and.  fot(j.tgt)))  THEN 

IF  ((aod(nrd(j),nipods(arayf))  .gt.  0)  .or. 

1  fot(j.tgt))  THEN 

c  More  rds  in  pod 

IF  (fot(j.tgt))  THEN 
call  cancel (j, 'select' .NULL) 
busy(j)  «  .false. 

Tot  (j, tgt)  a  .false. 

ENDIF 

C  if  (tgt.ne.FLS  TGT)  fot(j,tgt)  =  ,f:.ise. 

thuaan  a  0.0«exp(rol I n (0 . S) ) 
call  selecs(t, j , thuaan) 
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bjse 

c  Treat  eapty  aissile  pod 

eapty(j)  «  .true, 
call  cancel (j, ’fire  ’,tgt) 
call  cancel  (j/select', NULL) 
busy(j)  *  .false. 
nrot(j)  ■  0 

c  shod  htf  that  is  slowing  to  engage  speed  up  now? 

call  skedul  (t*trelod(ariyf) ,j, 'reload' ,NILL) 
if  (keyd(l).ge.2)  print  2,t,color(areyf),j 
ENDIF 
ENOIF 

nchan(j)  =  nchan(j)  -1 
nrtgt(j)  =  8 
fot(j,tgt)  *  .false. 

ENDIF 

ENDIF 

if  (seen(j , tgt) )  ndet(j)  =  ndet(j)  -  1 
seen(,j,tgt)  =  .false. 

20  CONTINUE 

IF  (.not. repeat)  THEN 
repeat  =  .true. 

call  skedul  (t*. 01,0, 'search' .NULL) 

ENDIF 

c  Following  reaoved  cause  duplicates  code  in  deaths  routine 

c  and  is  incorrect, 

c  ndead=0 

c  IF  (arayt.eq.RED)  THEN 

c  DO  30  i=nblu«l,nblu«nred 

c  if(life(i).gt.FKILL)  ndead=ndead*l 

c30  CONTINUE 

c  if  (ndead.eq.nred)  call  skedul (t»8. .NULL, 'finish', NULL) 

c  BJSE 

c  DO  40  i=l,nblu 

c  if  (I ifa(i) .gt.FKILL)  ndead*ndead*l 

c40  CONTINUE 

c  if  (ndead.eq.nblu)  call  skedul (t«5. .NULL, ’finish* .NULL) 

c  ENDIF 

if  (trace)  print  e/Cnewtgt' 

END 
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SUBROUTINE  NEXTEV  (I , act, it, t) 
e  0  Nextev:  Find  the  next  scheduled  event, 
include  ‘clock. h’ 
character«6  act 
c 

c  Fill  arguments 

I  *  eho(nxevnt) 
act  *  what (nxevnt) 
it  3  ehoe(nxevnt) 
t  -  ehen(nxevnt) 

c  Drop  storage  unit  froe  active  storage  chain 

n  =  nxevnt 

nxevnt  *  next (nxevnt) 

c  Add  storage  unit  to  inactive  storage. 

next(n)  »  nxidle 
nxidle  *  n 
END 

c  V7. 1 
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c  9 

CHANCED 


c  V7.2 


FUNCTION  NRCF  (rg.rgincr) 

Nrgf :  find  which  rgincr  eeter  rg  band  range  is  in. 

1  Apr  86.  Next  line  changed. 

nrgf  *  eax0(l,int(  (250 . *rg) /500 .  )) 

nrgf  =  eax0(l, int(0.S*rg/rgincr)) 

END 
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SUBROUTINE  NXWAVE 

c  9  Nx  wave:  Simulate  all  reps  for  the  nth  engageaent. 
include  ’coeeon.h1 
characters  str(3) 
d ieens ion  i statb (8) 
data  str  /'Utg  ’ , ’Ratk’ , ’Batk1/ 

c 

if  (trace)  print  «,,>Nxwavel 
nrg  *  rg0/irginc 
nreps2  *  nsurv/nblu 
neval  =  iod(n5urv,nblu) 
noaiao  -  0 
loaaeo  *  0 
nsurv  =  neval 
DO  20  i =1 , 8 
statb(i)  *  0.0 

20  CONTINUE 

c  Changed  by  HLReed  for  printing 
DO  21  i=l,4 
eystat(i)  =  0 

21  CONTINUE 

DO  30  nrep  =  l,nreps2 
krep=nrep 
call  events 

30  CONTINUE 

c  Update  statistics  after  all  reps  of  nth  engageaent. 

c  Changed  by  HLReed  for  printing 

teapl  =  f loat(aystat(l))/f loat(nreps2) 
teap2  =  f loat(aystat(2))/f loat(nreps2) 
teap3  »  f loat(aystat(3))/f loat(nreps2) 
teap4  >  f loat(aystat(4))/f loat(nreps2) 

DO  40  1=1,0 

statc(i)  a  statc(i)  «  statb(i) 

40  CONTINUE 

noaeo2  a  noaeo2*noaaao 

loaao2  a  loaao2«loaaao 

DO  50  i«l,4 

istatb(i)  =  0.5  ♦  100*statb(i)  /  nreps2 
statb (i -4)  =  statb (i *4)  /  nreps2 

50  CONTINUE 

c  Changed  by  HLReed  for  printing 

2  ^18^5.0^8.3^6.3/6.3/6.3, i4, i 4 , f 6 . 2 , f 8 . 3, f 8 . 3) 

exchc  a  0 

if  (statc(S)  gt.0)  exchc  a  statc(6)/state(5) 
print  2 , rg0 , teapl , teap2 , teap3 , teap4 , i statb (1) ,  i statb (2) , 
1  exchc, statb (7) ,  statb(8) 
nreps3  a  nreps3»nreps2 
if  (trace)  print  a/Cnxwave’ 

END 
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SUBROUTINE  PATH  (firer, t,  aotio2,  dalt,  x,  y,  vx,  vy) 
c  *  Path:  search  path  table  for  position  and  vel  at  tiee  t. 
include  'coeeon.h' 
logical  is  atkr,  kan  go,  old 
integer  firer 
c 

if  (trace)  print  «, ’>path’ 
narey  =  arey  (firer) 

is  atkr  *  (score. eq.RATTAK  .and.  narey. eq. RED)  .or. 

1  (scene. eq.BATTAK  .and.  narey. eq.BLU) 
kan  go  •*  (aotio2.ne.STATNY  .or. 

1  I ife(f i rer) .eq.ALlVE  or.  I ife(f irer) .eq.FKILL) 
dt  *  t-t0(firer) 

old  =  dt  .gt.  dalt  _ 

IF  (is  atkr  .and.  kan  go  .and.  old)  THEN 
c  Update  positions  and  velocity. 

t0(f  i  rer)  =  t 
if  (eot i o2 . eq . SLOW NC)  THEN 
dv  =  decel (nareyi*dt 

y0  (firer)  =  y0(f  i  rer)-dt*(vy0(f  i  rer)-0 . 5*dv) 
v  =  vy0(f  i  rer)-dv 
if  (abs(v). It. 0.001)  v  =  0.0 
vy0(firer)  =  v 

ELSEIF  (eot i o2 . eq . STATNY)  THEN 
vy0(firer)  =  0.0 
E.SEIF  (eot  i  o2 .  eq .  ACCB.G)  THEN 
dv  =  accel (narey) *dt 

y0(f  i  rer)  =  y0(f  i  rer}*dt«(vy0(f  i  rer)-0.  S»dv) 
vy0  (firer)  =  vy0(f irer)*dv 
B.SEIF (eot i o2 . eq . UAXVL)  THEN 
y0 (firer)  *  y0(f  i  rer)-vy0(f  i  rer)vdt 
vy0(f i rer)  *  speed(narey) 

asE 

print  v/PATH:  no  such  aotion.  aotio2=,  *  ,eoti^2 
STOP 
ENDIF 

IF  (acceax(naray)  .ne.0.)  THEN 
c  Add  sinusoidal  aotion 

oaegs  a  T10PI/evlth(naray) 

x0(firer)  =  aepl (narey)»sin(oaega»y0(f irer))* 

1  xp(firer) 

vx0(f  irer)*  aepl  (narey)»cos(oeega> y0(f  i rer))* 

1  oaega*vy0 (firer) 

ENDIF 

ENDIF 

x=x0(f  irer) 
y=y0(firer) 
vy=vy0(f i rer) 
vx=0.0 

if  (trace)  print  «,,<path' 

ENO 
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SUBROUTINE  PINPNT  (t.firer) 

e  8  Pinpnt:  Simulate  firing  signature  (pinpoint)  detection  by  soae  foes, 
include  'coeeon.h' 
integer  first,  firer 
logical  eilsee 

1  foreat  (f8.2,lx,a4, i3, 1  sees  ',84, i 3 , 1  euzzle  flash’) 

c 

if  (trace)  print  «,’>pinpnt’ 
first  =  1 

if  (f i rer. le.nblu)  first  =  nblu»l 
last  3  nblu 

if  (firer. le.nblu)  last  =  nblu*nred 
pinpxx  =  pinp(arey(f  irst)) 

DO  20  i=first,  last 
wilsee  =  pinpxx. gt. ranu(0.0) 

IF  (I ife(i) . It.FKILL  .and.  eilsee  .and. 

2  ndet(i) . It.ndets(arey (i))  .and. 

1  I o* ( i .firer)  .and.  .not. seen (i ,f i rer) )  THEN 

if  (keyd(l) .ge.2)  print  1, 

1  t,  color(arey(i)),  i,  color(arey (f i rer)) ,  firer 

seen (i ,fi rer)  *  .true. 
ndet(i)  *  ndet(i)  ♦  1 

c  Change  by  HLReed  added  nee  variable  for  pinpoint  tiae  and 

c  transfer  of  pinpoint  detection 

thuean  =  pntiee(arey(i))  »exp(rol I n (0 . 5) ) 
call  selecs(t, i .thuean) 

if (xxfer(arey(i)))  call  skedul (t»thuean,f i rst, 'xfer  ’ , f i rer) 

ENDIF 

20  CONTINUE 

if  (trace)  print  »,’<pinpnt’ 

END 
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SUBROUTINE  ^OP  DN  (t, finer) 
c  0  Pop  dn:  Have  defender  pop  down  to  reload? 
include  ’coeaon.h’ 
integer  firer 
c 

if  (trace)  print  «,'>pop  dn’ 
call  vanter(t, firer, NULL) 
if  (trace)  print  «,’<pop  dn1 
END 
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SUBROUTINE  PR  GAME 

c  9  pr  gaae:  print  gaae  control  constants, 
include  ’coaaon.h' 

1  f  oraat (21* , 1 f 81 ues  f  Reds ’ ,  / , 

1  '  Meeting  eng:;»«ent:  ',3i2,4x,3i2,/, 

2  *  Red  attack:  * , 3 i 2, 4* , 3 i 2. / , 

3  ’  Blue  attack:  ',3i2,4x,3i2) 

2  foraat  ('  DO  rg=’,3i6,’  (opening  ranges)') 

3  foraat  (20 i 2) 

4  foraat  (’  nreps  * ’ , i 5 , ’  nwaves  *',13,'  iangd  =’,i2) 

5  foraat  (/,20(’='),’RUN  DESCRIPTION', 20 ('=*)) 

6  foraat  (55('*'),/) 

c 

if  (trace)  print  *,’)prgaae’ 
print  5 

print  l, ((ntanks(i , j) , j=l,C) , i=l,3) 
print  2,ainrg,aaxrg, incrg 
print  3,  keyd 
print  3,  key* 

print  4,  nreps,  nwaves,  iangd 
IP  (aeth  sa.eq.l)  THEN 

print  *,'AMSAA  US  error  treataent. ' 
a.SEIF  (aeth  sa.eq.2)  THEN 

print  »,'8RL  US  error  treataent.' 

ELSE 

print  », 'Error  treataent  value  'aeth  sa*  wrong.’ 
ENDIF 
print  6 

if  (trace)  print  «,’<prgaae‘ 

END 
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INTEGER  FUNCTION  PRIORN  (t,  firer,  lev  old) 
c  6  Priorn:  select  tgt  with  highest  priority, 
include  'coaaon.h1 
logical  better,  ck  tgt 
logical  pick 
integer  firer,  areyf 
c 

if  (trace)  print  »,’>priorn' 
areyf  =  arey  (firer) 

c  'eake'  dueiy  tgt  for  comparison 

rg  old=l.e3S 
t  old=l.e3S 
lev  old=1000 
priorn  *  NULL 
last  *  nblu«nred 
DO  30  itgt=l,last 

c  Coipare  all  possible  targets 

pick*. true. 

c  Don't  select  this  target  if  I'a  already  servicing  it. 

i f (sot (f i rer , itgt) .or  f ot (f i rer , itgt) ) p i ck= . f  a  I se . 
rg  tgt  *  rgf  (t,f i rer, itgt) 

ck  tgt  =  seen (firer, itgt)  .and.  I ife(itgt) . It.IKILL 
I  .and.  .-gtgt. I e. 4000.0  .and.  pick 
IF  (ck  tgt)  THEN 

c  Firer  sees  tgt,  it's  threatening,  t  he's  not  firing  at  it. 

call  priort(f i rer,  itgt,  rg  tgt,  t,  level) 
c  Change  by  HLReed  for  say  share  is  handled 

if  (share (anyf))  then 
III  «  0 

do  20  jjj  *  1, last 

if (lot (jjj, itgt) -or.fot(jjj,itgt))  1 1 1*30 
20  continue 

level  *  level  *111 
end  if 

c  Noe  pick  the  tgt  with  highest  priority 

rg  tgt  *  rg  tgt  «(l».0S«rol ln(1.0)) 
t  tgt  *  tf i re (f i rer, itgt) 
better  *  level  .It.  lev  old 
IF  (lev  old.eq. level)  THEN 
c  Sale  priority  class;  now  break  ties 

c  if  new  tgts  pick  closer 

if  (t  tgt.le.  0)  better  *  rg  tgt  .It.  rg  old 
c  if  old  tgts,  pick  older  (least  recently  fired  on) 

if  (t  tgt.gt.  0)  better  *  t  tgt  .It.  t  old 
ENDIF 

-  IF  (better)  THEN 
lev  old  *  level 
t  old  *  t  tgt 
rg  old  *  rg  tgt 
priorn  *  itgt 
ENDIF 
ENDIF 

30  CONTINUE 

if  (trace)  print  •, '(priorn1 

END 
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SUBROUTINE  PRIORT (firer,  tgt,  rg  tgt,  t,  L) 
e  0  PRIORT:  find  priority  of  tgt  (w/  preference  to  old  tgts) 
include  ' coaeon. h' 
integer  firer,  tgt 
diaension  I ev (21 , 2) 
save  lev 

data  lev/1, 2,3, 4, 6,6,7,8,9,10, 11, 12, 13, 14, IS, 16, 17, 18, 19, 20, 21, 
.  1,2,3,4,5,8,7,8,9,16,16,17,18,19,20,21,10,11,12,13,14/ 

1  foraat(’  PRIORT:  ’ ,a4, i 3, “  considrs  ’,a4,i3,’  with  sriority', 

1  <4.  *  C,  i2,  *)*) 

c 

if  (trace)  print  •,’>priort’ 
j  =  nprior(aray(f  irer)) 
a  s  aotion (tgt) 
t  activ  s  l.e35 

if  (tf i re2 (tgt) . gt . 0 . )  t  activ  =  t-tf ire2(tgt) 

IF  (tf i re (f irer, tgt) . gt . 0)  THEN 
c  firer  has  already  shot  at  this  target  previously 

IF  (know (f irer, tgt) .eq.l)  THEN 
c  Missed  target  with  last  round  fired  at  it 

IF  (rg  tgt. It.recknz(aray(f irer)))  THEN 
c  Target  is  within  recognition  range 

L  *  5 

if  (a.eq.STATNY  .or.  a.eq.SLOING)  L  =  3 
if  (nrtgt (tgt) .ne.0)  L  *  2 
if  (t  activ  .It.  30.)  L  *  1 

a.sE 

c  Target  is  beyond  recognition  range 

L  «  7 

if  (a.eq.STATNY)  L  *  6 
if  (t  activ  .It.  30.)  L  «  4 
ENDIF 
B.SE 

c  Hit  target  with  last  round  fired  at  it 

IF  (rg  tgt.  It.reckni(aray(f  irer)))  THEN 
c  Target  is  within  recognition  range 

l  -  14 

if  (a.eq.STATNY  .or.  a.eq.SLOING)  L  *  12 
if  (nrtgt (tgt) .ne.0)  L  *  11 
if  (t  activ  .It.  30.)  L  *  10 
B-SE 

c  Target  is  beyond  recognition  range 

L  =  16 

if  (a.eq.STATNY)  L  =  15 
if  (t  activ  .It.  30.)  L  «  13 
ENDIF 
ENDIF 
ELSE 

c  Target  is  a  new  target 

IF  (rg  tgt. It. reeknz(ariy(f i rer) ) )  THEN 
c  Target  is  within  recognition  range 

l  *  19 

if  (a.eq.STATNY  .or.  a.eq.SLOING)  L  *  18 
if  (nrtgt (tgt) .ne.0)  L  *  17 
if  (t  activ  . It.  30.)  L  «  8 
B_SE 

c  Target  is  beyond  recognition  range 

L  »  21 

if  (a.eq.STATNY)  L  ■  20 
if  (t  activ  . It.  30.)  L  *  9 
ENDIF 
ENDIF 

L  *  la»(L,  j) 

if  (trace)  print  •,'<priort' 

END 
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SUBROUTINE  PR  MISC (flirty) 
c  3  Pr  aisc:  print  aisc  tank  characteristics, 
include  'coaaon.h' 
integer  irg(8) 

10  foraat(/, '  ==^=«g===BLUE  SYSTBJ  DESCRIPTION’, 

1 

U  foraat (/,’  *xx=xx=xxx=x=x=x=xxx=RED  SYSTEM  DESCRIPTION’ , 

1  ’XXXXXXXXXXXXXXXXXXX,’) 

12  foraat(/, 

1  '  SYSTEM  DIMENSIONS', 21a, 'MOTION  CHARACTERISTICS’,/, 

1  ’  Distance  (a)  froa  center  of’, 

1  llx, ’Acceleration’, f6. 2, '  a/s»»2',/, 

1  '  turret  ring  to: ',23x, 'Deceleration' , f 6 . 2, ’  a/s«»2',/, 

1  ’  Turret  top  ’,f6.2,’  Ground  ’,f6.2, 

1  ’  Tine  to  hide’, f6. 2,’  sec’,/, 

1  ’  Turret  Side  *,f8.2, *  Hull  side  ’ ,fS.2 ,/, 

1  ’  Turret  front’ ,f6. 2, ’  Hull  front’, f8. 2,/, 

1  ’  Turret  back  ’  ,f6.2,  ’  Hull  back  \f6.2 ,/) 

13  foraat(/, 

1  9x, ' - DETECTION  CAPABILITY - ’, 

2  '  - FIRING  CYCLE - ’ 

I  ,/,'  Rg  Psense  P-detect  T-aedian  ', 

1  '  Tfirst  Tfixed  Tf ly ’ ,/, 

1  '  (a)  HD  FE  FE-M  HD  FE  FE-M’, 

1  '  (sec)  (sec)  (sec)’,/, 

1  8  ( i  7 , 2f  7 .2, 2f  8 . 2 ,  f  3 . 2 , 2f  7 . 2,  f  7 . 1 ,  f  8 . 1 ,  f  8  -  2 ,  /) ) 

15  foraat('  Systea  aaao  load  is’,i3,'  KE  rounds') 

16  foraat(’  Systea  aaao  load  is',i3,’  HEAT  rounds’) 

17  foraat(’  Systea  aaao  load  i * ' , i 3 , ’  wait-ti l-iapact  aissiles’) 

18  foraat(’  Systea  aaao  load  is',i3, 

1  '  siaultaneous-f I ight  aissiles  with' , i 2, ’/pod  I’,/, 

2  '  reload  tiee«‘,f7.1) 

19  foraat(’  Systea  aaao  load  is’ , i3, ’  STAFF- 1  ike  rounds') 

20  foraat(*  Systeas  can  engage', i 2, 1  targets  siaul taneously . ’) 

21  foraat( 

1  ’  Switch  targets  after:',/,’  1.  A  IC-k ill’, 

2  1  (i  don”t  re-engage)’,/, 

1  ’  2.  An  MtF-kill  and’ , i2, *  Hits  or’.fS.l,’  sec.’, 

2  ’  (I  don”t  re-engage)') 

22  foraat(’  3.  After  scoring  a  hit.') 

23  foraat(*  3.  After', i3 , '  shots.') 

24  foraat('  Return  to  partially  serviced  target  after  vainly',/, 

1  '  searching' ,f7. 2, '  sec  for  a  new  target.’) 

25  foraat('  Systea  halts  to  fire.') 

26  foraat('  Systea  fires  on  the  aove.') 

27  foraat('  Median  tiaa  between  rounds  is’,f6.2,'  sec.') 

28  foraat('  Systea  fires’, i2,'  round  bursts  with  rounds’, 

1  ’  spaced', f6. 2,'  sec  apart.’) 

29  foraat('  There  are',  i3, '  decoys, ’ , i 3 , 

1  ’  of  which  are  flashing.') 

30  foraat('  Miniaua  tiae  to  fire  next  round  i » ’ , f 6 . 2 , ’  sec.’) 

31  foraat('  Probability  of  firing  signature  detection  is’,f6.2) 

32  foraat('  Reliability  is ' , f8 . 2) 

33  foraat('  Recognizes  targets  inside' , f 8 . f, '  aeters.') 

34  foraat('  Probability  of  picking  falsa  HO,  FE  tgts  are',2f6.3) 

35  fsraat(’  Selects  old,  hit  targets  before  new  targets.') 

36  foraat('  Selects  new  targets  before  old,  hit  targets.’) 

37  foraat(‘  Tgt  sharing  is  set.  Ion"t  pick  a  tgt  being  serviced  ') 

38  foraat(’  Tgt  sharing  is  off.  fill  pick  a  tgt  being  serviced.') 

39  foraat('  Systeas  can  detect’, i2,'  targets  siaul taneously.  ’) 
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50 

c 


c 


c 


c 


c 


if  (trace)  print  », 'Jpraise’ 

00  50  i*l,8 

irg(i)  «  ieirginc 
CONTINUE 
Write  header 

if  (naray. eq.l)  print  10 
if  (naray. eq.2)  print  11 

Write  systea  dieensions  and  notion  characteristics 
print  12,  accel (naray) ,decel (narey) , 

1  sysd ie (narey ,1),  sysd i* (narey, 5),  thide(narey), 

1  (sysd i a (na  ray ,  i ) , sysd i a  (na  ray , i *4) , i=2 , 4) 

Writs  range  dependent  values 
print  13, 

1  ( < rg ( i ) .  psense(naray.i),  (pinf in(naray, j, i),j=l,3), 

2  (tbar(naray, j , i) , j=l,3) , 

1  tf irst(naray, i),tf ixed (naray, i), 

1  tof  (naray, i),i*l,8) 

Target  saitching  policy 
i  *  tactic (naray) 

print  21,  nbuap(naray) ,  tbuap(naray) 
if  (i .eq.2)  print  22 

if  (i.eq.3)  print  23,  nrpt(naray) 

if  (i.gt.l)  print  24,  tlook(naray) 

print  39,  ndets(naray) 
print  20,  nchans(naray) 
if  (npri or (naray) .eq.l)  print  35 
if  (nprior(naray) .eq.2)  print  36 
print  • 

Write  projectile  inforaation 

if (kind  rd(naray) ,eq.l)print  IS,  nrds(naray) 
if (kind  rd(naray) ,eq.2)print  16,  nrds(narey) 
if  (kind  rd(naray) ,eq.3)print  17,  nrds(naray) 
if (kind  rd(narey) ,eq.4)print  18,  nrds(naray), 

1  n ipods (naray),  trelod(naray) 

if (k ind  rd(naray) ,eq.5)print  19,  nrds(naray) 

if ( ishtfs (naray) .gt.0) print  25 

if (ishtfs(naray) .eq.0) print  26 

print  27,  taed in(naray) 

print  30,  tain(naray) 

print  32,  rel iab(naray) 

print  31,  pinp(naray) 

print  33,  recknz(naray) 

print  34,  (pfa Isa (naray, i), i =1 , 2) 

if (nrpb (naray) .gt. 1) print  28,  nrpb(naray),  rof(naray) 

if (ndecoy (naray) .gt.0)  print  29,  ndeeoy (naray) , 

1  nf I  ash (naray) 

if  (share(naray))  print  37 
if  (.not.share(naray))  print  38 
if  (trace)  print  «, '(praise1 
END 
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FUNCTION  RANU  (da) 

c  Rami:  A  version  of  uran31  randoa  unifora  nr  generator, 

coaaon  /cranda/  i,  j 
real  al 

j»« 

j=j*2S 

j=j-(j /67108864) <67108864 
j=j<2S 

j = j -  ( j /67 108864) *67108864 
j=j<5 

j = j -  ( j /67108884) <87108864 

*l*j 

1=j 

ranu=  al/67108884 
END 
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SUBROUTINE  RD  MISC  (dbnaae, naray) 
c  9  Rd  aisc:  raid  aiscel laneous  tank  characteristics, 

include  ’coaaon.h’ 
character  dbnaae*32,  line«72 
real  high (2) 

2  foraat(lai) 

3  foraat(a) 
c 

if  (trace)  print  *,’>Rdaiscl 
open (4,  file-dbnaee,  status='old') 
rewind  4 

read (4, •)  (sysdia(naray, i) , i=l,8) 
read(4,«)  (psense (naray,  i) ,  i  =1 , 8) 
c  Read  nvl  outputs. 

read(4,«)  (pinf in(naray,l, j), j=l,8) 
read(4,«)  (pinf in{naray,2,.i) , j=l,8) 
read (4,*)  (pinf in(naray,3, j) , j=l,8) 
read (4 , «)  (tba r (na ray , 1 , j 1 , j =1 , 8) 
read(4,«)  (tbar(naray,2Jj),j=lI8) 
read (4 , »)  (tba  r (na  ray , 3 , ; ) , j = 1 , 8) 
c  Change  by  KLReed  added  pinpoint  tiae  input 

read(4,»)  recknz (naray) , (pfalse(naray, i), i=l,2) , 

1  tlook(naray) ,p inp(narey) , rel iab(narey) , trelod(naray) , 

2  pntiae(naray) 

read (4,«)nrds(naray) 1nrpt(naray),nrpb(naray) , 

1  tact ic (naray) ,k ind  rd(naray) , npr ior(naray) , 

2  n i pods (na ray) ,nchans (naray) 
read (4,*)  (tof (naray, i) , i=l,8) 
read (4 , •)  (tf i rst (naray , i ) , i =1 , 8) 

read(4,»)  taed in (naray),  tain(naray),  rof (naray) 

read(4,«)  (tf ixed (naray, i), i*l,8) 

read(4,«)  accel (naray) ,  decel (naray) ,  speed(naray) , 

1  angle (naray) ,  thide(naray) 

read (4,*)  i shtfs (naray) ,  nbuep(naray) ,  ibuap 
tbuap(naray)  •  ibuap 
read(4,«)  ndecoy (naray) ,  nf I  ash (naray) 
c  Change  by  HLReed  added  xxfer  for  control  of  target  transfer 
read(4,e)  share(naray) ,  xxfer(naray) 
read (4 ,2)  kview(naray) 
read (4 , •)  ndets(naray) 
close  (4) 

if  (keyd(2) .gt.0)  call  pr  aisc  (naray) 
c  Convert  tbar  to  detection  probability  /  second. 

00  30  i*l, 8 
DO  20  j*l,3 

tbar (naray , j , i)  ■  1.0-exp(-l .0/tbar(naray , j , i)) 

20  CONTINUE 

30  CONTINUE 

if  (trace)  print  •/Crdaisc' 

END 
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SUBROUTINE  RDPKH  (dbnaee,  naray) 
c  3  Rd  pkh:  read  probabi I ity-of-ki 1 1  data, 
e  Changed  for  siaplified  hit  and  kill  aodel  klay  19,1989,  HL 

include  'coaaon.h' 
character*32  dbnaae 
coaaon  /cpkh2/  pki 1 1  (2,3,2, 5,9) 
save  /cpkh2/ 

if  (trace)  »rite(»,«) ’>rdpkh* 
open  (4,  file=dbnaae,  status=’old’) 
rewind  4 

00  100  ncase  =1,3 
DO  70  nhdfe=l,2 
00  30  i=l,S 

read  (4,*)  nl,n2,n3, (pk i 1 1 (naray , ncase, nhdfe, i , j) , j=l,9) 

30  CONTINUE 

70  CONTINUE 

100  CONTINUE 

close (4) 

90  if  (trace)  >rite(«,«) '<rdpkh’ 

END 

e  V7.1 
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SUBROUTINE  RELOAD  (t,  firer) 
c  6  Reload:  siauiates  coapletion  of  reloading 
c  30  Oct  8S  Fixed  statement  printing  error  aesaage 

include  ’coaaon.h* 
integer  firer 
logical  defndr 

1  foraat(f8.2,lx,a4, i3, 1  finiahea  reloading’) 

2  foraat(f8.2,lx,a4, i3, 1  popa-up’) 

c 

if  (trace)  print  */>reload' 
naray  *  aray  (firer) 

if  (keyd(l) .ge.2)  print  l,t,eolor(naray) , firer 
nrtgt(firer)  =  0 
eapty (firer)  *  .false. 

defndr  *  ( scene. eq. BA TTAK  .and.  naray. eq. RED)  .or. 
1  (scene. eq. RATTAN  .and.  naray. eq.BLU) 

IF  (defndr)  THEN 

c  Defender  pops  back  up  end  will  start  searching. 

if  (keyd(l) .ge.2)  print  2,t,color(naray),f irer 
call  aprter(t,f irer,tgt,HD) 

ELSE 

c  Attacker  or  ‘aeeter’  never  popped  down, 

thuaan  «  0.0*exp(roi ln(0.5)) 
call  selecs(t  f i rer, thuaan) 

ENDIF 

if  (trace)  print  (/(reload' 

END 
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SUBROUTINE  RESET  (prflg) 
c  0  Reset:  Initialize  the  clock  to  ti 

include  'clock. h’ 
logical  prflg 
e 

prflag  =  prflg 
nxevnt  =>  0 
niidle  =  1 
00  10  j=l , NE 
next(j)  =  j*l 
10  CONTINUE 

next(NE)  =  0 
END 
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FUNCTION  RGF  (t,  firer,  tgt) 

c  3  Rgf:  find  the  position  of  the  firer  w.r.t.  the  tgt. 
include  'coeeon.h' 
integer  firer,  tgt 
coeeon  /pathc  /  xf,  yf,  xt,  yt 
save  /pathc  / 

1  foriat  (9x, 'Firer  x,  y,  vx,  vy  =  ’,  4fl0.1,/ 

*  9x,’Target  x,  y,  vx,  vy  =  ',  4f 10 . 1} 
c 

if  (trace)  print  »,'>rgfl 

call  path  (f i rer,t,eotion(f irer) ,0.0,xf ,yf , vf (1) , vf (2)> 

call  path  (tgt,t,iotion(tgt),0.0,xt,yt,vt(l),vt(2)) 

s(l)  =  xf-xt 

s(2)  =  yf-yt 

s(3)  =0.0 

vt(3)  =  0.0 

vf (3)  =  0.0 

teep  =  vabs(s) 

if (teep.GT.4000)  teep  =  4000 

nrg  =  nrgf (teep.rgincr) 

rgf  *  teep 

rg  =  irgincvnrg 

if  (keys (20) .gt  0)  print  1, 

•  xf,  yf.  »f(l).  »f(2).  *t,  yt,  vt(l),  vt(2) 
if  (trace)  print  «, *<rgf ’ 

END 
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FUNCTION  RNDANG(iangd) 

c  RNOANG:  Draw  a  randoa  angle  froe  a  eardioid/other  distribution, 

c 

PI=3. 1415928538 
denoe  *  0.5/PI 
psranu(duaay) 

c  Do  binary  search  to  find  theta  associated  with  randoe  draw 

tlo  *-PI 

if  (iangd.gt.l)  tlo  =  -PI/3, 
thi  *  PI 

if  (iangd.gt.l)  thi  =  PI/3. 

DO  20  i si , 10 
theta  *  0.5»(t lo-thi) 

if  (iangd.eq.l)  px  =  (theta*sin(theta)*PI)«denoe 
if  (iangd.gt.l)  px  =  (3.«theta»sin(3.*theta)»PI)«denoa 
IF  (px.lt.p)  THEN 
tlo  =  theta 
B.SE 

thi  =  theta 
ENDIF 

20  CONTINUE 

rndang=theta 
theta=theta» 180 . /PI 
END 
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FUNCTION  ROLLN (sigaa) 

c  6  Rolln:  find  a  randoa  nuaber  froa  a  noraal  distribution, 
c  Box -Uu liar  aethod 

save  j,  z 

data  j/0/  , 

e 

IF  (j.eq.0)  THEN 
x  =  sqrt(-2.«alog(ranu(da))) 
y  «  2. *3. 1415926535«ranu(da) 
rolln  *  x*cos(y)«sigaa 
z  *  x*sin(y) 

E3E 

j  a  1-j 

rol  In  *  z'sigaa 
ENDIF 
END 
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SUBROUTINE  SEARC2  (t, firer, tgt,naray,eond,dt) 
e  ?  Se»rc2:  see  if  a  tank  detects  a  target  during  this  second, 
include  ’coMon.h’ 
integer  firer,  tgt,  cond 
c 

if  (trace)  print  «,’>seare2’ 
teep  »  rg/rgincr 
indx  »  int(teap) 

IF  (indx  .It.  1)  THEN 
tlo  *  1.0 

thi  x  tbar(naray,cond,l) 

BLSEIF  (indx  .It.  8)  THEN 
tlo  >  tbar(na ray, cond, indx) 
thi  x  tbar(naray,cond, indx»l) 

ELSE 

tlo  s  tbar(naray,cond,8) 
thi  x  0.0 
ENDIF 

frac  x  teap-aint(teap) 
pdetct  x  tlo  *  fraex(thi-tlo) 

IF  (ranu(0.0) .gt. pdetct)  THEN 

c  The  firer  doesn’t  detect  the  target  in  the  next  second, 

repeat  x  .true, 
dt  x  1.0 
E.SE 

c  This  firer  detects  the  target  in  this  second. 

call  skedul (t*ranu(0.0), firer, ’detect’,  tgt) 

ENDIF 

if  (trace)  print  «,’<searc2’ 

END 
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SUBROUTINE  SEARCH  (t} 

c  3  Search:  see  if  any  targets  are  detected  ini  the  next  second, 
include  'coaaon.h* 
logical  ignore 

coeeon  /cserch/  il, in, jl,jn,v, rgtbl (NN,NN) , ignors(NN) ,rgvs(NN) , 

1  yaax (NN) , i a  ray , j  arey , ndet i , ndet j , n i , n j 

save  /cserch/ 
rss(x,y)  *  sqrt(x«x»y»y) 
e 

if  (trace)  print  *,’>search' 
repeat  *  .false, 
c  Update  status  of  tanks. 

c  (Next  line  shud  eventually  be  updated  in  daaage.f,  I tk ill.) 

DO  20  i=l,nblu»nred 

ignore(i)  =  ignore(i) .or. I ife(i) .ge. IKILL 
IF  (.not. ignore(i)  .and.  aotion(i) .ne.STATNY)  THEN 
cal  I  path ( i ,t,aotion(i) ,0.0, da, da.de, da) 

DO  10  j*jl,jn 

rgtbl(i.j)  =  rss(x0(i)-x0(j)  ,y0(i)-y0(j)) 
rgtbl (j , i)  =  rgtbl(i.j) 

10  CONTINUE 

ENDIF 

20  CONTINUE 

c 

DO  40  i*i 1,  in 

c  Loop  thru  Southern  tanks. 

IF  (.not. ignore(i))  THEN 

c  Consider  tank  i  (It  is  alive  and  can  detect  or  be  detected.) 

c  Change  aade  March  20,  1989  by  H.L.Reed  to  alio*  the  individual  condition 
c  of  each  target  tank  to  be  used  to  define  the  probability  of  acquisition 
icond  *  2 

if  (not ion (i). ne.STATNY)  icond  *  3 
if  (knceal (i) .eq.HD)  icond  *  1 
DO  30  jsjl.jn 
IF  (.not. ignore (j))  THEN 

c  Consider  tank  j  (Also  alive  and  can  detect  or  be  detected.) 

jcond  «  2 

if (aot ion (j) .ne.STATNY)  jcond  =  3 

if  (knceal  (j) .eq.HD)  jcond  =  1 

rgi  *  rgv is (jcond, i) 

rgj  *  rgvis(icond, j) 

rgaax  *  aaaxlfrgi ,rgj) 

rg  *  rgtbl  (i,j) 

IF  (rg.lt. rgaax)  THEN 

c  At  least  one  is  in  detection  rg  of  the  other. 

IF  (los(i.j))  THEN 

c  Line-of -sight  exists  between  thes. 

c  Treat  Southern  tank  as  searcher 

IF  (rg.lt. rgi  .and.  ,not.seen(i , j)  .and. 

1  ndet(i) . It.ndeti)  THEN 

cal  I  searc2(t, i , j , iaray, jcond, dt) 

&SE 

repeat  *  .true. 

ENDIF 
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e  Treat  Northern  a*  searcher 

IF  (rg.lt.rgj  .and.  ,not.seen(j, i)  .and. 

1  ndet (j ) . I t . ndet j )  THEN 

call  searc2(t, j , i , jarny , icond.dt) 

ELSE 

repeat  *  .true. 

ENDIF 

ajsE 

repeat  *  .true. 

ENDIF 

ELSEIF  (y0(i)  .gt.yeax(i))  THEN 
c  Neither  are  in  detection  rg  of  the  other, 

call  cancel (i .’all  '.ALL) 
ignore(i)  *  .true. 

if  (keyd(l).ge.2)  print  Cancel  all  events  for', 
GOTO  41 
ENDIF 
ENDIF 

30  CONTINUE 

ENDIF 

40  CONTINUE 

if  (repeat)  call  skedul (t«l. 0,0, 'search',  NILL) 

if  (trace)  print  •, '(search1 

END 

c  V7.1 
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SUBROUTINE  SB-ECS  (t,firer,dt) 
include  'coeaon.h' 
logical  loaded 
integer  firer,  arayf 

1  foraat  (f 8 . 2 , lx , »4 , i 3 , "  doea  not  select;  selecting  already.’) 

2  forest  (f 8 . 2 , lx,a4, i 3 , 1  does  not  select;  channels  full.1) 

3  foraat  (f8.2,lx,a4, 13, 1  does  not  select;  pod  eapty.*) 

4  foraat  (f8.2,lx,a4, i3, ’  begins  selection.') 

c 

if  (trace)  print  *,’>se lees’ 
arayf  *  aray (firer) 

if  (kind. eg. 4)  leaded  *  nchan(f irer) .ge.nehans(arayf) 
if  (kind.ne.4)  loaded  *  nrtgt(f i rer) .ne.0 
IF  (busy(firer)  .or.  eapty(firer)  .or.  loaded)  THEN 
c  Wait  cause  busy  selecting,  pod  eapty,  or  channels  full. 

IF  (keyd(l) .ge.2)  THEN 
IF  (busy (firer))  THEN 
print  1,  t,  eolor(arayf) ,  firer 
B^EIF  (loaded)  THEN 
print  2,  t,  color(arayf) ,  firer 
B-SEIF  (eapty (firer))  THEN 
print  3,  t,  color(arayf) ,  firer 
ENDIF 
ENDIF 
ELSE 

c  Start  selection:  none  in  progress  and  a  channel  is  free, 

busy (firer)  »  .true. 

if  (kind. eg. 4)  nchan(firer)  *  nchan(f i rer)«l 
call  skedu I  (t*dt, firer, ’select’,  NULL) 
if  (keyd(l).ge.2)  print  4,  t,  color(arayf),  firer 
ENDIF 

if  (trace)  print  «,’<se lees’ 

END 
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SUBROUTINE  SEECT  (t,  firer) 

c  6  Select:  gunner  chooses  aost  dangerous  target  he  sees, 
include  ’coeaon.h’ 
character*4  colort 
logical  tgt  fls,  f  alive,  can  go 
integer  firer,  tgt,  priorn,  areyf 

1  forsat(f8.2,lx,a4,  i 3 ,  ’  selects  \a4,i3,'  with  priority’ ,  i4, 

1  '  ftgts*' , i 2) 

2  forsat(f8.2,lx,a4, i3, 1  selects  ’,a4,'  -1', 

1  '  l  discards  ',a4,i3,  ’  ftgts=’,i2) 

3  forsat(f8.2,lx,a4, i 3 , '  selects’ ,8x, *-  (espty  target  set)’) 

4  foreatC  SEECT:  * , a4 ,  i 3 ,  *  selects  ’ , a4 ,  i 3 ,  ’  with  priority*,  i4) 
c 

if  (trace)  print  «,’>select’ 
areyf  *  arny (firer) 
kind  =  kindrd(areyf) 
f  alive  a  I ife(f i rer) . It . FKIUL 
IF  (f  alive)  THEN 

c  Firer  can  shoot,  so  have  hie  select, 

tgt  »  priorn(t,f irer, level) 

IF  (tgt. eq. NULL)  THEN 

c  Firer  has  no  targets  to  select  so  he  aoves  if  possible 

if  (keyd(l). ge.2)  print  3,  t,color(areyf) ,  firer 
busy(firer)  *  .false. 

if  (kind.eq.4)  nchan(firer)  *  nchan(f i rer) -1 
IF  (can  go (firer, t)  .and.  (kind. I e . 2  .or. 

1  kind.eq.S  .or.  nchan(f irer) .eq.8))  THEN 

call  cancel (firer, ’halt  ',  NILL) 
call  eaneel(f irer, 'aecel  ',  NILL) 
call  skedul (t, firer, 'acce I  *,NILL) 

ENDIF 

ESE 

c  Tgt  has  been  selected 

colort  >  color(arey(tgt)) 

IF  (tf ire(f irer, tgt). Ie.8.)  THEN 
c  Tgt  is  new;  replace  with  false  tgt  randoely. 

i  a  knceal (tgt)-l 
pf  a  ranu(B) 

tgt  fls  a  pf  .It.  pfalse(areyf , i) 

IF  (tgt  fls)  THEN 
seen(f irer, tgt)  a  .false, 
if  (keyd(l) .ge.2)  print  2,  t,  color(areyf) , 

1  firer,  colort,  colort,  tgt,  nchan(firer) 

tgt  a  FLS  TGT 

c  Restart  search  if  it  is  turned  off 

IF  (.not. repeat)  THEN 
repeat  a  .true. 

call  skedul (t, 8, ’search' .NULL) 

ENDIF 

ESE 

fot (firer, tgt)  ■  .true. 

if  (keyd(l) .ge.2)  print  1,  t,  color(areyf) , 

1  firer, colort, tgt, level ,nehan (firer) 

ENDIF 

ESE 

c  Firer  has  previously  serviced  this  target, 

fot (firer, tgt)  »  .true, 
if  (keyd(l) .ge.2)  print  1,  t,  eolor(areyf), 

1  firer, colort, tgt, level ,nchan(f irer) 

ENDIF 

call  engage  (t,  t,  firer,  tgt) 

ENDIF 

nrtgt(firer)  «  tgt 
ENDIF 

if  (trace)  print  *,'<select’ 

END 
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SUBROUTINE  SERCH0 

c  SerchB:  Find  useful  constants  for  search, 

include  ’coaaon.h’ 
logical  ignore 

coaeon  /cserch/  i 1 , i n , j 1 , j  n , v , rgtb I (NN , NN) , i gnore (NN) , rgvs (NN) , 
1  y Bax  (NN) , i a  ray , j  a  ray , ndet i , nde  t j , n i , n j 
save  /cserch/ 
integer  ncols(3) 
data  ncols  /2,3, X / 

if  (trace)  print  *,  ’Jserchl* 
c  Find  1st  and  last  in  Southern  k  Northern  forces. 

IF  (scene. eq.BATTAK)  THEN 
il  *  1 
in  3  nblu 
jl  =  nblu»l 
jn  =  nblu«nred 
v  *  speed (BLU) 

ELSE 

il  x  nblu«l 
in  x  nblu«nred 

jl  =  1 

jn  x  nblu 
v  x  0.0 

if  (scene. eg. RATTAK)  v  *  speed (RED) 

ENOIF 

c  Find  actual  detection  ranges  for  targets, 

nb  x  ncols (scene) 

DO  20  ixl,nblu 
rgvs(i)  ■  rgvis(nb.i) 
ignore(i)  »  .false. 

20  CONTINUE 

nr  x  4-nb 

DO  25  ixnblu*l,nblu«nred 
rgvs(i)  »  rgvis(nr.i) 
ignore(i)  »  .false. 

25  CONTINUE 

iaray  x  aray(il) 
jaray  x  aray(jl) 
ndeti  x  ndets(iaray) 
ndetj  x  ndets (jaray) 
ni  x  2 
nj  *  2 

if  (scene. ne.UEETNG)  ni  •  3 
if  (scene. ne.UEETNG)  nj  *  1 
if  (trace)  print  «,  '<serch0’ 

END 
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SUBROUTINE  SERCH1 

c  Find  whether  8  when  search  should  be  started, 

include  'coaeon. h* 
logical  ignore 

coason  /eserch/  i 1 , i n , j 1 , jn , v , rgtb I (NN , NN) , i gnore (NN) , rgvs (NN) , 
1  yeax(NN) , iaray, jarey ,ndeti ,ndetj ,ni ,nj 
save  /cserch/ 

if  (trace)  print  «,  '>serchl' 

ca 1 1  serch0 

dtain=1.0el0 

c  Loop  thru  Southern  force  and  Northern  force. 

00  40  isil,in 
yain  *  l.el0 
yaax(i)  =  -I.elB 
00  30  j=jl,jn 
x  «  x0(j)-x0(i) 

1  *  y®(j) 

d  *  sqrt(x«*2  »  y*»2) 
rgtbl(i.j)  *  d 
rgtbl  (j.  i)  =  d 
r  «  aaaxl(rgvs(i) ,  rgvs(j)) 

IF  (r  .gt.  d)  THEN 

c  At  least  one  is  in  detection  range  at  tiae  zero, 

yain  «  aainl(yain,0.0) 
yaax(i)  *  aaaxl(yeax(i) ,y*sqrt(r**2-x*»2)) 

B.SE 

c  Neither  is  in  detection  range  at  tiae  zero. 

IF  (v  .gt.  0.0  .and.  abs(x)  .It.  r)  THEN 
c  At  least  one  sill  enter. 

q  ■  sqrt(r»»2  -  x*«2) 

yain  *  aaaxl(0.0,aainl(yain,  y  •  q)) 

yaax(i)  *  aaaxl(yaax(i) ,  y  ♦  q) 

ENDIF 

ENDIF 

30  CONTINUE 

c  See  if  tank  i  should  be  ignored  and  update  start  tiae  (dtain) 

dt  =  1.0el0 

if  (yain. eq. 0.0)  dt  =  0.0 

if  (v.ne.0.0)  dt  s  yain  /v 

ignore(i)  *  dt.gt.taax  .or.  yaax(i) . le.0.0 

if  (ignore(i))  call  cancel  (i , '  a  1 1  ’.ALL) 

IF  (keyd(l) .ge.2)  THEN 

if  (ignore(i))  print  «,i,'  Never  enters  detect  rg  of  foe. 
if  (.not. ignore(i))  print  *,  i, 

I  '  Enters  detect  rg  after  traveiing’.yain, ’  aetres.’ 

ENDIF 

dtain  ■  aainl(dtain,dt) 

40  CONTINUE 

if  (dtain. It. taax)  call  skedul (dtain, ALL, 'search' , ALL) 
repeat  *  dtain. It. taax 
if  (trace)  print  «,  '<serchl' 

END 
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SUBROUTINE  SKBXJL  (t, I, act, it) 

Schedule:  Schedule  an  event  for  later  execution, 
include  ’clock. h' 
charactered  act 

foraat(9x, ’skedul  ',i3,’  ’,a0,i3,’  at  tiee',f8.2) 

if  (prflag}  print  1,  I,  act,  it,  t 
IF  (nxidle.eq.0)  THEN 
If  storage  a  1 1  used  stop 

print  e,'  Storage  overloaded  with  too  aany  events.’ 
STOP 
ELSE 

Store  the  event 
Cut  storage  unit  froa  eapties 
n  a  nx idle 

nxidle  *  next(nxidle) 

Then  find  where  to  insert  this  event  in  the  event  list. 
IF  (nxevnt. I e.0)  THEN 
New  event  is  only  event 
next(n)  »  0 
nxevnt  3  n 
B-SE 

Then  find  where  to  insert  it. 

Point  to  first  2  events 
I  3  nxevnt 
a  3  next(l) 

Find  where  to  insert  thee 
IF  (t.ge.when(l))  THEN 
See  if  between  2  scheduled  events. 

Loop  till  found. 

IF  (a.no.l  .and.  t.ge.when(a))  THEN 
I  >  ■ 

■  3  next(e) 

GOTO  20 
B-SE 

Splice  new  event  into  list 
next(n)  3  ■ 

next(l)  ■  n  * 

ENDIF 

B.SE 

Place  new  event  as  aost  iaainent 
next(n)  3  nxevnt 
nxevnt  *  n 
ENDIF 
ENDIF 

Finally  store  event  info 
when(n)  3  t 
what(n)  3  act 
who(n)  3  I 
whoe(n)  3  it 
ENDIF 
ENO 
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SUBROUTINE  SLOW  UP  (t,  f irer) 
c  8  Slow  up:  siaulate  tank  starting  to  slow  down, 
include  'coaaon.h' 
integer  firer 

1  foraat  (f8 . 2, lx , a4 , i 3 ,  ’  continues  to  slow  up.’) 

2  foraat  (f8.2,lx,a4,  i3,  ‘  would  slow  up  if  it  seren"tl, 

1  ’  already  stopped.') 

3  foraat  (fS.2,lx,a4,  i3, '  brakes’  ,llx,'  (was  accelerating)') 

4  foraat  (fB.2,lxla4l i3, '  brakes’.llx, ’ (was  cruising)’) 

c 

if  (trace)  print  *,’>slowup' 
kind  nv  *  aotion(f irer) 
naray  =  aray(firer) 

IF  (kind  av.eq. SLOVNG)  THEN 
c  Previous  aotion  was  slowing 

i f (keyd (l).ge.2)print  1,  t,  color(naray) ,  firer 
ELSE  IF  (kind  av.eq.STATNY)  THEN 
c  Previous  notion  was  stationary 

if  (keyd (l).ge. 2) print  2,  t,  color(naray),  firer 
ELSE  IF  (kind  av.eq.ACCH.G)  THEN 
c  Previous  notion  was  accelerating 

if  (keyd (1) .ge. 2) print  3,  t,  color(naray) ,  firer 

call  path  (f irer, t,aotion(f irer), 0.0,x,y,vx,vy) 

dt  =  vy/decel (naray) 

aotion (firer)  *  SLOVNG 

call  skedul (t«dt,f irer, 'halt  NULL) 

ELSE  IF  (kind  av.eq.UAXVL)  THEN 
c  Previous  notion  was  cruising  at  aax  vel 

if  (keyd (1) .ge. 2) print  4,  t,  eolor(naray) ,  firer 
call  path  (firer,t,aotion(f irer), 1.0, x,y,vx,vy) 
c  schedule  halt  tiae 

dt  *  vy/decel (naray) 
call  skedul(t«dt, firer, 'halt  ’,  NULL) 
not ion  (firer)  a  SLOVNG 
ENDIF 

if  (trace)  print  «,’<slowup‘ 

END 


c  V7.2 
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c  e 


i 


c 


c 


70 

80 


c  V7.1 


SUBROUTINE  SUOKE 

Seoka:  Find  path  lengths  where  attack* r  ia  hidden  by  saoke. 
include  1 coaeon. h* 

data  ptbl  /0. , .05, .1, .16, .2, .’a, .3, .35, .4, .45, .5, .SS, .8, .65, .7, 
.75, .8.. 85, .9, .95, 1.0/ 
data  rtbl  /0. ,  1000., 2000., 3000., 4000./ 

if  (trace)print  », ’Dseoke’ 

DO  80  nb*l,nblu 
DO  70  nrEnblu«l,nblu*nred 

Find  first  tiae  window  for  LOS  between  tanks  nb,  nr. 
p=ranu(da) 
r=rg  0 

if  (kview(RS)) .eq. ‘V’)  dt=tdintp(ptbl,rtbl,toutvl,p,r,21,S) 
if  (kview(RED) .eq. ’I')  dt=tdintp(ptbl,rtbl,toutil,p,r,21,5) 
call  skedu I  (dt.nb, ’appear', nr) 

if  (kview(BLU) .eq. ’V’)  dtstdintp(ptbl,rtbl,toutvl,p,r,21,5) 
if  (kview(BLU)  .eq. ’I’)  dt=tdintp(ptbl,rtbl,toutil,p,r,21,5) 
call  skadul (dt, nr, ’appear', nb) 

CONTINUE 

CONTINUE 

i  f  (trace) pr i nt  * , 1 <snoke ’ 

END 
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FUNCTION  TDINTP(xla,  x2a,  y,  *1,  *2,  ixla,  ix2i  ) 
c  TOINP:  Interpolates  in  a  two  diiensional  utrix. 

integer  ill*,  ix2a 

real  y(ixla, ix2a) ,  xla(ixla),  x2a(ix2a) 
integer  j,  k 

real  yl,  y2,  y3,  y4,  t,  u 
c 

j  *  INDEXX(  xla,  ixla,  xl  ) 
k  *  INDEXX(  x2a,  ix2a,  x2  ) 

IF(k.eq.0)  THEN 

PRINT-,,  ’TDINTP:  p, r, j ,k=’ ,xl,x2, j ,k 
print  «,  xla,  x2a, ixla, ix2a 
ENOIF 

yi  =  y(J.k) 
y2  ■  y(i*V.k) 

y3  =  y(j-l,k*l) 
y «  -  y(j,l<*l) 

t  *  (xl-xla(j))/(xla(j-l)-xla(j)) 
u  *  (x2-x2a(k))/(x2a(k»l)-x2a(k)) 
c 

TOINTP  =  (l-t)*(l-u)»yl  ♦  t*(l-u)«y2  ♦  t»u»y3  ♦  (l-t)«u»y4 
END 

c  V7.2 
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SUBROUTINE  TERAIN  (if irst, last) 

c  0  Uask  at:  find  path  lengths  where  attacker  is  sasked  by  terrain 
include  'coaion.h1 

coaaon  /terane/  d  (40) ,  so  I d (20) ,  yold(20),  dist(20),  i seg (20) 

1  forsat  ('  visible  for* ,fS.0, then  hidden  for* ,fS.0, 

c 

if  (trace)  print  *1,>terain’ 

c  Find  segsent  length  at  start  of  each  engagement. 

00  20  i=l,39,2 

c  Hunfeld  terrain  constants 

f  *  -alog(ranu(0.0)) 
d(i)  »  300. *f **1.2 
f  =>  -alog(ranu(0.0)) 
c  d(i-l)  *  7S0.*f**2.0 

d(i-l)  *  100. *f 

if  (keyd(l) .ge.2)  print  1,  d(i),  d(i«l) 

20  CONTINUE 

c  Initialize  data  for  each  tank 

DO  30  i=ifirst, last 

call  path  (i ,0. ,aotion(i) ,0.0.x, y,vx,vy) 

xold(i)  =  x 

yold(i)  =  y 

dist(i)  =  d(l) 

iseg(i)  *  1 

call  skedul  (0. ,  i ,  'vanish' .NULL) 

30  CONTINUE 

if  (trace)  print  *,'<terain' 

c  NOTES: 

cl.  If  eeeting  engagement  skip  this  routine  entirely 
c  2.  If  side  not  moving  skip  at  least  part 
END 

c  V7.1 


FUNCTION  VABS  (») 

Vabs:  find  abslute  value  of  a  vector  (magnitude), 
dieension  a(3) 

vabs  *  sqrt(a(l)«»2  ♦  a(2)«»2  •  a(3)«»2) 

END 
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SUBROUTINE  VANISH (t, tgt, f irer) 

c  0  Vanish:  if  tgt  vanishes  treat,  otherwise  reschedule  vanish 
include  'coeaon.h' 
integer  tgt,  firer 

coeeon  /terane/  d (40) ,  xold(20),  yold(20),  dist(20),  iseg(20) 
!•**(*,  y)=sqrt(x«x*y«y) 
c 

if  (trace)  print  *,‘>vanishl 
narey  =  aray(tgt) 

IF  (invisb.eq.l)  THEN 

if  (speed(naray) .  le.0.)print  *, 'VANISH:  naray.speedz1 , naray, 
1  speed (naray) 

IF  (speed(naray) . le.0.)  STOP 
cal  I  path (tgt ,t,aot ion (tgt) ,0.0,x,y,vx,vy) 
c  Terrain  causes  intervisibility 

travel  =  rss(x-xold(tgt) ,  y-yold(tgt)) 

IF  (travel .gt.dist(tgt))  THEN 
c  Tgt  is  now  Basked  by  terrain 

xold(tgt)  =  x 
yo I d  (tgt)  3  y 
iseg(tgt)  *  iseg(tgt)«l 
if  (iseg(tgt) .gt.40)  iseg(tgt)=iseg(tgt)-40 
dist(tgt)  *  d(iseg(tgt)) 
call  vanter(t, tgt, firer) 
dt  *  dist (tgt) /speed (naray)  ♦  0.01 
call  skedul  (t«dt, tgt, 'appear* , NULL) 

ELSE  IF  (I ife(tgt) .eq. ALIVE)  THEN 
c  Not  yet  aasked  by  terrain,  so  reschedule 

dt  *  (dist (tgt)  -  travel)  /  speed (naray)  *0.01 
call  skedul  (t*dt, tgt, 'vanish* .NULL) 

ENDIF 

H.SE 

c  Tgt  is  now  aasked  by  saoke 

call  vansak(t, tgt, firer) 

ENDIF 

if  (trace)  print  «,’<vanish’ 

END 


c  V7.1 
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SUBROUTINE  VANSUK(t, tgt, firer) 
c  0  Vansak:  Treat  tgt  that  vanished  behind  saoke. 
include  ’eoaaon.h’ 
integer  tgt,  firer,  arayf,  arayt 

1  foraat(fB.2,lx,a4,  i 3 , '  LOS  to  * , a4 , 53,  *  broken  by  seoke.’) 

c 

if  (trace)  print  •1’>vansek‘ 
arayt  ■  a ray (tgt) 
arayf  «  3-arayt 

if  (keyd(l).ge.2)  print  1,  t,  color(arayf) ,  firer, 

1  color (arayt) , tgt 

c  Cancel  I ine-of -sight  between  tgt  and  firer. 

los(f irer, tgt)  =  .false. 

if  (seen (firer, tgt))  ndet(firer)  »  ndet(f irer)-l 
seen  (firer, tgt)  *  .false, 
tf ire  (firer, tgt)  =  0.0 
CHANCED  2  Oct  88.  Next  line  is  new. 

if  (busy(f irer) .and. nrtgt (firer) .eq. tgt)  busy(f irer)=. false, 
c  Abort  firer  aissile  on  tgt. 

IF  (aot(firer,tgt).or.fot(f irer, tgt))  THEN 
cal  I  diseng(t, firer, tgt, .true. ,.true.) 
if  (aot (firer, tgt))  call  abort (t,f i rer,tgt) 

ENDIF 

c  Accelerate  tgt  that  was  halting  to  fire. 

IF  (aot ion (tgt) .eq.SLOING  .and.  I ife(tgt) .eq.l  .and. 

1  fot(tgt, firer))  THEN 

call  skedul  (t,tgt, ’accel  ' , NULL) 
call  cancel  (tgt.'halt  1 , NULL) 

ENDIF 

if  (trace)  print  •,'<vansak’ 
c  NOTE:  shouldn't  halted  tgt  accelerate  too? 

END 


c  V7.3 
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SUBROUTINE  VANTER(t,tgt,f irer) 
c  0  Vanter:  Treat  tgt  that  vanished  behind  terrain, 
include  ’coeeon.h' 
integer  tgt,  firer 

1  foraat(f8.2, lx,a4,  i3,  *  vanishes’, 9x,'(xz',f8.1,'  yz' , f 8 . 1 ,  ’)  ’) 

c 

if  (trace)  print  «,'>vanter' 
naray  =  arey(tgt) 

if  (keyd(l).ge.2)  print  1,  t,  color(narey) ,  tgt, 

1  x0(tgt),  y0(tgt) 
knceal(tgt)  =  FD 
nrtgt(tgt)  =  0 

CHANGED  22  Sep  86.  Next  line  added. 
ndet(tgt)  *  0 

c  Cancel  all  I ines-of -sight  and  sightings  involving  tgt 

DO  20  isl,nblu*nred 
I  os (tgt, i)  *  .false. 
los(i,tgt)  *  .false, 
if  (seen(i ,tgt))  ndet(i)=ndet(i)-l 
seen(tgt.i)  ■  .false, 
seen (i, tgt)  =  .false, 
tf ire (tgt, i)  z  0.0 
tf  ire(i ,tgt)  =  0.0 
fot(tgt,i)  =  .false. 

c  Change  by  HLReed  1-12-90.  Seees  to  be  needed  to  clear  busy  flag 

c  GROUNDVARS  has  it  also. 

busy (tgt)  *  .false. 

20  CONTINUE 

c  Abort  outgoing  a issi I es 

call  abort(t,tgt, ALL) 
nchan(tgt)  z  • 

c  Abort  incoaing  rounds  k  disengage  tanks  firing  at  tgt 

if irstzl 

if  (naray.eq.l)  ifirst  z  nblu-1 
c  kind  z  kindrd(3-naray) 

call  newtgt(t, if irst.tgt) 
call  cancel  (tgt,’fire  ’ .NULL) 

CHANGED  27  Jun  85.  Added  next  line. 

call  cancel  (tgt, 'select' , NULL) 
c  Accelerate  tgt  that  *as  halting  to  fire. 

IF  (aot i on (tgt) . eq . SLOWNG  .and.  I ife(tgt) .eq.l)  THEN 
call  skedul  (t,tgt, 'accel  ’.NULL) 
call  cancel  (tgt, 'halt  '.NULL) 

ENDIF 

if  (trace)  print  «,  '(vanter* 

END 

c 

c  V7.1 
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SUBROUTINE  WAVES (n  scans) 
c  7  lavas:  loop  thru  vivas  of  rad  tanks, 
includa  ’coason.h' 


di sans  ion  i state (8) 

2 

forsat  (8 i 10) 
forsat(’  Total *, i4, i8, 

,  2 i 3 , 17x , i9) 

S 

iS,4fS.l 

7 

forsat(f0.2) 

CHANGED 

31  Mar  86.  Next  foraat  added. 

8 

1 

foraat(*  Shots  by: 

•  Fired 

.256,/ , * 

Hue  Red’,/, 
Wasted 

'<2i0,/i 

’.2.6,/, 

2 

1  Aborted 

,216,/,* 

False  tgts 

3 

'  Hidden  tgts 

,2i6./.* 

lapacting 

*.2i6 ,/, 

4 

'  Uisses 

,2i0./i ’ 

Hits 

’i2i0,/, 

5 

'  Duds 

,218,/,' 

No  daaage 

* .  2  i  6 ,  / , 

6 

•  U-kill  only 

,2.6,/, * 

F-kill  only 

* 1 2i6 , / , 

7 

’  lUF-kill  only 

i2i8,/, ’ 

K-kill 

‘,216,/) 

if  (traca)  print  *,'>Iavas’ 

IF  (nreps*nblu.gt.3000  .and.  nsavas.gt.l)  THEN 
print  «, 'WAVES:  Too  sany  raps  or  blues. * ,nreps,nblu 
ELSE 

c  Initial izo  scenario  statistics 

CHANGED  31  Uar  80.  Next  4  linos  added. 

DO  S  i*l,20 
kshot(l, i)  >  0 
kshot(2, i )  >  0 
S  CONTINUE 

nwavo  >  I 
nreps3  >  • 

loaso2  >  I 

noaao2  >  I 

scans  »  n scans 
DO  22  i-1,0 
state (i)  *  I. 

22  CONTINUE 

CHANGED  1  Apr  80.  Nsxt  tins  changed. 

c  nrg  »  rg0/500 

nrg  >  rg0/irginc 
DO  25  i >1,3000 
nosed (ij  >  0 
25  CONTINUE 

nsurv  >  nrepsvnblu 
c  call  hsadrl 

c  Loop  thru  up  to  10  vavss  of  red  tanks 

30  CONTINUE 

nwavs  >  n«ave*l 
ca 1 1  nxaavs 

IF  (nsurv.gs.nblu  .and.  naavs. It.nvaves)  GOTO  30 
c  Calculate  sueeary  statistics 

rpbs  ■  statc(S)/(nreps*nblu) 

IF  (nwavs. gt.l)  THEN 
DO  60  i«l,4 

istatc(i)  *  0.5  *  100»stitc(i)/nreps3 
state  (i  *4)  >  statc(i*4)/nreps3 
50  CONTINUE 

print  6,  nreps3, 

1  (istatc(i) , i>l,2) , (statc(i) , i>5,8) , 

2  loaso2,  noawo2,  i rands 
ENDIF 

CHANGED  31  Uar  88.  Next  line  added 

kshot(l,2)  ■  kshot(l,3)-kshot(l,4)»kshot(l,S) 
kshot(2,2)  >  kshot (2,3) »kshot (2,4) «kshot (2,5) 

ENDIF 

if  (trace)  print  •,'<wavas' 

END 
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SUBROUTINE  XFER(t,il,j) 

c  Change  HLReed.  Added  Xfer  for  the  transfer  a  detection  of  target  j 

c  to  all  vehicles  on  side  starting  with  il 

include  'coeaon.h1 
i2  *  nblu 

if  (il  .ne.  1)  i2  *  nblu*  nred 
DO  18  i  «  il,i2 
IF  (life(i)  .It.  FXILL  .and. 

1  ndet(i)  .It.  ndets(arey(i))  .and. 

2  los(i,j)  .and. 

3  .not. seen (i , j)  )  THEN 
seen(i,j)  *  .true. 
ndet(i)  ■  ndet(i)  *  1 

thuean  *  1.6  *  exp(rol In (0.5)) 
call  selecs(t, i .thuean) 

ENDIF 

10  CONTINUE 
END 


APPENDIX  B 


BASIC  PROGRAM  FOR  VULNERABILITY  TABLE 


Intentionally  left  blank. 


Appendix  B 


BASIC  prograa  to  produce  vulnerability  table  froe  BRL  vulnerability 
data  and  weapon  accuracy. 

100  DIU  PKH(8,2,11,4,8),  SIGMAX(3,3),  SIGMAY(3,3),  PK(3,2,4,8) ,  CARD(8) 

10S  DIM  PKILL (3 , 2 , 4 , 8) 

108  OIU  PH(8, 3, 2 , 6)  1 (range, case, exposure, angle) 

110  OPEN  'v6808.doc'  FOR  INPUT  AS  ft  'typical  naee  of  file  with  BRL  vui  data 

115  OPEN  'v4_basea.doc'  FOR  OUTPUT  AS  |3 

118  NS1C  *  0~  ’nsig  ■  0  for  base,  nsig  *  1,2,3  for  var  1,2,3 

120  INPUT  f2,  RRX,  EX,  OISP,  TYPE,  PKH(RRX/500,EX,DISP,TYPE,0) 

130  FOR  N  =  1  TO  6 

140  INPUT  #2,  PKH(RRX/500, EX, OISP, TYPE, N) 

150  NEXT 

160  INPUT  |2,  DUMMY 

170  IF  -(1*EDF(2))  GOTO  120 

171  ’  the  following  set  of  data  are  horizontal  (sigaax)  and  vertical 

172  ’  (sigeay)  dispersions  in  ails.  The  first  index  is  the  variation 

173  ’  as  picked  by  nsig,  the  second  is  the  case  where 

175  ’  case-1  is  stat  /  stat;  case=2  is  aoving  firer;  case=3  is  aoving  tgt 

180  SIGMAX(0, 1)  =  .52  :  SIGUAX(0,2)  =2.59  :  SIGMAX(0,3)  =2.54 

181  SIGMAX(l.l)  =  .41  :  SIGMAX(1,2)  =  .63  :  SIGUAX(1,3)  =1.07 

182  SIGMAX(2, 1)  =  .40  :  SICMAX(2,2)  =  .45  :  SIGUAX(2,3)  =  .88 

183  SIGMAX(3, 1)  =  .21  :  SIGMAX(3,2)  =  .24  :  SIGUAX(3,3)  *  .44 

190  SICUAY(0, 1)  =  .52  :  SIGUAY(0,2)  =2.28  :  SIGMAY(0,3)  =1.57 

191  SIGMAY(1,1)  =  .41  ;  SIGUAY(1,2)  *  .63  :  SIGUAY(1,3)  =1.07 

192  SIGUAY(2,1)  =  .40  :  SIGMAY(2,2)  =  .45  :  SIGUAY(2,3)  «  .88 

193  SIGMAY(3, 1)  =  .21  :  SIGUAY(3,2)  =  .24  :  SIGMAY(3,3)  =  .44 

195  VK  =  1!  1  a  factor  that  can  be  used  to  look  at  reduced  lethality 

200  HT  =  .375  '  half  height  of  turret 

210  WT1  =1.175  '  half  width  of  turret 

220  LT1  =1.475  'half  length  of  turret 

230  HH  =  1.5  '  height  of  hull 

240  INI  =  1.775  ’  half  width  of  hull 

250  LH1  =  3.375  '  half  length  of  hull 

255  '  define  the  values  of  the  cardioid  distribution  for  30  deg  increaents 

260  CARD(0)  =  .1657 

270  CARO (6)  »  .001 

280  FOR  N»*  1  TO  5 

290  THETA  =  .5236  «  N 

300  CARD(N)  *  .16687  «  .16477  *  COS (THETA) 

310  NEXT 

320  FOR  CASE  *  1  TO  3 

330  FOR  EX  =  1  TO  2  'hull  defilade  or  fully  exposed 

340  FOR  RANGE  =  1  TO  8 

350  SIGMAX  =  SIGMAX (NSIG, CASE)  «  RANGE  *  .5 

351  SICUAY  *  SIGMAY (NSIG, CASE)  •  RANGE  •  .5 

360  0  =  3. 28*SQR (SIGMAX  *  SIGMAY)  '  dispersion  in  feet  for  vul  table 

361  IF  0  >=  11!  THEN  0  *  10.999:  ELSE  IF  0  <  1!  THEN  0*1! 

365  OX  =  0  -  .5 

370  IF  OX  >  10  THEN  OX  =  10 
380  IF  DX  <  1  THEN  OX  =  1 

385  D  *  0  -  OX 

386  01  =  1!  -  0 

387  D1X  =  OX  ♦  1 

390  IF  EX  =  1  THEN  G0SU8  630:  B-SE  COSUB  720 

400  FOR  TYPE  *  0  TO  4  *  type  of  ki 1 1  where  TYPE  »  0  is  Ph 

410  TEMPPK  =  0! 

420  FOR  N  =  0  TO  8 

430  IF  TYPE  »  0  THEN  FACTOR  *  l!: 

ELSE  FACTOR  »  (  D1*PKH(RANGE,EX,DX  ,TYPE,N) 

.  0  «PKH (RANGE, EX , D1X, TYPE, N) )  •  VK 
440  TEMPPK  *  TEMPPK  .  CARO(N)  *  PH(RANGE,CASE,EX,N)  «  FACTOR 
450  NEXT  '  loop  on  N 
460  PK (CASE,  EX,  TYPE,  RANGE)  »  TEMPPK 
470  NEXT  '  loop  on  TYPE 

480  NEXT  '  loop  on  Range 

49 0  NB(T  '  loop  on  EX 

500  NEXT  ’  loop  on  case 

505  ’  convert  pk’s  to  the  intervals  called  p()  in  section  2.5. 

510  FOR  CX  *  1  TO  3 


B  -  2 


520  FOR  EX  *  1  TO  2 
530  FOR  RX  *  1  TO  8 
535  PKILL(CX,EX,0,RX)  -  PK(CX,EX,0,RX) 

540  PKILL(CX,EX,1,RX)  ■  PK(CX,EX,0,RX)-PK(CX,EX,4,RX) 

545  PKIU. (CX , E3S , 2, R*)  n  PK(CX,EX,0,RX)-PK(CX,EX,1,RX)-PK(CX,EX,2,RX).PK(CX,EX,3,RX) 

550  PKILL(CX,EX,3,RX).  =  PK(CX,EX,0,RX)-PK(CX,EX,2,RX) 

555  PKILL(CX,EX,4,RX)  -  PK(CX,EX,0,RX)-PK(CX,EX,3,RX) 

570  NEXT  1  loop  on  RX 

590  NEXT  1  loop  on  eX 

595  NEXT  ’  loop  on  cX 

596  G0SU8  1510 
600  END 

610  * 

620  '  Hull  defilade  hit  probability  given  sigaa  x,  sigia  y,  and  n  =  theta/30 
630  FOR  N  =  0  TO  8 
640  THETA  *  .5236  *  N 

650  WT  =  WTl  »  ABS (COS (THETA))  .  LT1  *  ABS (SIN (THETA)) 

660  GAUSS ARG  =  HT/SIGUAY:  GOSUB  860:  PHTEUP  *  2  <  GAUSS  -1 

670  GAUSSARG  =  IT/SIGUAX:  GOSUB  860:  PH(RANGE,CASE,EX,N)  =  (2  •  GAUSS  -1)»  PHTEUP 
680  NEXT 
690  RETURN 
700  ’ 

710  '  Exposed  target  hit  probability  given  sigaa  x,  sigaa  y,  and  n  «  theta/30 
720  FOR  N  =  0  TO  6 
730  THETA  *  .5236  •  N 

740  IT  =  WTl  •  ABS (COS (THETA))  -  LT1  «  A8S(SIN(THETA)) 

750  GAUSSARG  *  (.3  ♦  21  *  HT)/SIGUAY:  GOSUB  860:  PHTEUP  *  GAUSS 
760  GAUSSARG  =  ,3/SIGUAY:  GOSUB  860:  PHTEUP  *  PHTBIP  -  GAUSS 

770  GAUSSARG  =  WT/SICUAX:  GOSUB  860:  PH (RANGE, CASE, EX, N)  =  (2  »  GAUSS  -1)*  PHTEUP  ’  ph  turret 

780  GAUSSARG  *  ,3/SIGUAY:  GOSUB  860:  PHTBIP  *  GAUSS 

790  GAUSSARG  «  (.3  -  HH)/SICUAY:  GOSUB  860:  PHTBIP  -  PHTEUP  -  CAUSS 

800  WH  =  WH1  *  ABS (COS (THETA) )  ♦  LH1  *  ABS (SIN (THETA)) 

810  GAUSSARG  *  WH/SIGUAX:  GOSUB  860:  PH (RANGE, CASE, EX, N)  »  PHTEUP  .  (2  .  GAUSS  -1)  -  PH (RANGE, CASE, EX ,N) 
820  NEXT 
830  RETURN 
840  1 

850  '  Noraal  Distribution  Subroutine 
860  TEUP  *  ABS (GAUSSARG) 

870  GAUSS  .  .398942  •  EXP(-.S  •  GAUSSARG  «  GAUSSARG) 

880  IF  TEUP  >  4.6844  THEN  GOTO  940 
890  TEUP  *  1 ! / (1 !  ♦  . 2316419«TEUP) 

900  GAUSS  *  1!  -  GAUSS  •  TEUP  »  ((((1.33027  »  TB4P  -  1.821256)  »  TEUP 
♦  1.781478)  *  TEUP  -  .3565638)  .  TB4P  «■  .3193815) 

910  IF  GAUSSARG  <  0  THEN  GAUSS  =  1!  -  GAUSS 
920  RETURN 


930  ’Approxiaation  for  large  values  of  the  arguaent 
940  GAUSS  =  1!  -  GAUSS  .  (ll/TEUP  -  1!/TEUP*3  .  3!/TB4P*5) 

950  IF  GAUSSARG  <  0  THEN  GAUSS  *  1!  -  GAUSS 

960  RETURN 

1510  FOR  CX  *  1  TO  3 

1520  FOR  EX  x  1  TO  2 

1530  FOR  TX  *  0  TO  4 

1540  PRINT  #3,  USING  *f  »  f  f .f»f';CX,EX,TX,PKIH(CX,EX,TX,l); 


1550  FOR  RX  «  1  TO  8 
1560  PRINT  |3,  USING 
1570  NEXT 
1575  PRINT  |3,  •  • 
1580  NEXT 
1590  NEXT 
1595  NEXT 
1600  RETURN 


'  M|f;PKIU.(CX,EX,TX,RX); 
loop  on  RX 


loop  on 
loop  on 
loop  on 


tX 

eX 

el 


range  0  is  range  1 


No  of  No  of 

Copies  Organization  Copies  Organization 


1  Office  of  the  Secretary  of  Defense 
OUSD(A) 

Director,  Live  Fire  Testing 
ATTN:  James  F.  O’Bryon 
Washington,  DC  20301-3110 

2  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 

1  HQDA  (SARD-TR) 

WASH  DC  20310-0001 

1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

1  Commander 

US  Army  Laboratory  Command 
ATTN:  AMSLC-DL 
Adelphi,  MD  20783-1145 

2  Commander 

US  Army,  ARDEC 
ATTN:  SMCAR-IMI-I 
Picatinny  Arsenal,  NJ  07806-5000 

2  Commander 

US  Army,  ARDEC 
ATTN:  SMCAR-TDC 
Picatinny  Arsenal,  NJ  07806-5000 

1  Director 

Benet  Weapons  Laboratory 
US  Army,  ARDEC 
ATTN:  SMCAR-CCB-TL 
Watervliet,  NY  12189-4050 

1  Commander 

US  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  SMCAR-ESP-L 
Rock  Island,  IL  61299-5000 

1  Commander 

US  Army  Aviation  Systems  Command 
ATTN:  AMSAV-DACL 
4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120-1798 


1  Director 

US  Army  Aviation  Research 
and  Technology  Activity 
ATTN:  SAVRT-R  (Library) 

M/S  219-3 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 

1  Commander 

US  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R  (DOC) 
Redstone  Arsenal,  AL  35898-5010 

1  Commander 

US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-TSL  (Technical  Librarv) 
Warren,  MI  48397-5000 


1  Director 

US  Army  TRADOC  Analysis  Command 
ATTN:  ATAA-SL 

White  Sands  Missile  Range,  NM  88002-5502 

(a***-  <»iy)  1  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD  (Security  Mgr.) 

Fort  Benning,  GA  31905-5660 

(Undue,  only)  \  Commandant 

US  Army  infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905-5660 


1  Air  Force  Armament  Laboratory 
ATTN:  AFATL/DLODL 
Eglin  AFB,  FL  32542-5000 


Aberdeen  Proving  Ground 


2  Dir,  USAMS.aA 

ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 

1  Cdr,  USATECOM 

ATTN:  AMSTE-TD 

3  Cdr,  CRDEC,  AMCCOM 

ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-MSI 

1  Dir,  VLAMO 

ATTN:  AMSLC-VL-D 
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Commander 

US  Army  Laboratory  Command 
ATTN:  AMSLC-TP-PB,  I.  Bartky. 

F.  Ostovic 

Adelphi,  MD  20783-1145 

Commander 
US  Army,  ARDEC 
ATTN:  SMCAR-TD 
SMCAR-TDT 
SMCAR-SCF,  V.  Calgano 
E.  Del  Coco 
M.  Barbarisi 

SMCAR-CCS 

SMCAR-CCL-FW,  H.  Kahn 
Picatinny  Arsenal,  NJ  07806-5000 

Commander 

US  /umy  Wateivliet  Arsenal 
ATTN:  SMCWV-QAR,  Building  44 
W.  Jarrett 
Watervliet,  NY  12189 

Commander 

US  Army  Tank  Automotive  Command 
ATTN:  AMSTA-RY,  R.  Beck 
AMCPM-BFVS 
AMCPM-ABMS 
AMCPM-ABMS-SW 
AMCPM-AB  MS-SI, 

G.  Vander  Waerdcn 
Warren,  MI  48090 

Commandant 

US  Army  Infantry  School 

ATTN:  ATSH-TSM-FV 

ATSH-CD-MLS-M 
Fort  Benning,  GA  31095 


3  Commander 

US  Army  Armor  Center 
ATTN:  ATSB-CD 

ATZK-CD-MS,  Mr  Falkovich 
ATZK-AE-PD,  Mr.  Wells 
Fort  Knox,  KY  40121 

1  Commander 
TRAC  RPD 
ATTN:  ATRC-DCS 
Fort  Monroe,  VA  23651 

1  Director 

Development  Center 
ATTN:  MCDEC/D092, 

Firepower  Division 
Quantico,  VA  22134 

1  Director 

Benet  Weapons  Laboratory 
US  Army,  ARDEC 
ATTN:  SMCAR-CCB 
Watervliet,  NY  12189 

2  PM-TMAS 

ATTN:  DRCPM-TMA-PA,  K.  Rubin 
DRCPM-TMA,  F.  Steinberg 
Picatinny  Arsenal,  NJ  07806 

1  OUSD(A) 

ATTN:  Executive  Director,  DSB 
Room  3D1020,  Pentagon 
Washington,  DC  20301 

1  OUSD(A) 

ATTN:  DUSD(R&AT) 

Room  3E114,  Pentagon 
Washington,  DC  20301 

1  OUSD(A) 

ATTN:  DUSD(TWP) 

Room  3E1044,  Pentagon 
Washington,  DC  20301 


No.  of 
Copies 

2 

1 

1 

1 

1 

1 

I 

1 

1 


Organization 

ADUSD 

ATTN:  TWP/LW,  Room  3D1049/Viilu 
TWP/AW,  Room  3E1049/Loder 
Pentagon 

Washington,  DC  20301 

Assistant  Secretary  of  Defense  of  C1 
Room  3E172,  Pentagon 
Washington,  DC  20301 

OUSD(A) 

ATTN:  Director,  Program  Integration 
Room  3E1034/Christie,  Pentagon 
Washington,  DC  20301 

Director,  PA&E 
ATTN:  Director,  LFD 
Room  2B256,  Pentagon 
Washington,  DC  20301 

Director 

DARPA 

1400  Wilson  Blvd. 

Rosslyn,  VA  22209 

HQDA  (SAUS-OR,  Mr.  Hollis) 

WASH  DC  20310-0001 

HQDA  (SARD-ZD) 

WASH  DC  20310-0001 

Commander 

USAMC 

ATTN:  AMCCG 

5001  Eisenhower  Avenue 

Alexandria,  VA  22333-0001 

Commander 

AVSCOM 

ATTN:  AMCPM-AAH 
4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120-1798 

Commander 

CECOM 

ATTN:  AMSEL-RD-EW-D 
Fort  Monmouth,  NJ  07703-5303 


No.  of 

Copies  Organization 

1  Commander 
CNVEO 

ATTN:  AMSEL-RD-NV-T 
Fort  Bel  voir,  VA  22060-5166 

1  Commander 
ASL/LABCOM 
ATTN:  SLCAS-D 

White  Sands  Missile  Range,  NM  88002-5501 

1  Commander 
VAL/LABCOM 
ATTN:  SLCVA-D 

White  Sands  Missile  Range,  NM  88002-5513 

2  Commander 

US  Army  Missile  Command 
ATTN:  AMCPEO-FS 
AMSMI-RD-AS 

Redstone  Arsenal,  AL  35898-5001 
1  Director 

Survivability  Management  Office 
ATTN:  SLCSM-GS 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 

1  Commander 
USA  FSTC 
ATTN:  William  Rich 
220  Seventh  St.,  NE 
Charlottesville,  VA  22901-5396 

1  Commander 
TRACRPD 
ATTN:  ATCG 

Fort  Monroe,  V  A  23651-5000 

1  Commander 

USA  Air  Defense  Artillery  School 
ATTN:  ATSA-CD 
Fort  Bliss,  TX  79916 
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No.  of 
Copies 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


Organization 

Commander 

USA  Armor  Center  and  Fort  Knox 

ATTN:  ATSB-CD 

Fort  Knox,  KY  40121-5000 

Commander 

USA  Combined  Arms  Combat  Development 
Activity 

ATTN:  ATZL-CA 

Fort  Leavenworth,  KS  66027-5300 

Commander 

USAAVNC  and  Fort  Rucker 

ATTN:  ATZQ-TDS 

Fort  Rucker,  AL  363*7-5 163 

Commandant 
USA  Infantry  School 
ATTN:  ATSH-CD-MLS 
Fort  Benning,  GA  31905-5000 

Commander 

USA  Intelligence  Center  and  School 

ATTN:  ATSI-CD-TE 

Fort  Hauchuca,  AZ  85613-7000 

Director 
TRAC- WS  MR 

White  Sands  Missile  Range,  NM  88002-5502 

Assistant  Secretary  of  the  Navy 
Research,  Engineering  and  Systems 
Room  4E736,  Pentagon 
Washington,  DC  20350 

Commandant 

USMC 

ATTN:  CodeRD 
Washington,  DC  20380-0001 

Deputy  Commanding  General 
MCRDAC 
ATTN:  Code  D032 
Quantico,  VA  22134-5080 

Director 
USMC  OTEA 
Quantico,  VA  22134 


No.  of 

Copies  Organization 

1  Director 
DIA 

ATTN:  RTS-2,  Col  MacNicoU 
Pentagon 

Washington,  DC  20301-6111 

2  Director 
CIA 

ATTN:  Mr.  John  Ewen/Mr.  Robert  Gomez 
P.O.  Box  1925 
Washington,  DC  20505 

1  Joint  Electronic  Warfare  Center 
ATTN:  MAJ  Harry  Ladewig 
San  Antonio,  TX  78243-5000 

3  Los  Alamos  National  Laboratory 
ATTN:  CDT  (J.  Immele) 

AT  AC  (F.  Day,  P.  Howe) 

Los  Alamos,  NM  87545 

1  Lawrence  Livermore  National  Laboratory 
ATTN:  Dr.  Milton  Finger 
Livermore,  CA  94550 

1  Geometric  Solutions 

ATTN:  Mr.  Harry  L.  Reed,  Jr. 

100  Custus  Street 
Aberdeen,  MD  21001 

Aberdeen  Proving  Ground 

1  Dir,  USAMSAA 
ATTN:  Mr.  T.  Ruth 

1  Dir,  USAHEL 
ATTN:  J.  Waugh 

1  Cdr,  USACSTA 

ATTN:  STCS-CC-P,  P.  McCall 

1  Commander 

Project  Manager  Smoke/Obscurants 
ATTN:  AMCPEO-CNS-CT 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number  BRL-CR-638  _ Date  of  Report  APGPST  1990 - 

2.  Date  Report  Received _ _ _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ — — — — - - 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ _ _ - _ _ _ — 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.)  _ 


Name 


CURRENT 

ADDRESS 

Organization 

- 

Address 

City,  State,  Zip  Code 


7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


FOLD  HERE 


Department  of  the  Army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  2101.  -5066 

OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  No  0001,  APG,  MD 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 

% 


Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-9989 


FOLD  HERE 


